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FOREWORD 


The research described in this report was conducted by Mechanical 
Technology Incorporated under NASA contract NASw-1713. Mr. Henry B. 
Tryon of the Lewis Research Center Space Power Systems Division was the 
NASA Project Manager. 

The computer program used to calculate transverse vibration response 
of the rotor-bearing system was developed by Dr. Jorgen Luhd of the 
Technical University of Denmark, Copenhagen, in his capacity as consultant 
to MTI. The program is an extension of one of MTI’s basic rotor- response 
programs previously developed by Dr. Lund. The axial response computer 
program was developed by Mr. Paul Spencer of MTI. 

The metallurgical analysis reported in Appendix D was conducted by 
Mr. Charles Nolan of the Watervliet Arsenal, Watervliet, New York. 

The report was originally issued as Mechanical Technology Incorporated 
report MTI-70TR19. 
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SUMMARY 


Testing of a high-speed gas-bearing-supported rotor assembly has been performed 
to obtain an initial assessment of the effects of environmental shock and vi- 
bration on the mechanical performance of Brayton Cycle space-power turbomachin- 
ery. The rotor-bearing assembly consisted of a 10.5-pound, 39,000-rpm rotor, 
supported by two self-acting, pivoted-pad journal bearings and a self-acting, 
spiral-grooved thrust bearing. Sinusoidal vibration testing was conducted 
over a frequency range of 5 to 2000 Hz. Shock testing was performed with 10- 
millisecond half-sine pulses of approximately 20-g peak acceleration. The 
vibration and shock conditions were applied in three orthogonal directions: 
along the rotor axis (axial response), and transverse to the rotor axis in 
two directions (transverse response). Testing was conducted with and without 
shaft rotation, and with and without the rotor-bearing assembly mounted on 
vibration isolators. Dynamic displacements of the assembly were measured 
throughout the tests with capacitance probe equipment. 

Further objectives of the investigation were as follows: to identify the design 
factors which most strongly influence shock and vibration response of gas lub- 
ricated space-power turbomachinery, and to validate analytical methods for 
predicting dynamic response of such turbomachinery. To this end, an analytical 
study of the shock and vibration response of the high-speed rotor-bearing test 
assembly was performed*. For axial response calculations, a nonlinear three- 
degree-of-freedom analytical model was developed. The equations of motion 
were programmed for digital computer solution using a fourth-order Runge-Kutta 
numerical integration procedure. Solutions can be obtained for various types 
of excitations: sinusoidal vibration, half-sine or haversine shock pulses, or 
any arbitrary vibration or shock characteristic expressed in terms of displace- 
ment (or acceleration) and time coordinates. 

Analysis of transverse response was performed only for the case of sinusoidal 


* The analytical study was performed only for the condition of shaft rotation 
(i.e., for the condition of hydrodynamic gas-film lubrication). 



vibration excitation. A digital computer program was used to calculate steady- 
state vibration amplitudes of the rotor-bearing . system based on a finite- element 
model of both the rotor and the rotor casing. A total of eight, frequency- 
dependent, linearized coefficients were used to represent the stiffness and 
damping characteristics of each journal bearing. 

The significant results of this investigation are summarized below. From a 
quantitative standpoint, it must be remembered that the results apply to one 
specific rotor-bearing system. From a phenomenological standpoint, however, 
the system response characteristics are indicative of those which would be ex- 
hibited by similarly configured gas- lubricated rotor-bearing systems. 

Shock Tests 

1. Gas- lubricated rotor-bearing systems of the type tested can satisfac- 
torily survive at least a limited number of externally imposed 10- 

-k 

millisecond, 20-g shock impulses , both with and without shaft rota- 
tion, and with or without isolation mounts. For the system tested, 
momentary contacts between the rotating and stationary bearing sur- 
faces were observed under all test conditions. However, there was no 
surface damage or degradation of bearing performance as a result of 
these contacts 

Sinusoidal Vibration Tests — Shutdown Condition (Nonrotating Shaft) 

1. Under shutdown (nonrotating) conditions, gas-bearing- supported rotor 
systems of the type tested can satisfactorily survive, at least for 
a limited period of time, axial imposition of the specified sinusoidal 
vibration conditions, either with or without isolation mounts. 


A total of 73 shock pulses were accumulated during the test program. Four of 
the axial pulses were inadvertently conducted at 60-g peak acceleration. 

Based on previous work, it is known that the ability of gas bearings to suc- 
cessfully survive repeated high-speed contacts is strongly dependent upon 
the bearing materials used. The experimental results reported herein were 
all obtained using bearing parts whose mating surfaces were coated with 
plasma-sprayed chrome oxide. Based on considerable testing, chrome oxide is 
the optimum gas-bearing surfacing material thus far identified for operation 
at temperatures up to 600^F. 



2. For transversely imposed vibration under shutdown (nonrotating) con- 
ditions, vibration isolation may be required to survive the specified 
sinusoidal excitation conditions. During the nonrotating transverse 
vibration tests without isolators, the objective input excitation 
levels could not be achieved in the frequency range from 190 to 235 
Hz, this being the region of one of the critical speeds (resonant 
frequencies) of the nonrotating shaft. The limiting factor in this 
frequency range was deflection (overstressing) of the flexures used 
to support the individual journal bearing pads. With isolators in- 
stalled, excitation of the flexures was greatly reduced and the ob- 
jective input vibration levels could be safely imposed in the two 
transverse directions. 

3. A total of 27 minutes of vibration testing was accumulated under non- 
rotating conditions. No damage to the bearing surfaces nor degradation 
of bearing performance was detected as a result of the nonrotating 
tests • 

Sinusoidal Vibration Tests — Normal Operation (Rotating Shaft) 

1. Under normal rotating conditions, gas- lubricated rotor-bearing systems 
of the type tested may have to be vibration isolated to survive the 
specified sinusoidal vibration conditions. During the shaft rotation 
tests without isolation, rotor- to- bearing contacts became imminent, 
over broad regions of the frequency spectrum, at input vibration 
levels considerably below the test objectives. However, with isola- 
tors installed, the objective vibration input levels were achieved, 
without bearing contacts, over most of the frequency range. Only in 
the vicinity of the isolator resonant frequency (approximately 12 Hz), 
and in the vicinity of the rotor critical speeds (130, 166, and 225 
Hz), did bearing contact become imminent at less than the objective 
vibration input levels. 

2. Input vibration levels were carefully controlled during the rotating 
vibration tests to prevent the occurrence of bearing contacts. Con- 
sequently, it is not known whether' bearing performance would have 


been degraded under vibration-induced contact conditions. 


Calculated Shock and Vibration Response 

1. Accurate predictions of the axial shock and vibration response of gas- 
bearing machinery can be obtained using a three-degree- of freedom 
analytical model of the rotor-bearing system with nonlinear representa- 
tion of the load and damping characteristics of the thrust bearing gas 
film. An approximate, simplified method of representing the gas- film 
nonlinearities can be used which yields both accurate and economic 

(in terms of computer cost) solutions. The only drawback of the herein 
presented axial response analysis is the modeling of the vibration 
isolators in terms of linear stiffness and damping coefficients. Under 

shock conditions (i.e., large displacement conditions) the isolators 

•k 

are quite nonlinear and must be so modeled . 

2. Accurate predictions of the transverse resonant frequencies of the 
rotor-bearing test system were obtained using the finite-element 
analysis model of the rotor and rotor casing, together with the linear- 
ized journal bearing stiffness and damping coefficients. There was, 
however, considerable discrepancy between the calculated and measured 
amplitudes of vibration, and of the rotor- to-casing orbit shapes. Cal- 
culated amplitudes were considerably higher than the measured ampli- 
tudes. The discrepancy in calculated amplitudes may be due to "large 
amplitude" nonlinear stiffness and damping effects in the journal 
bearings which were not represented in the analysis. There was also a 
small amount of horizontal vibration (cross talk) in the vibration 
table which was neglected in the calculations . Finally, there may have 
been significant damping from the thrust bearing and from the shaft 
labyrinth seals. These sources of damping were also neglected in the 
analysis . 


The axial- response computer program has, in fact, been modified during Task 
3 of the present contract to permit nonlinear representation of the isolator 
deflection and damping characteristics. The results of additional Task 3 
calculations with the nonlinear isolator model will be reported in Part II 
of this Final Report. 
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3. The analysis results clearly show the beneficial effects of gas-bearing 
squeeze- film damping on minimizing the severity of shock-induced bear- 
ing contacts. Gas-bearing systems which must operate under shock 
conditions should be optimized to take maximum advantage of this 
effect. 

The above listed results indicate that 10-millisecond, 20-g shock pulses should 
not be a problem for gas- lubricated Brayton Cycle space-power turbomachinery. 
Sinusoidal vibration could be a problem, particularly if imposed during shaft 
rotation. The use of vibration isolators appears to be a solution to the vibra- 
tion problems during shutdown (nonrotating) conditions. During rotation, the 
use of isolators greatly reduces, but may not eliminate, the possibility of 
vibration-induced bearing contacts. Further testing will be required to 
identify if such contacts would cause bearing damage or degradation of bearing 
performance . 
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INTRODUCTION 


As America’s space capabilities grow, it becomes more and more apparent that a 
long-life, reliable space-power source is required for the 10-15 kilowatt range- 
One likely candidate for a space-power system of this range is the Brayton Cycle 
System- This consists of an inert gas in a closed loop which, after being heated 
in a recuperator and a heat source, drives one or more turbines turning an 
alternator and a compressor before being cooled again in the recuperator and 
radiator. To avoid the problems associated with an oi 1- lubricated bearing 
system, the rotating equipment can be supported by gas bearings utilizing the 
working gas as a lubricant. These bearings have the added advantage of long 
life, since bearing and shaft are separated during rotation by a film of gas. 

Several items of Brayton Cycle rotating equipment have been built during recent 
years. This equipment has been tested by the NASA Lewis Research Center and by 
the associated contractors under both design and off-design conditions with 
excellent results. However, an assessment of the effects of environmental shock 
and vibration was not included in these initial tests. Although some low- 
frequency vibration and shock testing for Naval applications has been success- 
fully conducted on a 62-pound, 8,000-rpm gas-bearing rotor [l], it was not felt 
that this experience was applicable to the dynamic environment associated with 
space-power applications. In particular, it was not known whether the latter 
environment might cause sufficient displacements of the shaft or bearing support 
assembly to result in contact between moving parts, fatigue, deformation, or 
fracture of parts with resulting failure of the machine. 

In order to gain insight into the effects of shock and vibration upon gas-bearing 
machinery for space- power applications, the NASA Lewis Research Center contracted 
with Mechanical Technology Incorporated to conduct a combined analytical and 
experimental investigation using a high-speed, gas-bearing machine representative 
of the type under consideration. This was an especially opportune time, since 
a turbocompressor gas-bearing simulator of the general size, weight, and speed 
required had recently been made available from another NASA contract pj. With 
minor modifications and some added instrumentation, the machine was well suited 
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for the experimental portion of this investigation. 


The analytical portion of the investigation was devoted to developing calcula- 
tion techniques and computer programs to the point where the effects of shock 
and vibration could be analytically predicted for a given rotor-bearing system. 
The experimental portion of the investigation consisted of subjecting the 
modified gas-bearing simulator to shock and vibration. This was done with the 
machine mounted in each of three mutually perpendicular positions, both with and 
without the shaft rotating, and with and without the use of isolators. The 
analytical results are compared with the measured results in this report, and 
the comparisons are used to validate the analytical procedures. 
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DESCRIPTION OF EQUIPMENT 


A schematic of the turbocompressor simulator mounted on the test stand is shown 
in Figure 1 . The rotor, shown in cross-hatching, consists of the following com- 
ponents: 

a) The drive turbine 

b) The turbine-end journal 

c) A center section (which would normally hold 
a six-stage, axial compressor) 

d) The compressor-end journal 

e) The thrust runner » 

Cross-sections AA and BB of Figure 1 show the configuration of the gas-lubrica- 
ted, compressor-end and turbine-end journal bearings respectively. Self- 
acting, pivoted-pad bearings are used. Each pad is individually supported by 
a mechanical flexure to permit the bearings to accommodate radial centrifugal 
growth of the journals and differential, radial, thermal expansions between 
the various bearing partSo Details of the journal-bearing designs are given 
in Reference 2. 

The gas -lubricated thrust-bearing assembly consists of a thrust runner, & 
reverse-thrust stator, a forward- thrust plate, and a support flexure and sup- 
port housing for the forward-thrust plate. The forward bearing is the primary 
thrust bearings the reverse bearing is used only during turbocompressor start- 
ip. The forward bearing is a self-acting, spiral-groove type. The reverse 
bearing is an externally pressurized type. Design details for the complete 
thrust -bearing assembly are likewise documented in Reference 2. 

Highly successful operation of the simulator gas bearings, up to 60,000 rpm, 
was demonstrated under static environmental conditions on a previous NASA 
contract [2], For the shock-and- vibration testing described herein, the 
simulator was operated between 38,000 and 39,000 rpm, this being close to 
the design speed of a Brayton-cycle , space-power system currently being 
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investigated by NASA [3]. 


Figure 2 is a schematic of the turbocompressor simulator oriented horizontally. 
Also shown in Figure 2 is the capacitance probe coding scheme, which applies to 
the probes shown in Figure 1 as well. 

Figure 3 is a photograph of the simulator and test stand mounted on the vibration 
shaker ready for testing. In this photograph, the simulator casing is rigidly 
mounted to the test stand. To* perform tests with the casing on isolators, the 
two mounting blocks shown in Figure 3 were removed, two additional legs were 
added to the test stand, and a shock-and-vibration isolator was installed be- 
tween the top of each leg and the simulator casing. The isolator mounting 
arrangement was symmetric about the rotor axis. 

Figure 4 shows the tripod mounting arrangement used to mount the simulator, in 
a horizontal attitude, on the vibration table (shaker). The rigid mounting block 
on the right side of the photo has a hidden counterpart away from the viewer. 
These mounting blocks were interchangeable with the isolators used when the sim- 
ulator was tested in the vertical as well as the horizontal attitude. 

Figure 5 is a photograph of the main components of the rotor-bearing assembly. 
Here we see the forward -thrust bearing in some detail. The self-lubricating 
action of this bearing may be visualized by imagining the face of the thrust 
runner being placed upon the thrust plate and rotating clockwise at high speed. 
The runner drags or pumps the ambient gas inwardly along the lighter-appearing 
spiral grooves causing a high-pressure region at the sealed end of the grooves. 
This pressure produces the load-carrying capacity required to support the steady- 
state axial loads (aerodynamic thrust and rotor weight, for example) as well as 
the dynamic loads resulting from externally-imposed axial shock and vibration. 

The simulator was designed to also support thrust loads with the rotor in a non- 
rotating condition. This is accomplished by supplying externally-pressurized 
gas to orifices in the thrust plate during startup and shutdown (hydrostatic 
operation). Figure 6 shows the size and location of these orifices, along with 
the spiral grooves which provide the hydrodynamic operation when the rotor 
achieves operating speed and the external pressure is withdrawn. A similar 
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arrangement of hydrostatic orifices was used in the reverse thrust bearing. 

Figure 7 shows the arrangement of a typical journal-bearing pad. Four such 
pads were situated around each of the two journal bearings. A capacitance- 
probe lead is shown along with the pivot assembly. The pivot itself is not 
visible in Figure 7, but is located between the convex surface of the pad and 
the pivot assembly. The two short beams protruding from the assembly consti- 
tute the flexure when the two lugs are bolted to the casing. 

Figure 8 shows the shock machine used in the experimental phase of the program. 
The simulator is shown mounted vertically on isolators. The operation of the 
shock machine is as follows: the table is raised pneumatically and then released 
to fall onto the elastomer pad shown resting on the anvil. The resulting im- 
pact produces a shock pulse that is roughly half-sine in shape. 

Figure 9 shows the vibration shaker with the simulator mounted horizontally on 
isolators. In this position it was necessary to add lead weights to the com- 
pressor end of the casing in order to properly load the isolators without 
nodification of the simulator casing. 

Instrumentation 


A permanent complement of capacitance probes was used to measure diaplacements 
of the rotating journals relative to the simulator casing, as well as film 
thickness in the journal and thrust bearings. Additional probes were used to 
measure relative displacements between various nonrotating parts of the bear- 
ing and casing assemblies. The capacitance probes are schematically indicated 
(in solid black) in Figure 1. The ends of the probes, mounted in the forward- 
thrust plate and the turbine-end journal-bearing pads, can be seen in Figures 
5 and 6 (similar probes were mounted in the compressor-end jounal-bearing pads). 
These probes, used to measure instantaneous film thickness, had a linear range 
of 5 mils. All of the capacitance- probe instrumentation channels had a flat 
frequency response to 8 kHz. 

The mcst significant transducer signals obtained in any given test were perma- 
nently documented on the multichannel tape recorder shown as part of the data 
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acquisition equipment in Figure 10. Dynamic responses were recorded on nine 
frequency-modulated (FM) channels, while a direct-record channel was used for 
voice recording. Most of the recording was done at 7.5 inches-per-second tape 
speed. Response of the FM channels’ at this speed was flat from DC to 2.5 kHz. 
Different combinations of transducer outputs to be recorded could be set up via 
the selector-switch panel. 

The semi-circular cutouts in two of the turbine-end journal pads, shown in 
Figure 5, provide clearance for two capacitance probes used to measure dynamic 
motions of the journal in orthogonal directions. (Similar probes were used in 
the compressor-end bearing). These probes had a linear range of 10 mils and 
were used to display the orbital motions of the journals. 

Crystal accelerometers were used at different points on the simulator casing 
and the test stand. Figure 3 shows several accelerometers mounted on the cas- 
ing, Frequency response of the accelerometers and related signal-conditioning 
equipment was flat from 2 Hz to 9 kHz, based on catalog data. 

Strain gages were mounted on several of the journal-bearing flexures, and on 
the forward-thrust-plate flexure, to measure dynamic strain in the flexures. 
Frequency response of the strain-gage system was flat to 3 kHz, based on cata- 
log information. 

Additional simulator transducers included pressure pickups to record simulator 
internal pressures, and a magnetic speed pickup to measure rotor speed. 

Figure 10 also shows the simulator control station and the instrumentation dis- 
play and recording equipment. At the simulator control station, rotor speed 
was indicated by an electronic counter located above the control panel. Speed 
was manually regulated by controlling air flow through the simulator drive 
turbine via a regulating valve located on the panel. Also located on the panel 
are two control valves for independent hydrostatic operation of the forward 
and reverse thrust bearings. During simulator startup, both bearings were 
operated hydrostatically (i.e., air pressure was supplied to the small orifices 
in the bearings). When sufficient speed was attained for hydrodynamic (self- 
acting) operation of the forward bearing, the hydrostatic pressure was shut off 


At this point the entire rotor became hydrodynamically supported, since the 
journal bearings are of the self-acting type. 

The test personnel normally monitored the eight dual-beam oscilloscopes shown 
in Figure 10 during simulator operation. Four of the scopes displayed the 
journal-bearing film clearances between each of the eight pads and the journals, 
while two scopes displayed the journal orbits relative to the casing. One scope 
monitored the clearance of the thrust runner with respect to both the forward 
thrust plate and the reverse-thrust stator. The eighth scope was generally 
used to monitor a shock or vibration input accelerometer, and a simulator 
response accelerometer. The locations of these accelerometers were arbit- 
rary, but they were usually chosen to be representative of the input, located 
near the test-stand base, and the output, located on the simulator itself. 
Although the above transducers were routinely monitored, certain circumstances 
required other sets of visual information. The selector-switch panel, shown 
in the center of the rack below the tape recorder, was used to obtain various 
combinations of transducer readout. A bank of charge amplifiers for the 
accelerometers was located just above the selector-switch panel. 

Test Method 


The vibration and shock tests were conducted independently. Each type of test 
was performed with and without isolation of the simulator casing, and with and 
without rotation of the shaft. Finally, the test combinations were conducted 
along each of three orthogonal axes of the simulator. 


The vibration testing was entirely sinusoidal in nature. For purposes of this 
first exploration with the turbocompressor simulator, the following schedule 
was used:: 


5-33 Hz 
33-140 Hz 
140-190 Hz 
190-2000 Hz 


140 mils double amplitude 
8 g ' s peak 

8 mils double amplitude 
15 g*s peak. 
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This vibration schedule was not a qualification test specification used to 
determine whether the equipment would fail or pass. It was, rather, a set of 



maximum vibration test objectives serving as a guideline while conducting the 
exploratory investigation of the vibration characteristics of the simulatorc 
The maximum vibration test objectives are outlined in Figures 11 and 12. The 
former is an acceleration plot which forms the basis of actual performance 
limits discussed in this report. Figure 12 is simply a restatement of the 
objectives in terms of displacement, which is a form commonly encountered in 
MIL-Specso The original specification was slightly modified to provide con- 
tinuity at the 190-Hz crossover point, which was previously set at 240 Hz. 

The shock-test objective calls for a half-sine pulse of 20-g peak acceleration 
and 10-milliseconds duration. This idealized shock pulse is shown in Figure 13^ 
together with the actual pulse obtained from the shock machine, and a haver sire 
approximation to the actual pulse. It is seen that the haversine is the better 
idealized approximation because it represents the case of elastic impact, which 
is nearly the situation of the shock table striking the elastomer pad on the 
anvil . 

The nonrotating vibration tests were conducted by sweeping through the 5- to 
2000-Hz range in accordance with the Figure 11 vibration amplitude schedule. 
During these tests, the signals from the strain gages mounted on the journal 
and thrust bearing flexures were continuously monitored to prevent overstress- 
ing of the flexures. For the rotating tests, an additional test criterion was 
imposed. The criterion required that the gas bearings always be maintained in 
a "contact-free*' operating condition such that the three-axis vibration survey 
could be completed with a minimum possibility of damaging the bearings. This 
criterion was satisfied by limiting the minimum operating film thickness of the 
bearings to 0.1 mil. 

To satisfy the above film thickness criterion during the rotating vibration 
tests, the vibration table was manually controlled such that the vibration 
conditions could be investigated slowly and carefully. If, as the 5- to 2000-Hz 
total frequency range -was manually traversed, minimum dynamic film thickness in 
any one of the bearings approached 0.1 mil, the vibration-table amplitude was 
reduced to prevent the minimum film value from dropping below 0.1 mil. In 
this manner, vibration amplitude schedules for "contact-free" bearing operation 
were determined as a function of frequency. These "contac t-f r ee" amplitude 
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schedules did not meet the vibration objectives, as is discussed in a later 
section of this report. 

Vibration-table amplitudes were also limited during the nonrotating tests to 
avoid excessive flexure stress. This problem would have to be corrected by 
redesign if the machine were to be flight rated. 

The shock testing was performed on the pneumatically raised vertical-shock 
machine shown in Figure 8. A typical shock pulse was obtained with a table 
drop height of about one inch. 
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AXIAL RESPONSE 


The shock and vibration tests were conducted along each of three orthogonal axes^ 
one axial and two transverse. There was very little observed interaction between 
the axial and the transverse vibration modes; therefore, they are considered to b 
separate and distinct throughout the remainder of this report. 

Analytical Models Used 

To be easily usable as a design tool, the analytical model of a physical system 
should be kept as simple as possible. At the same time, the model must be capabl 
of realistically describing, within acceptable limits of engineering accuracy, 
the effects that are of interest to the designer. The simplest possible model fo 
describing axial response of the simulator is the single-degree-of -freedom model 
shown in Figure 14. This model is comprised of the rotor mass, m, and the non- 
linear gas force, G(h,h), while the remainder of the simulator (the "base**) is 
considered to be entirely rigid. This one-degree-of -freedom model is adequate to 
describe the response of the rigidly-mounted simulator to axial shock, but it is 
inadequate (as discussed later) to describe the response to axial vibration for 
rigid-mount conditions. Furthermore, the model cannot be used to describe eithzrr 
the axial shock or vibration response when the simulator is mounted on isolators. 

Based on further consideration of the simulator design, as well as examination of 
the initial axial vibration test results, two three-degree-of -freedom analytical 
models were evolved which could satisfactorily describe the axial responses for 
all test conditions. These models are shown in Figures 15 and 16. 

Figure 15 shows the three-degree-of-freedom model used to calculate the axial 
response of the simulator when mounted on isolators. The base is taken to be a 
rigid platform to which the shock or vibration excitation is applied. The shock 
and vibration isolators are assumed to have axial symmetry, as is the entire 
simulator, including symmetrical application and distribution of all applied 
forces. The isolators are also assumed to be representable by a linear spring, 
k^ , and a viscous damper, c^ , connected in parallel as shown. The entire casing 
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is considered sufficiently rigid to be represented by the lumped mass, m^ „ 
thrust-bearing housing and the thrust-bearing flexure are considered as two simple 
springs in series having an equivalent spring constant, k^. Since these components 
are all metallic, they are considered to have negligible dampings The mass of the 
thrust plate, together with the effective combined mass of the thrust-plate flexure 
and the thrust-plate housing, are lumped into m^. The mass of the rotor is lumped 
as m^, which is connected to m 2 by the nonlinear gas-film force, G(h,h). G is a 
function both of the film thickness, h, and its time derivative, h. 

In order to use three lumped masses to represent the simulator on rigid mounts, a 
slightly different set of assumptions is required. The method of modeling the 
simulator on rigid mounts is shown in Figure 16. The base is considered to include 
the entire test fixture and simulator casing. The input shock and vibration is 
thus applied directly to the mounting flange of the thrust-plate housing. The 
stiffness of the thrust-plate housing is represented by k^ is now considered 

to be negligibly small), and m^ is taken as the effective mass of the housing. 

The stiffness of the thrust-plate flexure is lumped into k 2 , while the mass of the 
thrust plate is combined with one-third of the mass of the flexure to produce 
The rest of the model is identical of the previous case; that is, G represents the 
nonlinear gas forces while m^ represents the mass of the rotor. Again, complete 
axial symmetry is assumed with regard to properties and forces. 

Method of So l ution 

A computer program Wr3s written to solve the equations of motion which describe the 
r-zisponses of the systems shown in Figures 14, 15, and 16. The derivation of these 
equations is presented in Appendix A. The computer program is described in a 
gf-neral way by the flow chart of Figure 17, while the details of the program, 
me l ading program listing and operating ins tructions, ar e given in Appendix B. 

The first decision indicated on the flow chart of Figure 17 is whether to use the 
system representation shown in Figure 16 (rigidly mounted complete model). Figure 
15 (modt.l on isolators), or Figure 14 (the simple model). The user may then se- 
lect any of the following types of nonlinear functions to represent the load 
capacity and the damping coefficient characteristics of the thrust-bearing gas film 
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a) An exponential function of the form: 


where h = mean thickness of thrust bearing gas film 
e = base of natural logarithms, 

and A and B are constants selected by the user to give the best practi- 
cal approximation to the nonlinear characteristics of the gas film as 
predetermined by applying rigorous lubrication theory. 


b) 


A power function of the form: 


Ch 


-D 


where C and D are the arbitrary constants used to obtain a curve fit . 


The program will handle any of the four different types of forcing functions shown 
in Figure 17, which may be expressed either as base displacement or base accel- 
eration. These forcing functions include sinusoidal vibration and three types 
of shock pulses. 


The haversine pulse is the best simple idealization of the main shock pulse pro- 
duced by the type of shock machine used, which relies on quasi-elastic impact to 
generate the required levels of acceleration. The half-sine pulse, as called for 
in the specification, is closely approximated by the haversine pulse, since the 
latter falls within the half-sine shock pulse tolerance limits given in Figure 
516-2 of MIL-STD-810B , dated June 15, 1967. The half -sine pulse is physically 
unrealizable because the discontinuous derivatives at the end points imply a 
system with zero response time. The arbitrary shock-pulse excitation option adds 
considerable flexibility to the program. This option was used to describe the 
actual shock-machine pulse since it was necessary to include the acceleration 
values during the drop interval of the table. 

In order to calculate the response of the simulator using the computer program^ 
it is necessary to assign values to the parameters appearing in the equations of 
motion, as well as to define the nonlinear load capacity and damping character- 
istics of the thrust bearing gas film. Figure 18 shows curves of the static 
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load carrying capacity of the thrust-bearing as a function of film thickness. 

The solid curves x^ere calculated by gas lubrication theory using Reynolds equa- 
tion iwith appropriate boundary conditions for our particular bearing design 
operating under the conditions stated in the figure. The upper curve is based on 
calculations assuming flat and parallel surfaces. The effect of crooning (convex 
distortion) of the thrust plate results in a loss of load carrying capacity as 
is shown by the lower solid curve. Based on experience, it is not unreasonable 
to expect the amount of distortion represented by this lower curve, but we note 
that the curve still does not intersect a known load capacity point of the thrust 
bearing (10,5 pounds at a measured film thickness of 0.9 mils). Therefore, an 
exponential curve was assumed which was parallel to the original curve but dis- 
placed by a sufficient amount to agree with experimental evidence. The additional 
displacement is justifiable because of thermal effects. This dashed curve is 
the one used in the computer program. 

Figure 19 shows the coefficient of viscous damping for the thrust bearing as a 
function of frequency and film thickness. The solid curves, which are rigorously 
correct only for **small" amplitude sinusoidal motions, are also based on gas 
lubrication theory under the same conditions as noted on Figure 18, The dashed 
curve is the power function approximation to the damping characteristic that was 
used in the computer program to calculate shock response, as well as to calculate 
vibration response at 200 Hz and below. Both the damping coefficient and the 
load capacity approximations begin to diverge from the theoretical curves at 
bearing clearances in excess of 2 mils, a fact which should be realized in inter- 
preting some, of the shock results presented later in this paper. 

The above approach of treating the gas film has been made deliberately empirical. 
Alternatively, a more rigorous attack would have been to solve the Reynolds 
equation of gas dynamics [4], which is a nonlinear, partial-differential equation, 
to produce the thrust-bearing gas film force on a point-to-point basis. However, 
the rigorous solution would have required an increase in computer time by a con- 
siderable amount over the method evolved here. The intent here was instead to 
develop an economical design tool to determine the axial dynamic response of the 
rot or -bear i ng system. 
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In addition to the thrust bearing load capacity and damping characteristics, 
there are several other parameters required for axial response computations, 
such as the coefficients for isolator damping and stiffness, and for the housing 
stiffness. 

The isolator properties were initially determined by taking manufacturer's data 
in the vicinity of peak transmissibi lity and performing a curve fit with the 
isolator representation shown in Figure 20. The good correlation obtained with 
the simple linear isolator representation is deceptive, however, when applied to 
field test conditions, as we shall show later. The principal deficiency in the 
simple model lies in its linearity. 

Figure 21 presents the measured static- load characteristics of an individual 
shock-and- vibration isolator, presumed to be representative of the four that 
were used to isolate the simulator in the vertical attitude. It reveals the 
inadequacy of the simple linear representation. Although the isolator spring 
constant is essentially linear within its design static- load range, it exhibits 
strong nonlinearities under the conditions of shock previously stated. The total 
static load of the simulator was about 120 pounds, which initially loaded each 
isolator to 30 pounds. As the isolators go into tension shortly after shock- 
table release, they exhibit a spring-hardening characteristic. 

This nonlinear load characteristic is due to two factors “ the construction of 
the isolators, and the nonlinear modulus of the rubber used. The construction 
of an individual isolator is shown in Figure 22. The initial preload is carried 
by the concentric springs E and F. This load is applied through the specimen 
support mount G. For small deflections there is little compliance of the rubber 
spring inherent in the cylindrical section of A. The initial damping is primarily 
due to the relative motion between the split ring, C, which is loaded against 
the cup, B, via the split spring, D. After springs E and F have compressed, 
the rubber column of A provides some compliance and gives rise to a spring 
hardening characteristic when the deflection exceeds 0.5 inches. The rubber 
also provides the damping mechanism at these large deflections. The snubbers 
protruding above the retaining plate A transmit the load directly through the 
rubber column at very large deflections. 
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The stiffness of the thrust-plate housing is an important parameter required to 
calculate the response of the three- degree- of- freedom representation of the simu- 
lator on rigid mounts o This stiffness was initially calculated with the aid of 
Reference 5 for the three models shown in Figure 23* The spread of stiffness 
values shown was sufficient to require actual measurement to obtain sufficient 
accuracy and confidence* The methods of stiffness measurement along with the 
results obtained are presented in Figure 24* The method shown in the lower half 
of the figure is thought to be more representative of the actual attachment situa 
tion (the housing was bolted to the casing with three No. 10-32 screws) and there 
fore the lower stiffness value of 533300 Ib/in was used in the response calcula- 
: i ons * 


I nitial Results 

The static measurements and design values established in the foregoing section 
provide sufficient data to utilize the axial response computer program to gener- 
ate some initial results* 

A typical computer result is shown in Figure 259 which was produced by the com- 
puter program on an automatic plotter* The legend and axis scaling are entirely 
performed by the program* This particular response was obtained with a three- 
degree- of- freedom representation of the system* The lower right hand corner of 
the figure indicates the time of contact 9 at which point the rotor bottoms 
against the. thrust plate* No subsequent calculations are performed beyond this 
point and the computer program automatically terminates because the analysis has 
no provision for bearing contact* 

Figure 26 is also a typical example of the computer program output* This timeg 
the excitation is a sinusoidal vibration applied to the single-degree-of- freedom 
representation of the system* Note that the nature of the solution required 
calculating the transient response all the way up through the steady-state solu- 
tion* Scratiny of this steady-state response after 0*04 seconds in Figure 26, 
reveals the calculated effect of the nonlinear gas film characteristic which 
amounts to an asymmetrical periodic motion about the equilibrium film thickness* 
As we shall see later, the measured asymmetry was even more pronounced than the 
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calculated value 


Figure 27 shows a typical set of data that were obtained during the vibration 
test program. The oscilloscope photos presented were taken at 110 Hz with an 
input vibration amplitude of approximately 0.1 g*s peak. This frequency, 
while slightly below the resonance point, was representative of the resonance 
region. The photos show the modal behavior which prevailed up to about 300 Hz. 
The upper oscilloscope photo reveals that the gaps on each side of the thrust 
runner are moving out of phase and have equal total excursions. The middle 
photo shows that the reverse thrust plate has virtually no response relative 
to the housing except for a slight, high-frequency ripple corresponding to the 
running speed of the rotor. 

The bottom photo shows that diametric positions of the thrust-bearing gas film 
are responding in phase and are of about equal amplitude. These three photos, 
taken collectively, lead to the conclusion that this modal behavior consists 
of axial motion of the rotor with no noticeable wobble and no detectable 
participation of other components, such as the thrust plate or the housing. 

The undamped first natural frequency of the simulator was determined experi- 
mentally by observing the Lissajous pattern formed by the signals of a gas- 
film- thickness probe and a probe measuring base motion. Figures 28 and 29 
indicate the method used. The Lissajous patterns ot Figure 28 were obtained 
from the setup of Figure 29. The latter figure shows the capacitance probe 
"looking” at the displacement of the casing relative to the rigid piece of 
channel iron attached to the stationary portion of the vibration shaker. The 
signal obtained was fed into one channel of a dual beam scope while the gas 
film thickness was fed into the other channel. The relative amplitudes were 
adjusted to produce a circle at the 90-degree phase position to facilitate 
identification of this phase condition, which is depicted by the center 
Lissajous pattern in Figure 28. All three patterns were obtained with the 
same amplitude setting. The departure from ideal ellipses is due to the wave- 
form departure from an ideal sinusoid. As noted on Figure 28, the undamped 
first natural frequency of the rotor was 120 Hz. 
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The computer program was used to predict both undamped and damped first natural 
axial frequencies of the rotor on the thrust bearing. An equivalent linear 
stiffness of the gas film was determined to be 27,000 Ib/in by taking the deriva- 
tive of the load curve at the point of static equilibrium corresponding to a 
rotor weight of 10.5 pounds. These values of stiffness and weight give a calcu- 
lated undamped natural frequency of 158 Hz for a single-degree-of-freedom system. 
The discrepancy between this and the measured, undamped natural frequency (120 Hz) 
indicates that a single-degree-of-freedom representation is inadequate. The in- 
adequacy cannot be ascribed to the gas- film nonlinearity because the axial re- 
sponse computer program yielded an undamped natural frequency of 150 Hz for the 

single-degree-of-freedom nonlinear representation. 

A t hr ee-degree-of- freedom representation of the system was therefore considered. 
The equivalent lumped weight of m 2 was determined from a design drawing to be 
1.4 pounds, and the equivalent value of m^ was found by weighing the parts to be 

0.4 pound. The measured value of was 53,300 Ib/in while the design value of 

670,000 Ib/in was used for k 2 * These values produced calculated undamped natural 
frequencies of 128 Hz, 685 Hz, and 4,654 Hz for a linear three- degree- of- freedom 
system. The first of these frequencies shows good agreement with the 120 Hz 
measured value. 

Beyond the first natural frequency, two other distinct axial vibration modes were 
observed during the vibration test program. One of these is shown in Figure 30 
which displays data from the same probes described in Figure 27. Figure 30 shows 
the modal behavior that occurred in the second distinct resonance region. These 
photos were taken at 605 Hz with about 0.5 g''s peak input. The three displays 
are not simultaneous; therefore, amplitude comparisons cannot be made from photo 
to photo because of slight variations in vibration- table amplitude. Again, the 
mode is representative of the frequency region, and not necessarily at the 
resonance point. The upper photo shows that the probe readings on either side 
of the thrust runner are out of phase, and that the forward thrust-bearing gas 
film exhibits the greater of the excursions. The middle photo reveals that there 
is now some relative motion between the housing and the thrust plate. It 
further reveals that this relative motion is approximately 180 degrees out of 
phase with the thrust-bearing gas film. It is also worth noting that the second. 
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undamped, natural frequency for the linearized model was 685 Hz and lies with- 
in this second, observed resonance region. Furthermore, the gas film's total 
excursion is about equal to the sum of the other two relative motions, which 
implies that the thrust-plate housing itself is not participating in this 
mode. The lower photo, again, shows no significant wobbling of the thrust 
runner. These three photos, taken collectively, imply that the second mode 
consists of axial motion of both the rotor and the thrust plate, with the 
rotor exhibiting the greater amplitude. 

A third mode, shown in Figure 31, was observed in the vicinity of 870 Hz. 
Oscilloscope photos obtained were similar to those of the second mode except 
that the reverse thrust bearing gap now displays an excursion one-third as 
large as the forward- thrust-bearing gas film variation. This implies that 
this third mode involved some flexing of the housing, or some relative motions 
between the housing, the casing, and the reverse- thrust plate. At any rate 
this third mode is not as readily explained as were the first and second modes. 

The waveform of the thrust bearing response seen in Figures 27, 30, and 31, 
although not a pure sine wave, is reasonably sinusoidal. This condition pre- 
vailed throughout the frequency spectrum above 80 Hz. At lower excitation 
frequencies, however, the wave form became asymmetrical due to the film non- 
linearity. This is shown in an idealized way in Figure 26. 

Figure 32 shows that the actual asymmetry is considerably more severe. In 
fact, the waveform is highly irregular indicating the influence of effects not 
calculated by the simple theory. This study does not delve into the causes 
of the erratic waveform at low frequencies, but merely documents observation 
of the effect . 

Results of 5- to 2.0 0 0-Hz Vibration Tests 

Figures 33 and 34 show the input vibration levels which could be sustained, 
without and with isolators, while maintaining contact-free operation of the 
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thrust bearing. It is seen that considerably higher levels of input vibration 
could be sustained when using isolators. However, even with isolators, the desired 
levels of input acceleration were not achieved at low vibration frequencies due 
to resonance of the isolator-casing system at 12 Hz, and resonance of the thrust 
bearing at 13C Hz. Because of our decision to maintain contact-free bearing oper- 
ation during the vibration tests with shaft rotation, we know only that bearing 
contacts would probably have occurred (with or without isolators) in the frequency 
regions where the objective input levels were not achieved. Further testing would 
be needed to determine whether sustained, vibration-induced, rotor -bearing contacts 
would significantly degrade bearing performance. 

A series of computer runs was conducted using the nonlinear load and damping char- 
acteristics of the gas film to explore the region around the first natural frequency. 
The results are presented in Figure 35, which compares the vibration response of 
I he three“degree-of-freedom system with that of the single-mass system and with 
actual measured values. In order to make these comparisons as consistent as pos- 
sibles all the measurements were taken with a 0.2 mil peak-to-peak shaker amplitude, 
the same value used in the computer runs.* The single-degree-of -freedom results 
predict a peak at about 160 Hz. Observation of the phase angle between the input 
and response indicated a calculated undamped natural frequency of about 150 Hz. This 
is consistent with simple linear theory for a single mass system which predicts 
that the effect of damping is to shift the peak response of the relative displace- 
m.-nt above the undamped natural frequency [6] c 

The thret-degree-of" freedom model results in considerably better agreement with 
measured values than does the single^mass approximation. The three-mass model pro“ 
daces a p<=ak response at 130 Hz which is consistent with the previously mentioned, 
m'.asur(-,d, undamped frequency of 120 Hz. The effect of the gas-film non-linearity is 
to lower the first natural frequency of the three-mass model by a few percent. The 
foregoing is summarized in Table I. 

The final Stories of axial vibration tests were performed with a nonrotating shaft, 
with and without vibration isolators. For either condition of mounting, the spe- 
cified vibration excitation could be imposed over the complete 5- to 2,000-Hz range 


*A common input amplitude was necessary because the natural frequencies of non- 
linear systems are amplitude dependent. 
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TABLE I 


Comparison of Values for the First 
Natural Frequency of the 
Simulator Rotor-Bearing System in the 
Axial Response Mode (hydrodynamic operation at 38,500 rpm) 



Undamped Natural Frequency 
(Hz) 

Resonant Frequency 
(Hz) 



Analytical Model 

Calculated 

Ca Iculated 

1 mass linearized 

158 

— 

•k 

1 mass nonlinear 

150 

160 

3 mass linearized 

128 

— 

k 



3 mass nonlinear 

125 

130 


Measured 

Measured 


120 

Between 120 & 150 


The undamped natural frequency for the nonlinear system is actually a func- 
tion of response amplitude; frequency decreases with increased amplitude » 
Values stated are based on 0.2 mils peak-to-peak input. 
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without exceeding safe dynamic stress levels in the thrust or journal bearing 
flexures. After each nonrotating vibration sweep, the simulator was started up 
and operating clearances of all bearings carefully checked. No significant change 
in bearing system performance was detected as a result of the nonrotating vibra- 
tion tests. 

Results of Axial Shock Tests and Calculation^ 

Figure 36 shows a typical set of test data displays for axially imposed shock with 
the simulator mounted on isolators and the shaft rotating at 39,000 rpm. The top 
traces of photos A and B show the relative motions between the rotor and journal- 
bearing pads. Photos C and E are accelerometer traces. Photo D depicts the 
casing-to-thrust-plate motion in the top trace, and the relative motion between 
the casing and the reverse side of the thrust runner in the bottom trace. 

Figure 37 shows a comparison of measured and calculated shock response of the 
thrust-bearing film thickness for the rigidly mounted simulator with the shaft 
rotating at 38,500 rpm. The calculated results were obtained using the axial 
response computer program described earlier. The initial response is a gradual 
increase in film thickness occurring as a result of the shock-table drop interval 
preceding the impact. As the impact occurs, however, the gas film undergoes a 
reversal, followed by a rapid decrease in film thickness. It is of interest to 
note that a slight hump appears in the calculated response curve during the final 
period of the gas-film collapse just prior to contact. It was initially thought 
that this hump was due to an irregularity of the shock pulse, which was modeled 
as a straight line during the drop interval. Therefore, the response was recal- 
culated using the closed-form expression for a haversine as a pulse approximation. 

Figure 38A shows the calculated undamped response of the thrust-bearing gas film 
using the haversine shock pulse approximation. Figure 38B is a repeat of the 
same calculation, but with the normal nonlinear thrust bearing damping restored. 

It is concluded from these two calculations that the inflection in the thrust 
bearing film thickness curve in Figure 38B results from the nonlinear property 
of the gas film rather than from a quirk of the shock pulse. This characteristic 
is of significant importance because, as illustrated by Figure 38B, the film 
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damping results in a slower approach velocity of the runner just prior to contact. 
This may result in far greater shock resistance for this type of gas bearing than 
one might otherwise expect. 

Figure 39 furnishes further experimental evidence of the precontact inflections 
together with further comparisons of the calculated shock response during collapse 
of the thrust-bearing gas film. Before film collapse, the three-degree-of-freedom 
model and the single-mass model yield similar calculated responses. The essential 
difference during the collapse interval, shown in Figure 39, is that the three- 
mass model predicts a simple inflection rather than the hump exhibited by the 
single-mass model. The simple inflection bears a closer resemblance to measured 
results and also appears for the three-mass model when a haversine approximation 
is used to represent the shock pulse. It is interesting to note that there is 
only a spread of about one millisecond in the contact time regardless of the model 
used to represent either the simulator or the shock pulse. 

Figure 40 gives further explanation of the low velocity thrust bearing contact. 

The nonlinear film stiffness and damping functions are such that the damping 
ratio of the rotor-gas-film system becomes very large at small values of film 
clearances . 

Figure 41 shows the measured shock response of the thrust-bearing gas film with 
the simulator mounted on shock-and-vibra tion isolators. For comparison, calcu- 
lated results are shown using two different values of isolator stiffness. The 
agreement between calculated and measured response is considerably worse than 
was indicated for the rigidly mounted case shown in Figures 27 and 39. There 
is only a general agreement with respect to the gross behavior when the stiffer 
value of k^ (16,000 Ib/in) is used. For example, an initial peak occurs, with a 
noticeable discrepancy in amplitude, in both measured and calculated responses at 
about 70 milliseconds after shock-table release. This comment also applies to the 
dips occurring at about 77 milliseconds and the subsequent rise. Beyond 82 milli- 
seconds, the calculated and measured responses diverge considerably. At about 
87 milliseconds, the experimental data indicates that the thrust runner contacts 
the reverse-thrust plate. The computed response, however, does not even approach 
a contact amplitude until 110 milliseconds after table release. 
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The poor correlation evidenced in Figure 41 stands in sharp contrast to the close 
agreement between measured and calculated results shown in Figure 37 „ The funda- 
mental difference in the measured behavior in the two cases is as follows: when 

rigid mounts are used (Figure 37), the initial collision occurs between the 
thrust runner and forward-thrust plate during the table-impact interval, with 
subsequent collisions between the thrust runner and reverse thrust plate; however, 
when shock-and -vibration isolators are used (Figure 41), the initial collision 
occurs between the thrust runner and reverse-thrust plate befor e the main pulse 
of the table-impact interval. 

The computer program was not able to predict this contact with the reverse- thrust 
plate even when using a meticulous representation of the shock pulse with measured 
values of acceleration obtained every millisecond during the table-drop interval. 
The manner in which the isolators produce this sudden rise in film thickness to 
cause reverse contact can be qualitatively discussed as follows: initially, the 

shock table is released and begins to fall; the inertia of the low-frequency 
system (consisting of the simulator on the isolators), along with release of the 
preload on the isolators, cause an inherent response delay. Eventually the iso- 
lators are fully extended and begin to "snap back" so as to impart a sudden down- 
ward acceleration to the simulator casing. The r ever se-thrust plate then collides 
with the thrust runner, which has its own inertia by virtue of the 10.5-pound 
rotor, as we shall see presently. An accelerometer attached to the casing veri- 
fied that this was indeed the sequence of events. 

The poor quantitative correlation shown in Figure 41 results from the inadequacy 
of the simple linear model used to represent the isolators. As Figure 41 indi- 
cates, adjustment of the linear spring constant, k^ , will not, by itself, give 
good correlation of results. 

When table Impact occurs at about 107 milliseconds after table release, the base 
acceleration reverses direction and tends to drive the isolators into compression. 
The; subsequent g-levels on the simulator become large enough to produce . compr es- 
si.ve loads in excess of those shown in Figure 21. These loads, which are of th'=; 
order of the product of the simulator weight (120 pounds) and the g-level (17 
g^s), place the operating point above the steep region near the upper right hand 
corner of Figure 21, In fact, the effective stiffness is in excess of the value 
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indicated because, as is shown in Reference 7, the dynamic modulus of rubber, having 
the hardness measured for the subject isolator (Durometer hardness = 67), is about 
twice the static modulus. 

Figure 42 provides verification that the linear isolator modeling is the primary 
cause for the discrepancy shown in Figure 41. In order to get the results shown 
in Figure 42, the computer program was run using the measured casing accelerometer 
as the shock input, thereby effectively bypassing the isolators. By so doing, it 
was possible to analytically predict the salient features of the thrust-bearing 
gas-film response up to the point where the thrust runner collides against the re- 
verse-thrust plate. Beyond this point the computer results would not be valid be- 
cause the program contains no provision for bottoming in the form of bearing sur- 
face contact. 
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TRANg VERSE RESPONSE 


As was previously mentioned, the axial and transverse responses exhibited very 
little interaction. Hence, the discussion of transverse response deals primarily 
with rotor and bearing motions in the vicinity of the journal bearings, with less 
relative significance attached to the thrust-bearing response. 

The shock and vibration testing was conducted in two transverse directions, which 
are distinguishable by the dead- load position of the rotor with respect to the 
journal bearing pads. Each journal bearing consists of four identical pads 
supported by individual flexures. Two adjacent pads are supported by substan- 
tially stiffer flexures than are the remaining two. The two former pads will be 

referred to as "rigid" pads, and the two latter as "compliant" pads. The pad 

arrangement is seen in cross section in Figure 2. The two transverse directions 

in which dynamic loading was applied are defined as follows: 

1. The "Tl" direction, this being transverse rotor excitation in a direc- 
tion parallel to the plane containing the rotor centerline and passing 
midway between the pivots of the twp "rigid" pads; 

2. The "T2" direction, this being transverse rotor excitation in a direc- 
tion peip^ndicular to the "Tl" direction (i.e., parallel to the plane con- 
taining the rotor centerline and passing midway between the pivots of 
one "rigid" pad and one "compliant" pad). 

Analysis of transverse response was performed for the Tl direction only, and the 
only forcing function considered was that of steady-state vibration. 

Analysis Model For Transverse Vibration Calculations 


A distributed- parameter model of the rotor-bearing assembly was used to analyz.e 
system response to transversely applied vibration. This model, implemented in 
a MTl computer program, describes the rotor and casing in much the same way as 
one would consider a beam of varying cross section having continuously distributed 
mass and modulus of elasticity. The basic analysis is thoroughly described in 
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Reference 8 and will not be discussed in detail here. 


The more significant capabilities of the computer program are pointed out in 
Figure 43 which specifies the parameters used to calculate transverse vibration 
response. In addition to the large number (theoretically limitless) of degrees 
of freedom arising from modeling the casing and rotor as continuous finite ele- 
ments, other vibration modes result from the isolators and from the motions of 
the journal bearing pads. There are a total of eight pads in the journal bearing 
system and each pad has four degrees of freedom as shown in Figure 44. The 
translational mode and each of the three rotational modes may be excited, there- 
by giving rise to a number of frequencies at which various resonances might 
occur. The various frequencies that are associated with the particular bearing 
design used in the simulator are presented in tabular form in Table II. 

A cross-section view of the rotor- bearing system simulator is shown in Figure 45. 
The details concerning each of the significant features shown are presented in 
Reference 2 and will not be dealt with here. We will consider primarily the 
analytical modeling of the rotor and casing as shown in Figure 46. This latter 
representation was obtained from the drawing of Figure 45 by dividing the rotor 
and casing into segments of an inch or two in length. At each segment, identi- 
fied by a station number, values were assigned to the equivalent physical dimen- 
sions of that particular segment. In this way, complex regions could be very 
closely approximated by short sections of uniform, cylindrical cross section. 

The journal-bearing pad mass and inertia properties were also specified along 
with the gas film characteristics. It should be emphasized that the computer 
program used to determine transverse response is a sophisticated tool that has 
evolved from previously successful programs to calculate unbalance response of 
flexible rotors. 


The gas film stiffnesses and damping characteristics of each journal bearing pad 

may be expressed in terms of eight linear coefficients. Thus, K refers to the 

xy 

change in force in the x-direction due to an incremental motion in the y-direc- 

tion. Similarly, B refers to the change in viscous damping force in the 

yx 

y-direction due to an incremental velocity in the x-direction. These coefficients 
are commonly treated in normalized form as follows: 
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i = 22S 

yx P LD 
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where C is the radial (machined) pad clearance, P is the ambient pressure, and 

a 

L and D are the nominal bearing length and diameter, respectively. These coeffi- 
cients are not constants, but vary both with journal eccentricity and with fre- 
quency of applied vibration. This frequency is commonly normalized by dividing 
it by the rotor speed. The resulting variation of the normalized coefficients 
is shown in Figure 47- 

Another set of properties required by the analysis consisted of the isolator 
stiffness and damping values. An attempt was made to obtain these experiment- 
ally by vibrating the simulator transversely, with isolators, in the vicinity 
of the natural frequency of the isolator system. Accelerometer readings were 
obtained simultaneously to represent the input and output of the isolators. The 
former accelerometer was located on a reasonably rigid region of che base while 
the latter was located on the isolated simulator. 

The results of these measurements, along with a series of curve- fit approxima- 
tions, are presented in Figure 48. Comparison with Figure 20 reveals that, 
again, the measured values are generally below the curves obtained by a linear 
system representation. The effect of increasing the damping coefficient is to 
lower the analytical transmissibility values. The damping coefficient was 
adopted as G = 33. It was felt that higher values tended to flatten the curve 
excessively, thereby diluting the pronounced peak that was actually observed. 

As in the vertical case, the nonlinearity of the isolator introduces consider- 
able error in attempting to describe the response with a linear model. At any 
rate, we have, at this point, determined the basic data that are required as 
input to the computer program. They are: 

1) Model of rotor and casing 

2) Journal bearing properties and gas film characteristics 

3) Isolator properties. 


CK 




P LD 
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Table II 


Calculated Natural Frequencies of 
Journal Bearing Pads 


Frequency 


Bearing 

F i exur e 

Mode 

(Hz 





Compressor 

Compliant 

Radial 

2100 

Compressor 

Compliant 

Pitch 

880 

Compressor 

Compliant 

Roll 

1170 

Compressor 

Compliant 

Yaw 

245 

Compressor 

Rigid 

Radial 

3580 

Compressor 

Rigid 

Pitch 

1300 

Compressor 

Rigid 

Roll 

1720 

Compressor 

Rigid 

Yaw 

370 

Turbine 

Compliant 

Radial 

1460 

Turbine 

Compliant 

Pitch 

625 

Turbine 

Compliant 

Roll 

860 

Turbine 

Comp liant 

Yaw 

177 

Turbine 

Rigid 

Radial 

2820 

Turbine 

Rigid 

Pitch 

1000 

Turbine 

Rigid 

Roll 

1400 

Turbine 

Rigid 

Yaw 

290 


I 


*Shaft rotating at 39,000 rpm in air at 14.7 psia and 70®F. 



Results of the transverse vibration response calculations are discussed start- 
ing on page 38. 

Transverse Vibration Tests 


The simulator transverse vibration tests were conducted in much the same way 
as were the axial vibration tests. That is, all bearing clearances were 
monitored while conducting a frequency sweep and continuously adjusting the 
amplitude to maximum values within the predetermined test envelope. 

The most pronounced effect of vibrating the system transversely was to induce 
orbital motion of the rotor that was observed at both journal bearings. Two 
distinct orbital resonances were observed during the tests. The first of these 
occurred at 166 Hz, and the precise nature of the rotor orbits is shown in 
Figure 49. Figure 49a shows the relative displacement between the rotor and the 
casing at the position of the compressor journal bearing while Figure 49b shows 
the simultaneously occurring orbit in the vicinity of the turbine journal 
bearing. 

The other distinct orbital resonance occurred at 225 Hz; the response is shown 
in Figure 50. As before. Figure 50a shows the rotor- to-casing displacement at 
the compressor journal bearing while Figure 50b shows the corresponding turbine 
bearing orbit taken simultaneously. Comparisons and further discussion of these 
orbits will be deferred to the next section. 

Discussion of Transverse Vibration Tests and Calculations 


It has been shown, in Figures 49 and 50, that rotor-bearing system resonances 
were excited at 166 and 225 Hz. These were subsequently identified as the 
rotor first and second rigid-body critical speeds, respectively. 

A previously developed MTI rotor lateral critical speed computer program was 
used to calculate the modal shapes corresponding to these two frequencies and 
the results are presented in Figure 51 which is based on Figure 11-21 of Ref- 
erence 2. Examination of Figure 49 indicates that the first rigid-body mode 
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has a noticeably greater amplitude at the turbine journal bearing compared with 
the compressor bearing* This is predicted qualitatively in Figure 51. 

The measured data of Figure 50 shows that, conversely, the second rigid-body mode 
produces a greater amplitude at the compressor bearing. This also is exhibited 
qualitatively in Figure 51. 

It is not surprising that, at each frequency, the ratio of measured maximum or- 
bital values is lower than what Figure 51 would indicate. The apparent discrep- 
ancy arises from the fact that pure modes are not excited in isolated ways. The 
two modal frequencies are quite close together and there is probably a substan- 
tial participation of the second mode at the lower forcing frequency and vice 
versa . 

When the previously described rotor-casing model and bearing properties were used 
in the computer program to calculate transverse response, the orbits shown in 
Figure 52 resulted. The input steady-state sinusoidal vibration is 12 g*s peak 
amplitude at 166 Hz, the same condition used to obtain the measured orbits of 
Figure 49. These measured orbits are presented in Figure 52 in an idealized form 
in order to make a direct comparison between measured and calculated paths. We 
observe only qualitative agreement in the sense that the major dimension, or 
orientation, lies roughly in the same general direction for the calculated and 
measured orbits. It should also be noted that the larger orbit occurs in the 
vicinity of the turbine journal bearing for both the calculated and measured 
responses. 

Figure 53 is similar to Figure 52 except that the applied vibration condition is 
now 12 g’s peak at 225 Hz, the critical speed of the second rigid body mode. 
Again, we get the relatively larger calculated orbit occurring at the compressor 
journal bearing as was observed in Figure 50 and predicted in Figure 51. A very 
significant discrepancy, however, is noted in Figure 53: the orientations of 

the measured orbits are off by about 90 degrees at the compressor journal bear- 
ing and by about 45 degrees at the turbine journal bearing. 

One would expect that the rotor orbits would generally have their major axes 
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aligned with the direction of applied vibration. A slight perpendicular com- 
ponent would arise from gyroscopic effects due to a slope along the shaft which 
might occur in the vicinity of a heavy disk such as the turbine or thrust bear- 
ing. Indeed such forces produce the slight elliptical broadening seen In the 
calculated orbits in Figures 52 and 53. Some horizontal vibration component 
was probably present, although no measurement of it was attempted. In each 
cascj, however, the measured response component perpendicular to the direction 
of motion was observed to be noticeably greater than that calculated. The dis- 
crepancy is documented here, but no further attempt is made to explain it. 

Although the computer program for transverse response yielded poor amplitude 
and orbital orientation data, it did yield remarkably good predictions of the 
resonant frequencies. This is evidenced in Figure 54 which shows the calculated 
peak amplitude of the relative displacement between the rotor and the housing 
as a function of frequency in the vicinity of the first and second rigid body 
critical speeds. The measured amplitudes at the observed critical speeds are 
spotted in on the graph, to show the poor amplitude agreement and the excellent 
frequency correlation. 

The first parameter that was examined in an attempt to account for the amplitude 
discrepancy was the damping coefficient of the isolator. A decrease in damping 
causes a decrease in the transmissibi lity of the isolator within the frequency 
range of interest. The sensitivity of the response to this parameter is shown 
in Figure 55. Comparison between Figures 54 and 55 reveals that isolator damp- 
ing is a relatively sensitive parameter for high damping values but it becomes 
less sensitive at the lower values. In fact. Figure 56 shows that regardless of 
how small we consider the isolator damping coefficient to be, it is impossible 
to explain the amplitude discrepancy at the turbine journal bearing solely in 
terms of isolator system damping. 

The next parameter to be investigated was the damping value of the journal bear- 
ing gas films. Figure 56 shows that film damping is a very sensitive parameter 
throughout a wide range of values. The figure shows that it was possible to 
account for the discrepancy in amplitude at the compressor bearing by doubling 
the damping in this bearing. The turbine journal bearing, on the other hand. 
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requires a damping increase by a factor of 4 or 5 to bring the calculated ampli- 
tude down to the observed experimental value. 

It thus appears that a source of discrepancy in the results is damping in the 
turbine region that is not properly accounted for in the computer program. This 
overlooked damping may possibly arise from the dissipation of energy in the 
turbine labyrinth seals. 

Another interesting effect that was encountered as part of the transverse re- 
sponse calculations was identification of the first bending mode of the casing, 
the resonance of which corresponds to a forcing frequency of about 1440 Hz. The 
mode shape of the casing is shown in Figure 57 which presents the rotor-to- 
housing relative clearance for a one-mil amplitude sinusoidal input to the simu- 
lator on rigid mounts. 

The other calculated natural frequencies, which are associated primarily with 
transverse response, are the individual journal bearing pad resonances. These 
frequencies are tabulated in Table II, and the terminology for the various vibra- 
tion modes was previously defined in Figure 44. The majority of the pad reson- 
ances were not severely excited during the test. In fact, the only pad resonance 
that was clearly identified was the radial mode of one of the compliant pads in 
the turbine journal bearing. This resonance was so severe that the flexure 
failed during testing at 1365 Hz. It is interesting to note that this frequency 
is about 6.5% below the theoretical value as presented in Table II. The de- 
tails concerning the failure are discussed in Appendix D. 

Figures 58 through 61 summarize the vibration test results for transverse re- 
sponse in much the same way as Figures 33 and 34 did for the axial testing. 

Figure 58 presents the maximum safe operating levels for the rigidly mounted 
simulator with vibration applied in the T1 direction. No attempt was made to 
identify the individual component resonances while taking this data, but there 
was a general tendency toward severe reaction by the compressor journal bearing 
at the lower frequencies while the turbine journal bearing participated in the 
higher frequency region (above 800 Hz). 

Figure 59 shows the response with vibration again applied in the T1 direction, 
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but this time the simulator was mounted on shock and vibration isolators « It is 
evident that the isolators performed their function well except for two narrow 
regions. The first of these occurred at the natural frequency of the isolators 
(between 6 Hz and 15 Hz). The second critical region occurred between 190 Hz 
and 235 Hz. Within this range of frequencies the rotor exhibited the same type 
of orbital motion discussed previously. 

Figure 60 presents the maximum safe operating levels when vibration was applied 
in the T2 direction without isolators. The same general response descriptions 
for the T1 direction are also applicable here. In fact, the only essential dif- 
ference is that there is some loading asymmetry inherent to the T2 direction. 

Figure 61 shows that, once again, the isolators performed their intended func- 
tion. In fact, even the orbital motions in the region between 190 and 235 Hz 
remained small enough to permit testing at the maximum test envelope. This is 
probably due to more of the orbital energy going into the perpendicular com- 
ponent because of the asymmetrical loading. The only region that could not be 
met was a one-Hz band width between 7 and 8 Hz. 

We see from Figures 62 and 63 that isolators are also required in the presence 
of transverse vibration with a nonrotating rotor. With the simulator rigidly 
mounted, vibration within the frequency range from 190 to 235 Hz produced large 
amplitude levels in the journal bearing flexures, regardless of whether the 
vibration was applied in the T2 or the T1 direction. However, when the simu- 
lator was mounted on isolators and the rotor was inert, the full test objective 
could be safely accommodated in any direction. 

Transverse Shock Testing 

Shock testing in the transverse direction required that transducers be more 
closely monitored than was the case for the axial direction, as may be seen by 
comparing Figure 6^ with Figure 36. The additional data presented in Figure 64 
results from two factors: 1) orbital motion between the rotor and casing 

is c/ccited, as indicated in oscilloscope photos b, c and d; and 2) the journal 
bearing pad-to-rotor motions are significantly induced by transverse shock, as 
IS seen in photos a, e, f, and h. 
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The remaining photo, g, presents an accelerometer reading on the input side of 
the isolators (lower trace) and one on the output side (upper trace). These 
simultaneous traces, taken collectively, indicate that the effect of the isola- 
tors was to reduce the amplitude of transversely applied shock by a factor of 
two. A secondary effect was to spread the shock pulse over a longer time, which 
indicates that much of the shock energy was transferred rather than immediately 
dissipated by the isolators. Figure 64 also gives a good indication of the type 
of data that was displayed on the screens of the oscilloscopes shown in the 
**Data-Displays” section of Figure 10. Although the nature of the responses were 
completely different for the axial and transverse excitations, the same displays 
were monitored. The only circumstance in which basically different monitoring 
was required was when the nonrotating tests were conducted, as these required 
paying considerable attention to strain-gage data. 
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PHYSICAL EFFECTS ON BEARING SURFACES AND COMPONENTS 


At the start of this study, the turbocompressor simulator, which had previously 
been used for training purposes, was disassembled and overhauled by replacing 
worn or defective components. Then the active bearing surfaces were refurbished 
and photographed for later comparisons ("before-and-after”) . The condition of 
these surfaces, at various intervals in the test program, is fully documented in 
Appendix C. 

After the machine was reassembled, adjusted, and brought up to operating speed 
(39,000 rpm) , testing was begun by applying vibration in the axial direction 
with the simulator rigidly mounted. Toward the end of these tests, one of the 
pad supports in the turbine-end journal bearing failed. The location of this 
support, referred to as a pad flexure, is shown schematically in Figure 2 
adjacent to probe PF22 in the upper-right-hand comer of Section BB. The 
flexure and its associated pad were in resonance at 1365 Hz, the frequency of 
applied vibration, when the failure occurred. 

Figure D1 of Appendix D presents two views of the failed turbine-end journal- 
bearing flexure. The failure was detected from the instrumentation signals which 
describe the position of the journal relative to the casing, as well as from the 
signals which describe the film thickness between each pad and the journal. 

These signals indicated a sudden displacement of the journal, followed by 
repeated contacting between the pads and journal. The simulator was shut down 
without difficulty. The broken flexure was located, removed, and replaced with 
a spare during a partial disassembly of the simulator. Visual inspection of the 
journal at this time did not reveal any. debris nor any damage to the journal or 
pad as a result of the failure; the simulator was reassembled and testing was 
resumed. 

Appendix D contains a complete report on the failed flexure. The report concludes 
that the failure resulted from the propagation of a crack which was formed during 
the quenching operation of the heat treatment process used when fabricating the 
part. Since the failure occurred when a resonance of this particular component 
was being explored, it seems reasonable to conclude that the amplitude of the 
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dynamic deflection at resonance resulted In a stress level of sufficient magnitude 
to permit accelerated crack propagation. It should be noted that materials, 
design techniques, and heat treatment techniques exist that permit the avoidance 
of this particular type of failure In future flexure designs. The following 
paragraphs document the sequence In which the vibration and shock testing was 
performed, and the results of the periodic checks and Inspections of the bearing 
parts made during the tests. * 

The first series of vibration tests, as Identified below, were all conducted with 
the simulator rtinning at 39,000 rpm: 

1. Axial vibration, rigid mounts 

2. Axial vibration, isolation mo\mts 

* 

3. Transverse (Tl) vibration, rigid mounts 

4. Transverse (Tl) vibration, isolation mounts. 

At the conclusion of the above tests, a complete disassembly and inspection of 
the bearing parts was performed. The results of this inspection are documented 
in Appendix C. The bearing surfaces appeared to have collected some dust and other 
fine debris as the result of moisture and incomplete filtering of the air used 
to drive the rotor. The surfaces were cleaned with soap and water, then with 
alcohol. The simulator was then reassembled and adjusted in preparation for 
the nonrotating (rotor inert) tests. 

The following series of vibration tests were next conducted with the rotor in a 
nonrotating condition: 


1. 

Transverse 

(Tl) 

vibration, 

isolation 

mounts 

2. 

Transverse 

(Tl) 

vibration. 

rigid mounts 

3. 

Transverse 

(T2) 

vibration. 

isolation 

mounts 

4. 

Transverse 

(T2) 

vibration. 

rigid mounts 

5. 

Axial vibration 

, isolation 

mounts 


6. 

Axial vibration 

, rigid mounts . 



*Recall that (Tl) refers to the transverse direction with the rotor supported 
by the two adjacent "rigid" pads, while the (T2) direction is also transverse 
but with the rotor supported by one "rigid" and one adjacent "compliant" pad. 
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At the conclusion of the transverse testing indicated above, one journal bearing 
pad was removed and inspected because a small, but noticeable, increase in the 
housing-to-rotor clearance had occurred. The pad surface was lightly covered 
with a small amount of chrome oxide (coating) dust and the pad pivot had slight 
wear marks. These indications were not considered significant, so the pad was 
replaced, adjusted, and testing continued. 

Upon completion of all the above indicated nonrotating tests, a complete strip- 
down and inspection of bearing parts was performed, including photographing the 
bearing surfaces. This stripdown is also documented in Appendix C. It includes 
photos of the pivot surfaces. Parts were thoroughly cleaned, and the machine 
was reassembled for resumption of tests. The vibration testing was concluded 
by conducting the following test series at a rotor speed of 39,000 rpm: 

1. Transverse (T2) vibration, isolation mounts 

2. Transverse (T2) vibration, rigid mounts. 

At the conclusion of these tests, another disassembly, inspection, and reassembly 
(without cleaning) was accomplished before shock testing. All pad pivots were 
replaced at this point in order that the shock effects could be assessed 
independent of the vibration history of these components. The pads, on the 
other hand, were left intact because of the expense involved in their replacement. 

As was mentioned earlier, the shock-test objective was to apply a 20-g, 10-milli- 
second, half-sine pulse to the simulator in each of three directions, using both 
rigid mounts and isolators, and with the rotor running and then inert. Early 
in the shock-test program, however, an error in the setting of a pressure 
regulating valve on the shock-test machine resulted in four axial impacts being 
imposed at a 60-g level, instead of the specified 20-g level. The shaft was 
rotating during these 60-g impacts and the simulator was mounted on isolators. 

No change in bearing performance was noted as a result of this high level of 
shock. During an inspection following the last of these high impacts, it was 
noted that the dynamic deflection of the simulator on the isolators was 
approximately one inch. Also, a fitting was nearly sheared off from the piping 
which provided external pressure to the thrust bearing for startup purposes. 

Since this fitting was one of two which were connected in parallel, it was 
removed, the hole was plugged, and testing was resumed. Clearance checks were 



performed and, as has been stated, no damage or deterioration was detected from 
the excessively high shock pulses. 

A total of 73 impacts were imposed on the simulator before the shock testing was 
completed and the machine was stripped down for the final inspection. Results 
of the final inspection are shown in Appendix C. Other than a few minor wear 
areas and rub spots, the bearing surfaces were unaffected by all of the testing, 
even though repeated rotor-bearing contacts occurred during the shock testing. 
Furthermore, there was no apparent degradation in performance as a result of the 
engine shock- and-vibration test program. A breakdown of the 73 impacts imposed 
on the simulator is as follows: 

1. 42 impacts in the axial direction 

2. 18 impacts in the T1 transverse direction 

3. 13 impacts in the T2 transverse direction. 

With respect to the nonrotating vibration tests (these being the only vibration 
tests during which bearing contact occurred), a total of 27 minutes of test time 
was accumulated; approximately 17 minutes in the transverse directions and 
10 minutes in the axial direction. 

Throughout the test program, the journal bearings were operated hydrodynamically 
during simulator startup and shutdown. This means that there was dry sliding 
of the journal bearing surfaces during the intervals of initial startup and final 
shutdown (i.e., during the intervals when shaft speed was not sufficient to 
generate a hydrodynamic gas film). In excess of 40 hydrodynamic startups were 
accumulated with the simulator horizontal (i.e., with shaft weight supported 
by the journal bearings). Many more startup cycles were accumulated with the 
simulator vertical. As has been demonstrated in previous material evaluation 
studies, the plasma-sprayed chrome oxide coating used on all of the bearing 
surfaces survived the start-stop sliding and the vibration- and shock- induced 
contacts in excellent condition. 
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CONCLUSIONS 


The objective of this study was to obtain an initial assessment of the effects 
of environmental shock and vibration on the mechanical performance of gas- lub- 
ricated Brayton Cycle turbomachinery for space-power applications. A three- 
fold approach was used to make this assessment: 

1. Experimental observation of the effects of shock and vibration on 
the performance of a gas- lubricated rotor-bearing system of known 
characteristics ; 

2. Periodic inspection of the bearing parts and bearing surfaces during 
the course of the shock and vibration test program, to determine the 
extent of rotor-bearing system damage and/or wear resulting from the 
shock- and vibration-induced dynamic stresses and bearing contacts; 

3. Development of analytical methods for predicting shock and vibration 
response of gas- lubricated rotor-bearing systems, and validation of 
the analytical methods via comparison of predicted and measured re- 
sponses for the experimental shock and vibration test system. 

-Observations and conclusions resulting from this threefold assessment are pre- 
sented in the following paragraphs. 

Experimental Observations 

Observations of the rotor-bearing system response during the sinusoidal vibra- 
tion tests with shaft rotation and rigid mounts, revealed a strong susceptibil- 
ity for bearing contacts over a broad frequency range encompassing the three 
rigid-body critical speeds of the rotor. These critical speeds, all occurring 
below the 650 Hz rotational speed of the simulator, were as follows: 

1. An axial critical speed of the rotor on the thrust bearing at about 
130 Hz which produced large dynamic amplitudes in the thrust bearing 
film; 
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2, Transverse critical speeds of the rotor on the journal bearings at 166 
and 225 Hz which produced large orbital motions of the rotor. 

At excitation frequencies above the rotor speed, there were likewise broad fre- 
quency ranges where the susceptibility for bearing contacts was high. These 
frequency regions were associated with the individual journal bearing pads, with 
the rotor casing, with the simulator test stand, and with the vibration test 
machine itself. 

The use of vibration isolation mounts greatly reduced the susceptibility for bear- 
ing contacts. At excitation frequencies above rotor speed, the vibration test 
objectives were fully achieved (for the conditions of shaft rotation and contact- 
free operation of the bearings) when isolators were used. However, except for 
the T2 transverse direction, the objective input vibration amplitudes still could 
not be achieved in narrow frequency bands around the three rigid-body critical 
speeds of the rotor. Furthermore, the isolators aggravated the rotor-bearing 
system response in the frequency region below 20 Hz because of the low-frequency 
casing resonance associated with the isolators. 

Application of 20-g, 10-millisecond shock pulses to the running simulator caused 
rotor-bearing contacts for all three directions of the applied shock loads. The 
use of isolators did not prevent the contacts, but did change the nature of the 
shock response. For example, axially applied shock without the isolators caused 
the initial bearing contact to occur between the thrust runner and the forward 
thrust plate, whereas the initial contact occurred between the runner and the re- 
verse thrust plate when isolators were used. 

A subtle, but noticeable and significant, effect was observed during axial shock 
testing without isolation mounts. Initially the runner would rapidly approach 
the thrust plate as the gas film collapsed. Then, quite suddenly, the velocity 
appreciably decreased just before contact. This braking action (which has been 
identified by the analysis results to be due to gas- film damping) greatly re- 
duces the thrust bearing susceptibility to shock damage. 
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Bearing System Damage and Wear 


Early in the vibration test program, there was a failure of one of the pad- 
support flexures in the turbine-end journal bearing. This failure occurred 
while purposely exciting the resonant frequency of the flexure. The cause of 
the failure is fully documented in Appendix D. Improved design and manufactur- 
ing procedures, as well as improved materials, are available which can eliminate 
this type of flexure failure. 

Aside from the flexure failure, there was no significant effect of the shock 
and vibration testing on the condition of the bearing surfaces and bearing 
parts. Excluding a few minor wear areas and polished rub spots, the surfaces 
were essentially unaffected by the testing even though several rotor-bearing 
contacts occurred during each of the 73 shock tests, and continuous contacting 
occurred during the nonrotating vibration tests. In addition, **dry** (hydro- 
dynamic) starts and stops of the simulator were performed without difficulty 
throughout the shock and vibration test program. 

It should be recognized that the ability of gas bearings to operate reliably 
under repeated contact and sliding conditions is strongly influenced by the 
surfacing materials used on the bearings. This has been demonstrated by tests 
reported in References 1, 2, 9 and 10. The bearing parts used in the simulator 
for the shock and vibration tests were all surfaced with plasma- sprayed chrome 
oxide. Chrome oxide is the optimum surfacing material thus far identified for 
gas bearings for temperatures up to 600°F. 

Analytical Study of Shock and Vibration Response 

Separate analyses were performed for the axial and the transverse responses. 

The axial analysis dealt with response of the rotor/ thrust-bear ing system, 
both with and without isolation of the rotor casing. Both axial shock and 
axial vibration excitation were considered. The transverse analysis dealt 
with response of the rotor/ j ournal-bear ing system, both with and without 
isolators. Only transverse vibration excitation was treated. 

The adequacy of the shock and vibration analyses were assessed by comparing 
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measured and calculated responses for the rt>tor-bearing system simulator. A 
single computer program was developed and used to perform all of the axial re- 
sponse calculations with the simulator considered either as a three-mass system 
or a one-mass system. Using the three-mass model for the rigidly mounted simu- 
lator, the first axial critical speed, and the corresponding response amplitudes, 
were accurately predicted. Accurate predictions of the axial critical speed 
could not be obtained with the single-mass model. 

The axial shock response, in terms of amplitude/time history, total time to con- 
tact, and the nature of the braking effect which occurs just before rotor bottom- 
ing, was accurately calculated using the three-mass model of the simulator on 
rigid mounts. Reasonably good agreement with the measured responses could also 
be obtained with the single-mass model of the rigidly mounted simulator. How- 
ever, for the case of the simulator mounted on isolators (for which the three- 
mass model must be used) , only qualitiative agreement was obtained between the 
calculated and measured axial shock response. The poor quantitative correlation 
was found to be due to the fact that the isolators were modeled as linear ele- 
ments. Subsequent measurements showed the isolator load characteristics to be 
highly nonlinear. This source of error was verified by bypassing the isolators 
and using the measured casing acceleration as the input shock function for the 
computer program. The resulting response calculation was in close agreement with 
the measured response. 

Simply expressed, the fundamental conclusion relating to the axial response cal- 
culations is that the masses, stiffnesses, and nonlinear-gas-film properties 
were all adequately modeled; the shock-and-vibration isolators were not. 

The most significant aspect of the analytical approach used to predict the axial 
response is the application of simple, frequency-independent, nonlinear functions 
to approximate the theoretical static load capacity and small-amplitude (sinusoi- 
dal) viscous damping characteristics of the thrust-bearing gas film. This 
simplified approach eliminates having to perform a step-by-step solution of the 
Reynolds equation simultaneously with the ordinary differential equations of the 
mechanical elements. The simplified approach can be solved much more rapidly 
than can the rigorous system equations. 
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Calculations of transverse response were made only for the case of steady-state 
vibration applied to the simulator through vibration isolators. The signifi- 
cant frequencies calculated were: 

lo All journal bearing pad resonances 

2, First resonant bending frequency of the rotor casing 

3. First and second transverse critical speeds of the rotor. 

The only pad resonant frequency that was experimentally identified was within 
seven percent of the calculated value„ The casing resonance was not clearly 
identified experimentally, but there were indications that it occurred in the 
vicinity of the calculated resonance. 

The first and second transverse critical speeds of the rotor-bearing system 
were calculated to within one percent of the measured values. However, ampli- 
tudes and orientation of the rotor- to-casing displacements were considerably 
in error. This could be due to several reasons, listed in the following 
order of likelihood: 

1, *'Large amplitude" nonlinear stiffness and damping effects in the 
journal bearings which were not represented in the analysis; 

2. The effect of damping from the thrust bearing and from the shaft 
labyrinth seals which was also neglected in the calculations; 

3o The effect of "cross talk" (a small amount of horizontal vibration) 
which was present in the vibration table, but not accounted-f or in 
the analysis. 

Summary of Significant Results 

The significant results of this investigation are summarized below. From a 
quantitative standpoint, it must be remembered that the results apply to one 
specific rotor-bearing system. From a phenomenological standpoint, however., 
the system response characteristics are indicative of those which would be 
exhibited by similarly configured gas- lubricated rotor-bearing systems. 
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Shock Tests 


1. Gas-lubricated rotor-bearing systems of the type tested can satisfac- 
torily survive at least a limited number of externally imposed 10- 
millisecond, 20-g shock impulses, both with and without shaft rotation, 
and with or without isolation mounts. For the system tested, momentary 
contacts between the rotating and stationary bearing surfaces were 
observed under all test conditions. However, there was no surface dam- 
age or degradation of bearing performance as a result of these contacts. 

Sinusoidal Vibration Tests ~ Shutdown Condition (Nonrotating Shaft) 

1. Under shutdown (nonrotating) conditions, gas-bearing-supported rotor 
systems of the type tested can satisfactorily survive, at least for 

a limited period of time, axial imposition of the specified sinusoidal 
vibration conditions, either with or without isolation mounts. 

2. For transversely imposed vibration under shutdown (nonrotating) con- 
ditions, vibration isolation may be required to survive the specified 
sinusoidal excitation conditions. During the nonrotating transverse 
vibration tests without isolators, the objective input excitation levels 
could not be achieved in the frequency range from 190 to 235 Hz, this 
being the region of one of the critical speeds (resonant frequencies) of 
the nonrotating shaft. The limiting factor in this frequency range was 
deflection (overstressing) of the flexures used to support the individual 
journal bearing pads. With isolators installed, excitation of the flex- 
ures was greatly reduced and the objective input vibration levels could 
be safely imposed in the two transverse directions. 

3. A total of 27 minutes of vibration testing was accumulated under non- 
rotating conditions. No damage to the bearing surfaces nor degradation 
of bearing performance was detected as a result of the nonrotating tests. 

Sinusoidal Vibration Tests ~~ Normal Operation (Rotating Shaft) 


1. Under normal rotating conditions, gas-lubricated rotor-bearing systems 
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of the type tested may have to be vibration isolated to survive the 
specified sinusoidal vibration conditions. During the shaft rotation 
tests without isolation, rotor- to-bearing contacts became imminent « 
over broad regions of the frequency spectrum, at input vibration 
levels considerably below the test objectives. However, with isola- 
tors installed, the objective vibration input levels were achievedg 
without bearing contacts, over most of the frequency range. Only in 
the vicinity of the isolator resonant frequency (approximately 12 Hz), 
and in the vicinity of the rotor critical speeds (130, 166, and 225 
Hz ) 3 did bearing contact become imminent at less than the objective 
vibration input levels. 

2. Input vibration levels were carefully controlled during the rotating 
vibration tests to prevent the occurrence of bearing contacts. Con- 
sequently, it is not known whether bearing performance would have 
been degraded under vibration-induced contact conditions. 

Calculated Shock and Vibration Res p onse 

1. Accurate predictions of the axial shock and vibration response of gas- 
bearing machinery can be obtained using a three-degree-of freedom 
analytical model of the rotor-bearing system with nonlinear representa- 
tion of the load and damping characteristics of the thrust bearing gas 
film. An approximate, simplified method of representing the gas-film 
nonlinearities can be used which yields both accurate and economic 

(in terms of computer cost) solutions. The only drawback of the herein 
presented axial response analysis is the modeling of the vibration 
isolators in terms of linear stiffness and damping coefficients. Under 
shock conditions (i.e., large displacement conditions) the isolators 
are quite nonlinear and must be so modeled. 

2. Accurate predictions of the transverse resonant frequencies of the 
rotor-bearing test system were obtained using the finite-element 
analysis model of the rotor and rotor casingj together with the linear- 
ized Journal bearing stiffness and damping coefficients. There was^ 
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however, considerable discrepancy between the calculated and measured 
amplitudes of vibration, and of the rot or- to-casing orbit shapes. Cal- 
culated amplitudes were considerably higher than the measured ampli- 
tudes. The discrepancy in calculated amplitudes may be due to "large 
amplitude" nonlinear stiffness and damping effects in the journal 
bearings which were not represented in the analysis. There was also a 
small amount of horizontal vibration (cross talk) in the vibration 
table which was neglected in the calculations. Finally, there may have 
been significant damping from the thrust bearing and from the shaft 
labyrinth seals. These sources of damping were also neglected in the 
analysis. 

3. The analysis results clearly show the beneficial effects of gas-bearing 
squeeze- film damping on minimizing the severity of shock-induced bear- 
ing contacts. Gas-bearing systems which must operate under shock 
conditions should be optimized to take maximum advantage of this 
effect . 

The above listed results indicate that 10-millisecond, 20-g shock pulses should 
not be a problem for gas- lubricated Brayton Cycle space- power turbomachinery. 
Sinusoidal vibration could be a problem, particularly if imposed during shaft 
rotation. The use of vibration isolators appears to be a solution to the vibra- 
tion problems during shutdown (nonrotating) conditions. During rotation, the 
use of isolators greatly reduces, but may not eliminate, the possibility of vibra- 
tion-induced bearing contacts. Further testing will be required to identify if 
such contacts would cause bearing damage or degradation of bearing performance. 
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RECOMMENDATIONS 


This report describes the results of work performed under Tasks I and II of 
NASA Contract NASw-1713. Additional work is currently being performed under 
Tasks III and IV of the contract as follows: 

Task III - Extension of the axial response computer program to permit 
nonlinear modeling of isolation mounts; Development of 
an axial response computer program for random vibration 
excitation. 

Task IV - Experimental test and evaluation of the Brayton Rotating 

Unit (BRU) simulator under vibration and shock conditions, 
including random vibration. 

As a result of the above on-going contract tasks it is somewhat premature to 
formulate a list of recommendations until these tasks are completed. There is, 
however, one important question identified by the herein reported work which 
will have to be investigated at some point. 

During all of the vibration tests performed while the simulator shaft was 
rotating, the input vibration amplitudes were always controlled to maintain 
contact-free operation of the simulator bearings. Because of this intentional 
desire to prevent bearing contacts, there were numerous instances when the 
objective input vibration amplitudes were not attained because of the imminence 
of contacts. Figures 33, 58 and 60 show that without isolation, there were broad 
frequency ranges over which the objective input amplitudes were not attained. 
Figures 34 and 59 show that even with isolators, there were a few narrow 
frequency ranges where, again, the objective input levels were not attained. 

An important question naturally arises from these results. If the vibration 
tests had been conducted at the objective input levels over the total 5- to 
2,000-Hz frequency range, there undoubtedly would have been bearing contacts 
within certain frequency regions. The question is, would there have been 
bearing damage, or degradation of bearing performance, as a results of the 



vibration- induced contacts? It appears that the effects of vibration- induced 
contacts during shaft rotation must be experimentally assessed in the near future. 
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APPENDIX A 


EQUATIONS OF MOTION FOR AXIAL RESPONSE 

If we represent the displacement of the three masses shown in Figure 15 as , 

and and represent the displacement of the base as x^ in the same inertial 

frame, then the three equations of motion may be written directly from Newton’s 

laws. If we introduce an aerodynamic force, F , acting on m«, (the rotor), and 

a j 

take g to be the Local gravitational acceleration, the equations are: 


"*1^*1 ” ~ ^1^ ■ ' ^l^^l ' ~ “l® 

(Al) 

^-^2 "G(h>h) - " ^1^ ' "*2® 

(A2) 

"'sS " ^a G(h,h) - m^g 

(A3) 

These equations can be expressed in a form more amenable to laboratory 
ments by some redefinition. Resolving the generalized gas force G into 
load component and a damping (dissipative) component, we get: 

measure- 
an elastic 

G(h,h) = W(h) - h B(h) 

(A4) 

Note that initial conditions may be obtained by considering the system 
static equilibrium for which (Al), (A2), and (A3) become 

to be in 

0 = k^Cx^ - x^) - k^(x^ - x^^) - m^g 

(A5) 

0 =-W(h) - ■ ®2® 

(A6) 

0 = + W(h) - m^g 

(A7) 

where x^, x^, x^r and are the values of the coordinates at static equilibrium. 
The equations of motion can also be written in terms of coordinates referenced to 
the static equilibrium position of the system by means of the transformation: 

yi = xi - 

(A8) 
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and similarly, for the coordinates locating and the base, b. We also 

define an initial film clearance ,, h, by 



The actual value of h is determined by the weight of the rotor, the force, F , 
and the thrust-bearing static-load characteristic to be discussed shortly. 
Applying the relationships (A4) through (A9) to Equations (Al), (A2) , and (A3) 
produces a set of equations of motion whose coordinates are referenced to static 
equilibrium: 

”1*4 " 

'”2*^2 ■ ^2 ” ^^2^^2 ' (All) 

” 3^3 ^I^y3 ~ y 2 ■ (y3 ■ y 2 ^ ®^y3 ■ y 2 (A 12 ) 

Note.*,* The expressions in square brackets immediately following W and B are argu- 
ments of the nonlinear functions W and B. For example, the value of B, in gen- 
eral, depends upon the value of the film thickness, h. The value of B determined 
in Equation (A12) is that which occurs when h = [y^ - y^ + h] . 

Equations (AlO) through (A12) were expressed as a set of six, simultaneous, first- 
order, nonlinear differential equations. These equations were programmed on a 
digital computer using a fourth-order, Runge-Kutta, integration technique supple- 
mented by an error control method which continuously adjusted the integration 
step size to maintain solution accuracy within prescribed limits. 

In certain instances, computer time is saved without significant loss in accuracy 
by representing the system with a single-degree-of -freedom model. This is per- 
missible when the thrust -bearing gas film provides substantially all of the com- 
pliance in the system. The analysis of such a simplified system results in a 
single, second-order, nonlinear, differential equation which describes the motion 
of a single lumped mass, m (the rotor), in terms of a displacement coordinate, x, 
in an inertial reference frame: 
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(A13) 


mx = F + W(h) - hB(h) - mg 

3 

As before, the motion is better expressed with respect to the position of static 
equilibrium. Using the same manipulations previously performed we arrive at: 

my - W(h) - W(h) - (y - y^) B(h) (A14) 

where y defines the position of the rotor with respect to static equilibrium. 

Since the equations of motion are written for all three representations of the 
system, (i.e., three-mass rigid mount, three-mass isolation mount, and one-mass 
rigid mount) in terms of general base motion, the computer program was written 
to handle any of the following four types of forcing functions expressed as either 
base acceleration or as base displacement: continuous vibration, half-sine shock, 

haversine shock, or arbitrary shock pulse read in as a table of amplitude and 
time values. 

As shall be shown in Appendix B, the computer program user can either elect to 
represent his system by Figure 15 or 16, which results in simultaneous solution 
of (AlO), (All), and (A12); or he may elect to use the representation of Figure 
14 which requires only the solution of (A14) . 


57 



APPENDIX B 


AXIAL RESPONSE COMPUTER PROGRAM 


This appendix presents a complete listing of the computer program, (written in 
the FORTRAN IV language), used to calculate axial response of the system as re- 
presented by Figures 14, 15, or 16. It also provides detailed instructions for 
preparing the input data cards, along with sample input and sample printed output 
Some of the typical graphical output has already been presented in the main text. 
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PROGRAM INPUT 


Card 1 Read NCASE, TITLE 
FORMAT (15, 75H) 

Set NCASE =* Any nonzero character to sequentially identify the case being 
computed (a blank here will stop the program) 

Title s Any message to describe case being computed. If computer 

graphs are requested, this message will appear at top of graph. 

Card 2 Read MODEL, LOADF, IDAMP, IFOFT, ICC, MOTION 
FORMAT (615) 

These are options as follows: 

Physical Model 

MODEL = 1. Complete system on isolators 

MODEL = 2. Complete system on rigid mounts with flexible lower casing. 

MODEL = 3. Simple model, gas film with nonlinear stiffness and damping; rigid 
casing and mounts. 

Gas film properties 
LOADF = 4 Exponential load 

LOADF = 5 Power function load 

IDAMP = 6 Exponential damping 

IDAMP = 7 Power function damping 

Type of Excitation Function 
IFOFT = 8 Continuous sine function 
IFOFT = 9 Shock function half sine 
IFOFT = 10 Shock function haversine 
IFOFT = 11 Shock function (table) 

Calcomp desired? 

ICC = 12 Yes 

ICC = 13 No 
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Is forcing function described in terms of displacement or acceleration? 

MOTION =0 if forcing function input represents displacement in inches. 

MOTION =1 if forcing function input represents acceleration in g*s. 

Card 3 - data is dependent upon "MODEL" option as follows: 

If MODEL = 3 
Read M3 

FORMAT (20X, E10.4) 

M3 = weight of rotor (lb) 

If MODEL = 2 

Read Ml, M2, M3, K1 , Cl, K2 
FORMAT (6E10.4) 

Ml = Weight of "lumped" mass assigned to lower region of casing (lb) 

M2 = Weight of thrust bearing stator (lb) 

M3 = Weight of' rotor (lb) 

K1 = Value of "lumped" stiffness assigned to lower region of casing (Ib/in.) 

Cl = 0.0 

K2 = Stiffness of thrust bearing flexure (Ib/in) 

If MODEL = 1 

Read Ml, M2, M3, K1 , Cl, K2 
FORMAT (6E10.4) 

Ml = Weight of casing (lb) 

M2 = Weight of thrust bearing stator (lb) 

M3 = Weight of rotor (lb) 

K1 = Total stiffness of isolators (Ib/in.) 

Cl = Total damping of isolators (lb sec/in;) 

K2 = Stiffness of thrust bearing flexure (Ib/in.) 

Card 4 - data is dependent upon "LOADF" and "IDAMP" options as follows: 

READ ALQAD, BLOAD, CLOAD, DLQAD, ADAMP, BDAMP, CDAMP, DDAMP 
FORMAT (8E10.4) 

If LOADF = 4 the expression used to represent the nonlinear load characteristic 
of the gas film is 
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W(h) = 
where; 

ALOAD = A 
BLQAD = B 
CLQAD = 0. 

DLQAD = 0. 

If LOADF = 5 the expression used to represent the nonlinear load characteristic 
is 

W(h) = Ch® 
where : 

ALOAD = 0. 

BLQAD = 0. 

CLOAD = C 
DLQAD = D 

I f IDAMP ^ 6 the gas film damping is given by 
B(h) = Ae®*^ 
where : 

ADAMP = A 
BDAMP = B 
CDAMP = 0, 

DDAMP = 0, 

If IDAMP «= 7 the gas film damping is given by 
B(h) = Ch° 
where : 

ADAMP = 0 
BDAMP = 0 
CDAMP ■= C 
DDAMP = D 

Card 5 - Read GL, FA, ELIM 
FORMAT (3E10.4) 

GL = Local gravitation acceleration in "g's” 

(GL = 1.0 on earth's surface.) 
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FA = Net aerodynamically induced axial force on rotor (positive, if force is up- 
ward) 

ELIM = Nonlinear load and damping functions are held constant at film thicknesses 
less than ”ELIM" inches 

Card 6 - Data is dependent upon options chosen for IFOFT and MOTION. 

If IFOFT = 8 Read AWPLD, FREQ 
FORMAT (2E10.4) 

AMPLD = amplitude of sinusoidal forcing function. 

If MOTION = 0, AMPrLD = displacement in inches 

If MOTION = 1, AMPLD = acceleration in g*s. 

FREQ = frequency of forcing function in Hertz (cps) 

If IFOFT = 9 Read PEAKS, PDURS 
FORMAT (2E10.4) 

PEAKS = peak value of half sine shock pulse 

If MOTION = 0, PEAKS = displacement in inches 

If MOTION = 1, PEAKS = acceleration in g*s 

PDURS = duration of half sine shock pulse (seconds) 

If IFOFT ^ 10 Read PEAKH, PDURH 
FORMAT (2E10.4) 

TEAKH = peak value of haversine (shock machine half sine) pulse 

If MOTION = 0, PEAKH = displacement in inches 

If MOTION = 1, PEAKH = acceleration in g*s 

PDURH = duration of haversine shock pulse (seconds) 

If IFOFT = 11 Read NINC, TINCP 
FORMAT (15, E10.4) 

NINC = number of evenly spaced intervals in pulse 
TINCP = time interval of each of the above increments 

This data card must be followed by as many data cards as are necessary to repre- 
sent the shock pulse. Call these cards 6A, 6B, etc. 

Cards 6A, 6B, etc: Read PUVAL 
FORMAT (8E10.4) 
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PUVAL values of the shock pulse at each increment specified above. Values must 
be contiguous starting with the first increment. 

If MOTION = 0, PUVAL displacement in inches. 

If MOTION = 1, PUVAL - acceleration in g*s. 

Note, if more than one increment size is desired, a different case may be run for 
each step size by using the conditions at the change point as restart values. 

Card 7 Read IPSTEP, M, MAXING, NCYCLS, TPRINT, TMAX, H, E, HM, RERR 
FORMAT (415, 6E10.4) 

This card controls the precision and duration of the internal computation and 
output . 

IPSTEP e ratio of printed to plotted output points. In most cases IPSTEP = 1. 

M ^ number of integrations performed internally before step size is auto- 

matically adjusted. Usually M «= 2. This will result in the best 
trade off between precision and economy. 

MAXING = the maximum number of internal integrations before program shuts it- 
self off. This is a safety control to prevent exorbitant computer 
charges due to unforeseeable misbehavior of the solution. 

Recommended MAXING = 1000 . 

NCYCLS = the number of cycles for which the solution will be calculated beyond 
the first minimum and maximum after the shock pulse ceases to act, 
NCYCLS has meaning only for shock pulses. Usually NCYCLS = 0. 

TPRINT s= output control which establishes the step size of printed and plotted 
output (seconds). May be any convenient interval into which we 
want to divide the time range of the solution, generally a few hun- 
dred points or so. Example: If we expect the range of interest to 
be between 0 and 2 seconds and we think 200 solution points are ade- 
quate to describe the system behavior in detail, we select 

TPRINT = ' ^200^^" = 0.01 

and will get results at 10 millis,econds , 20 milliseconds, etc. 
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IMiX = maximum time value (seconds for which solutions will be computed. 

This is another safety feature designed to prevent excessive compu- 
tational charges resulting from some unforeseeable characteristic 
of the solution. Normally this value should be twice as large as we 
expect the solution to run. 

H •- time increment (seconds for which to begin internal computations. 

Recommended value is 1/4 of TPRINT. 

E =--■ maximum absolute error permitted in internal computation. 

Recommended value is .01 

HM ~ minimum step size (seconds) for internal calculation. This parameter 

is another safety measure which keeps the step size from approaching 
zero in order to strive for greater accuracy. Recommended minimum 
value: l/lOO of TPRINT 

RERR “ comparable to E, except that RERR represents relative error whereas, 

E represents absolute error . Recommended value 0.01. 

Card 8 Initial value or RESTART CARDS 

See analysis section for definition of these coordinates 
READ T1NI.T, YlIN, Y2IN. Y3IN, VlIN, V2IN, V3IN 
EORMAT (7E10.4) 

TIMI initial or restart time (seconds) 
y I IN initial value of Y^ (inches) 

Y21N “ initial value of Y.^ (inches) 

Y3TN - initial value of Y^ (inches) 

VllN = initial value of (inches/sec) 

V^21N = initial value of (inches/sec) 

V3IN =■ initial value of (inches/sec) 

Rep-.-.u all above cards for each case. Insert a blank card as the Last data card. 
Ihis will stop t!'r.^ program gracefully. 
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S :: POWER FUNCTION LOAD CHARACTERISTIC 


damping CHARACTERISTIC 


& = EXPONENTIAL DAMPING COEFFICIENT 


7 r POWER FUNCTION DAMPING COEFFICIENT 


■fTESCRIPTION OF GAS FILM 


THE EXPONENTIAL LOAD CHARACTERISTIC IS OF THE FORM 
LOAD < LB ) ~ A* EX PI B* H) 


WHE PE 



1 

i 


H = FILH THICK 1C SS ( INCHES I 

A : . 1 05 n* 03 I lb ) 

B - -.2 55 5* DM 1/ IN I 


T HE _P OK EP_F UNCTION OA HPJN G CO EF FI Cl EN T IS OF THE FORM _ 

B (lb*se’c7in ) = C*H**0“ 

WHEPF 

H = FILM THICHNESSS I INCHES I 

• * 31 U -0 3 LB -SEC /I W 

"” 0 - V -•UBB+Ol rf/TNl” 

NON -L INEAp load A ND ~D AM pFnG F UNCTTo A RE~iTel D CO NS TA NT aT^ FI LH T HI CK NE SS ES L £S S TH AN .sllOO-O** INCHES. 


«« ** •• «• •• «• •• ** •« •« •« mm mm mm mm mm mm mm mm mm mm mm mm mm mm mm mm mm mmmm mmmm mm mm mm mm mm mm mm mm mm mm mm mm mm mm mmmmmmm 

EQUILIBRIUM film THICKNESS IS .90 126 HRS . 



ASSUMING STEADY-STATE AERODYNAMIC 

FORCES ON THE ROTOR OF .00 00 0 

L0 


♦ « < 

and a local gravitational field OF 

1. 00000 TIMES 

mm •* mm mm mm mm mm mm mm < 

S TA NO AP 0 




parameters of THE SYSTEM ON RIGID MOUNTS ARE AS FOLLOWS 





EFFECTIVE MASS 

OF lower housing 

EFFECTIVE HASS 
OF STATOR ♦ flexure 

MASS OF ROTOR 

OF 

STIFFNESS 
LOWER HOUSING 

flex (HE 
STIFFNESS 



( MA SS 1 ) • ( lb ) 
1 . «iD DO 

I HASS 21 *C lb 1 
.<♦0 00 

CMASS3 It (LB) 

1 0.500 

(LB/ IN ) 

. 58 00 ♦O 5 

(LB/ IN 1 
.G 7on«^n6 


«• i 

TYPE OF EXCITATION 

function 






8 

□ steady sinusoidal input 





9 

= half sinf shock 




.. 



■ 



10 r HAVERSINE SHOCK 








11 r arbitrary SHOCK (TABLE) 


* 

• ••««* 








IS CALCHSrp^ OUTPUl ULSIQEO- 


12 = YES 


C 7 N 



ON 1 3 = NO 

00 •« «• «* •« *« ^ •••* *^*4 •*♦•»• ••«• «« •« •* ♦* •• #• *♦«* *• *• •« *« «• «• •* •« «« «•*«•«*«•«•** 


INPUT HOTTON IS HAVERSTNE SHOCK PULSE 


PULSE DURATION IS -0 1*16 SECONDS 


PEAK amplitude IS 15 .700 0 G* S ACCELERATION 




program CONTROL PARAMETERS 


DESIRED TERMINATION 

TIME VALUE WILL NOT EXCEED -0 30 0 SECONDS 

AS A SAFETY MEASURE VE MILL NOT PERMIT MORE THAN 80 00 STEPS 

CALCULATION WILL CONTINUE I CY a ES AFTER PULSE DURATION 

STEP SIZE 

initial STEP SI 7E FOR CALCULATION PURPOSES IS .1 000-03 SECONDS 

values are PRINTED, out EVERY . 5000-03 SECONDS 

EVERY ITH POINT IS PLOTTED TF CALCOMP IS USED 

ERROR CRITERION F or INTEGRATION TECHNIQUE -2 00 0-02 

MINIMUM STEP-SIZE r, .1 00 00 -05 SECONDS 

NO- OF INTEGRATIONS BEFORE STEP-SIZE IS DOUBlEDr 2 

normalized ERROR LIMIT FOR INTEGRATION TECHNIQUE .0100 

INITIAL OR RESTART CONDITIONS 

INITIAL time - -00 00 SECONDS 

DISPLACEMENT OF MASS 1 r . 00 00 INCHES 

VELOCITY OF MASS I = .00 00 IN /SEC 

displacement of HASS 2 - .00 00 INCHES 

VELOCITY OF MASS 2 r . 00 00 IN/SEC 

displacement of MASS 3 r .00 00 INCHES 

VELOCITY OF HASS 3 = .00 00 IN/SEC 
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uuu 

T 

P 

uuu 

T 



TIIC 

FILM 

DISPLACEMENT 

OTSPLA CEMENT 

displacement 

VEL OC ITY 

VEL OCITY 

VELOCIT Y 

REL AT IVE 

FORCING 


THICKNESS 

OF NASS 1 

OF NASS 2 

OF ROTOR 

OF MASS 1 

OF MASS 2 

OF ROTOR 

DISP. BET. 

FUNCT ION 

IStCI 

fiNCiCSI 

CINCHESI 

C INCHES ) 

Cl NCHES) 

(IN/ SECJ 

(IN/SECI 

1 IN /SEC ) 

MASS? ♦MASSl 

CINCHES* 

•OODS 

-9 012-03 

- 13%%-OS 

• 9%36-07 

. 30 15-0 8 

. 15 11-0 2 

. 11 51-02 

.44 43 -04 

• 40 02-07 

• 14 59 -OS 

•0010 

.8 965-03 

.59%3-05 

.5 13 9-05 

. 40 07-0 6 

. 31 23-0 1 

.27 89-01 

-27 12-02 

• 00 42-06 

• 2324-04 

.001 5 

•86% 8-03 

• %6%0'0% 

.4 22 8-04 

. 58 25 -0 5 

. 15 03-0 0 

. 13 95-00 

.2500-0 1 

.41 26-05 

• 1168-0 3 

.0020 

.7 716-03 

. 17 8ii “0 3 

.165 7-03 

. 36 04 -0 4 

. 40 06-0 0 

.37 56 -00 

.11 25 *0 0 

.12 72-0 4 

>3651-03 

:oo25 ■ 

.^01^-03 

.%6 89-0 3 

-4 38 3-03 

. 14 37-0 3 

. 78 33-0 0 

. 73 51-00 

. 35 38-0 C 

.30 62-0 4 

.8790-0 3 

•0030 

-% 11 8-03 

• 96 98-0 3 

.928 0-03 

. 43 85-0 3 

. 13 42 ♦O 1 

. 12 65+01 

• 88 82-0 0 

.61 83t04 

.1792-0 2 

.0035 

-283^-03 

.18 73-02 

.1 76 6-02 

. 11 08-0 2 

. 22 70 1 

. 21 68 +0 1 

• 18 80 +0 1 

. 10 68-0 3 

.32 55-02 

•00«0 

.1952-03 

- 32 87-02 

.3 136-02 

. 24 30-0 2 

. 37 17*0 1 

. 35 11 +0 1 

. 34 94 +0 1 

.15 04 -0 3 

_.S4 26-D2 

-00% 5 

.1 7^0-03 

.56 63-02 

.5 44 8- 02 

. 47 18 -0 2 

. 60 69 ♦O 1 

• 58 06 +Q 1 

. 57 40 ^0 1 

.2156-03 

• 84 68 -0 2 

.0050 

-1 317-03 

.9302-0 2 

.904 3-02 

.82 74 -02 

. 87 21 1 

. 86 26 +0 1 

. 85 28 «Q 1 

• 25 84-03 

. 12 54 -0 1 _ 

.0055 

.8 279-» 

. 1% 38-0 1 

.1 412-01 

. 13 30-0 1 

. 11 61 *0 2 

- 11 70^02 

• 1159+0 2 

• 26 71 -0 3 

.1778-01 

.0060 

.3 330-0% 

• 20 98-01 

.2 07 3-01 

.19 86-0 1 

. 14 65 2 

. 14 77 +0 2 

• 14 68 +0 2 

.25 29-03 

. 24 32-0 1 

•0065 

-.1 201-0% 

• 29 09-01 

.2 886-01 

. 27 94 -0 1 

. 17 68 +0 2 

. 17 71 +0 2 

.176? +0? 

.23 63-0 3 

• 3226-01 


RUN COl*>tfTC *T THIS POINT. ROTOR mS BOTTOMED. 
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MECHAivilCAL Tt.ChNOLOGY INCORPORATED 
968 AL8ANY-SHAKER ROAD 
LATHAM* NEW YORK 12110 

THIS PROGRAM CALCULATES TmL RESPONSE OF A THRUST BEARING-ROTOR 
system to externally applied shock and VIBRATION. 

IT Treats the system a three-degree-of-freedom vibration model. 

THIS PROGRAM USES ThE FOLLOWING SUBROUTINES - 

INPUT. SIMINI. SOLVE. OIFtO. TLU. FUNCTION. FORMT. FILM. CALCOM. 
RONKU 

PROGRAM PNA57 ( INPUT . OuT PUT . T APtS= INPUT . T APE6=0UT PUT . TAPE22 ) 

COMMON TPRInT.TmAX. IPSTEP.NCYClS.MOT ION. TI TlE ( IS) . YBA. YBInT (200) 
Common NK.NW.NCASE.MoOtL.LoAOF. IDAM r. IFoFT . ICC. 

1 Ml .ME.M3.K1 .Cl .K2.GL.FA.AL0AD.ADAMP.8LUAD.B0AMP. 

2CLOAD.COAMP,L)LOAD.OnAMP.HlNlT » I F M7 . AMPLD . FR£(S . PE AKS . PUURS . PEAKH . 
3PDURH.NINCP.TINCH. TP (200).PUVAL(800) .PuLOOT (200).TINIT.Y1IN.Y21N. 
AY3IN.ViIN.v2IN.V3IN,TlME(bU00).Yl»Y2.Y3.Vl.V2.V3.N.VEC(6) .DVEC(6) . 
SMAxINC.H.hm.M.E.AKO.NE.T. ELIM.INDX. YB.YBUOT. 

60MFGA.WH0.WY3YPH.SY3Y8H 

COMMON FILMT(6U00).Y2MY1(6000) . V APEAK (buOO ) .YBPLUT (bOOO) 

COMMON ibUF(inOO) 

NR = 5 
NW = 6 
IFLAG=1 

5 READ(NR.IOIO) NCASE . ( T 1 TLfc ( I ) . I = 1 » 8) CUt 

IF(NCASE) 10.6.10 

6 IF (IFLAG.EQ.2) CALL PLOT (1 8 . . 0 . . 999 ) 

9 CALL EXl r 

10 WRITE (Nw, 1000 ) 

WRITE (NW. lOlS) 

call soate(date) cut 

WRITE (NW. lOOS) NCASE. ( T ITLE ( I ) » 1=1 . 8). DATE CD< 

CALL INPUT 

IF ( IFLAG.EU.2) GO TO 24 
20 IF ( ICC. EG. 13) GO TO 2S 
IFLAG=2 

CALL PLOTS ( IriUF . 1000.22) 

GO TO 2S 

24 CALL PLOT ( 18. .0. .-3) 

25 CALL SOLVE 
GO TO 5 

1000 FORMAT ( 7H1PN457 . iOX . M hRESPOnSE OF VERTICAL THRUST BEARING TO AXI 
lALLY-APPLIED SHOCK OP VIBRATION/) 

1005 FORMAT ( 15. 7A10.A5.2X.AB) CD 

1010 FORMAT ( IS. 7A10.A5) 

1015 FORMAK// 1 X . 1 1 9 ( 1 H<» ) ) 

END 
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SUBROUTINE INPUT 

THIS SUBROUTINE CALLS SUBROUTINE *FILM* WHICH FINDS EQUILIBRIUM 
FILM thickness. 

COMMON TPRINT»TMAX» IPSTEP ♦NCYCLS* MOTION* TITLE (IS) *YBA*YBINT (200) 
COMMON NR*NW*NCASE»MODEL»LOADF» IDAMP»IF0FT* ICC* 

1 Ml *M2*M3*K1 *C1 *K2*GL*FA*AL0AD*ADAMP*BL0AD*BDAMP* 

2CLOAD*CDAMP*DLOAD*DDAMP*HINIT* IFM7*AMPLD*FREQ*PEAKS*PUURS*PEAKH* 
3P0URH*NINCP*TINCP»TP (200) *PUVAL(200) *PULDO 
4Y3IN*V1IN*V2IN*V3IN*TIME (6000) »Y1*Y2*Y3*V1 
5MAXINC»H*HM*M*E*AKO*NE*T* ELIM 

60MEGA*WH0*WY3Y2H*BY3Y2h 
COMMON FILMT(6000) *Y2MY1 (6000) *V APEAK (6000 
REAL M1*M2*M3«C1*K1*K2 

LOGICAL 01 *02 *03*04 *US*06*07 *08*09*010 
WRITE (NW* IbOO) 

PRINT A ROW OF ASTERISKS 
WRITE (NW*1030) 

WRITE (NW*1020) 

WRITE (NW* 1031 ) 


CHOOSE OPTIONS 

READ (NK* 30) MODEL *LOADF * IDAMP ♦ IFOF T * ICC * MOT ION 
30 F0RMAT(6I5) 


SET LOGICAL VARIABLES AND TEST FOR CORRECTNESS OF INPUT 


01=M0DEL 

• LT 

.1 




02=M00£L 

• or 

.3 




03=L0ADF 

• LT 

.4 




04=L0ADF 

• GT 

.5 




05=IDAMP 

• LT 

.6 




06=IDAMP 

.GT 

.7 




07=IF0FT 

• LT 

• 6 




08=IFOFT 

.GT 

.11 




09=ICC.LT.l 

2 




O10=lCC. 

GT. 

13 




IF (01. OR 

.02 

) WRI 

T£ 

(NW 

*40) 

IF (03. OR 

.04 

) WRI 

TE 

(NW 

♦ 50) 

IF (OS. OR 

.06 

) WR I 

TE 

(NW 

*60) 

IF (07. OR 

.08 

) WRI 

TE 

(NW 

* 70) 

IF (09. OR 

.01 

0) WR 

IT 

£ (NW*80) 


T (200) *TINIT*Y1IN*Y2IN* 
*V2*V3*N*VEC (6) *0VEC(6) * 
*INDX* YB*YBDOT* 

) *YBPLOT (6000) 
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ERROR MESSAGES 






40 FORMAT ( 

35H0 

impkoper 

INPUT 

CARO 

2. COL 5. 

) 

50 FORMAT! 

3bH0 

improper 

INPUT 

CARD 

2.COL10. 

) 

60 FORMAT! 

35H0 

IMPROPER 

INPUT 

CARD 

2.C0L15. 

) 

70 FORMAT! 

35H0 

IMPROPER 

INPUT 

CARO 

2.COL20. 

) 

80 FORMAT! 

35H0 

IMPROPER 

INPUT 

CARD 

2.C0L25. 

) 


IF (01 .0i^.0^?.Ok.OJ. 0H,U<*,OH,O5.OK,O6.O«,O7.OK.O8.OH.O9.OH.O10) 
ICALL f.AlT 

IF (MODEL. Eu, 1) WKlTE('>h»»1040) 

I F‘ (model »E(J • 1 ) i»(K IT E (nW ♦ 1 ObO ) 

IF(MoOEL.EO. 1 ) WkITE (in,w, 1040) 

IF (MOOEL.EQ. 1 ) Wf<I re ('Nw» 1400) 

IF (MODEL. EO, 1 ) wJK I T E ( Nw » 1 ObO ) 

IF (MODEL. EO. 1 ) WKl ?e (Nw* 1400) 

IF (MOOEL.EQ. 1 ) WK 1 1 1 ( ima » 1 0 70 ) 

IF (MODEL. E(j.^) WKITE (Nw, lObO) 

IF (MODEL. £0.2) WRI TE (Nv* , 1^00 ) 

IF (MOOEL.EQ. 2) WK I T E ( Nw , 1 040 ) 

IF (MOOEL.EQ. 2) wK I TE ( Nw ♦ 1 060 ) 

IF (MOOEL.EQ. 2) WRITE (Ww» 1040) 

IF (MOOEL.EQ. 2) WR I T E ( Qw ♦ 140 0 ) 

IF (MOOEL.EQ. 2) WK I T E (Nw « 10 70 ) 

IF (MOOEL.EQ. 2) WK If E ( Nw ♦ 1 400 ) 

IF (MOOEL.EQ. 3) WRI TE (WW» 1050 ) 

IF(M00£L.EQ.3) WR I f E. ( Nw , 1400 ) 

IF (MOOEL.EQ. 3) WR I f £ ( Nw , 1 ObO ) 

IF (MOOEL.EQ. 3) WK 1 TE ( NW » 140 0 ) 

IF(M00EL.EQ.3) WR I TE ( NW ♦ 1 040 ) 

IF (MOOEL.EQ. 3) WRI TE ( NW ♦ 1 0 70 ) 

IF (MOOEL.EQ. 3) WR I TE ( NW » 1 040 ) 


PHYSICAL PARAMETERS I/O 

READ (NR.90 ) Ml »M2»M3.Kl tCl »K2 
90 FORMAT (8E10. 4) 

C CORRECTION TO CHECK FOR M3=0 
IF (M3) 902.900 .902 
C M3 = 0 

900 WRITE (NW.901 ) M3 

call exit 

901 format (24H0 IMPROPER INPUT M3 = .E10.4) 

902 WRITE(Nw.l020) 

WR 1 TE (nw ♦ 1 080 ) 

I F ( LQAOF . Eu .4 ) wR 1 T t (\vjw . 1 040 ) 

IF (LOADF.EQ.h) write (Nw. 1090) 

IF (LOAOF .EQ.*.) wR I T E ( Nw » 1 040 ) 

IF (L0ADF.EQ.4) wRllE (Nw.l400) 

IF (LOADF .Eq.4) wkITF, (Nw.llOO) 

I r ( LOADF .EQ.4) wkIIE (Nw»l‘*00) 


72 



IF(LOADF.EQ.S) WRITE (NW» 1090) 

IF (LOAOF.EQ.S) WR ITE (Nw * 1400 ) 

IF(LOADF.EO.S) WR 1 7E (NW ♦ 1 040 ) 

IF (LOAOF.EO.&) WRITE <NW» 1100) 

IF(LOADF.EQ.S) wR I TE (NW , 1040 ) 

WRITE(NW*1110) 

IF (I0AMP-.E0.6) WR1TE(NW*1040) 

IF(I0AMP.E0.6) WRITE(NW»1120) 

IF (IDAMP.EQ.6) WRI TE (NW» 1040) 

IF(IDAMP.E(J,6) WRITE (NW»1400) 

IF(I0AMP.EQ.6) WRITE (NW»1130) 

IF(I0AMP.E0.6) WRITE (NW*1400) 

IFdOAMP.EQ.T) WRITE (Nw*1120) 

IFdDAMP.EQ.?) WRITE (NW»1400) 

IF(IOAMP,E(J.7) write (NW. 1040) 

IF(I0AMP.E(J.7) write (^^W«1130) 

IFdDAMP.EQ. 7) WRITE (NW« 1040) 

WRITE<NW«1531) 

WRITE(NW»1020) 
write (N w» 1400) 

C READ LOAD AND DAMpjiNG PARAMETERS 

REAO(NR»90) ALOAD*BLOAU,CLOAD»DLOAD* ADAMP,BDAMP*CDAMP*DDAMP 
IF (L0ADF.EQ.4) WRITE (NW« 1230) ALOAD*BLOAD 
IF(LOADF.EQ.S) WRITE (NW* 1235) CLOAD*DLOAD 
IF (IDAMP.EQ.6) wR I 1 E ( Nw ♦ 1245 ) ADAMP.BDAMP 
IFdDAMP.EQ./) WRI IE (NW« 1240) CDAMP.DDAMP 
C READ THE LOCAL GRA v 1 7 A I lONAL ACCELERATION UNITS) AND AERODYNAMIC 
C FORCES 

REAO(NR.90) 6L»FA*EL1M 

C CALCULATE THE EQUILIBRIUM GAS FILM THICKNESS VIA A SUBPROGRAM AFILM# 
CAuL FILM 

C PRINT FILM THICKNESS AND STEADY SIATE ROTOR FORCES 
IF(ELIM.GT.O.O) GO TO 110 
ELIM-.0002 

110 WRITE(NW»1251) ELIM 
WRITE(Nw»1500) 

WRITE (Nw» 1020) 

C CONVERT HINIT TO MILS FOR PRINTOUT PURPOSES 
HINIT=KlNIT*10u0. 

WRITE (Nw* 1250)HINIT *FA*GL 
C RECONVERT TO INCHES FOR CALCULATION PURPOSES 
HINIT=HINIT/1000. 

WRITE(NWd020) 

IF(MODEL.EQ.1)wRI1E(NW» 1?20)M1»M2*M3»K1*C1*K. 

IF(M00EL.EQ.2) wRIT£(NW» 1340)M1»M2*M3*K1*K2 
IF(M0DEL.EQ.3) wRITE(Nw»1330)M3 
C CHANGE UNITS 

M1=M1 /386.069 
M2=M2/ 386.069 
M3=M3/ 386.069 
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WRITE (NW» 1020) 

WRITE(Nw»1140) 

IF ( IFOFT .EQ.8) WK I T E ( NW » 1 0^0) 
IF ( IFOFT.Ey.8) WRITE (NW» 1 150) 
IFdFOFT.EQ.S) WR I TE ( NW , 1 040 ) 
IF ( IFOF T,EQ,8) write (NW, 1400) 
IFdFOFT.EU.B) WK I T E ( NW d 1 60 ) 
IFdFOF T.EQ.8) wRI TE (Nw» 1400) 
IF(IF0Ft.EU. 8) WKlTt (N w, 11 70) 
IF(lFOFT.EGi.d) wRItE(Nw»1400) 
IF ( IFOFT ,EU,8) WR I TE ( imw , 1 1 80 ) 
IF (IFOFT,E(3.9) WR I TE (NW . 1400 ) 
IF(IF0FT.EQ.9) WK I TE ( NW d 1 50 ) 
IF (IFOFT. EQ.^) WRI TE (Nwd400) 
IF ( IFOFT .EQ. 9) WR 1 1 E ( NW d 040 ) 
IF(IF0FT,E0.9) WRlTE(Nwdl60) 
IFdFOFT.EU.9) WR I TE (NW ♦ 1 04 0 ) 
IF(IF0FT.EQ,9) write (NW» 1400) 
IF(IFoFT,Eu.9) WR I T E ( NW d 1 70 ) 
IF (IFOFT. EO, 9) WR I T E ( n w » 1 400 ) 
IF(IF0FT.EQ.9) WRlTE(NwdiaO) 
IFdFOFT.EO.lO) WHITE (NWd400) 
IFdFOFT.EO.lO) WHITE (NWdlSO) 
IF(IFOFT.EQ.IO) WH I TE ( NW d 400 ) 
IFdFOFT.EO.lO) wRlTE(NWdl60) 
IF ( IFOF T .EO. 10) WRITE (NWd400) 
IF (IFOFT ,EO. 10) write (NWd040) 
IFdFOFT.EU.lOIwRHF (NWtli /0) 
IF (IFOFT.ECJ. 10) WHITE (NW* 1040) 
IFdFOFT.EO.lO) WH1TE(NW*1400) 
IFdFOFT.EO.lO) WHITE (NW* 1180) 
IFdFOFT.EO.lO) WH1TE(NW*1400) 
IFdFOFT.EQ.ll) wH 1TE(NW* 1400) 
IFdFoFT.EQ.ll) WR1TE(NW»1150) 
IFdFOFT.EQ.ll) WHITE (NW* 1400) 
IF(IFOFT.EO.ll) WHITE (NW* 1 160) 
IF(IFOFT.EQ.ll) WRITE (NW* 1400) 
I F (I FOF T . F. 0 . 1 1 ) WR I TE ( NW ♦ 1 1 7 0 ) 
IF(IFOFT.EO.ll) WRITE (NW* 1400) 
IF(IFOF7.EQ.ll)WRITE(NW*10‘+0) 
IF (IFOFT.EO.il) WRITE (Nw*lldO) 
IF (IFOFT. EO. 11 ) WRITE (Nw* 1040) 
WRITE(NW»1531) 

WRITE (Nwd020) 

WRITE (NWdl90) 

IF(ICC.E0.12) wklTE(Nw*1040) 

IF ( ICC.EQ. 12) WRI IE (NW* 1200) 
IF(ICC.E0.12) WRI (E (NW* 1040) 
IF(ICC.EQ.12) WHITE (NW*1400) 

IF (ICC.E(J. 12) WHITE (NW* 1210) 
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IF(ICC.EQ.13) WRIT£(NW»1.200) 

IF(ICC.EQ.13) i«lRn£(NW»l400) 

IF<ICC.£Q.13) MRITE(NM«1040) 

IF(ICC.EQ.13) MRITE<NW*1210) 

IF(ICC.EQ.13) »<RITE(NW»1040) 

MRIT£(NW»1400) 

WRITE(NW»1020) 

C PRINT PARAMETERS ASSOCIATEO WITH FORCING FUNCTION 
IFM7=IFOFT-7 

GO TO (120«130tl40f ISO) «IFH7 
C STEADY STATE SINE FORCE 
120 REAO(NR»90) AMPLOfFREU 
WRlTE(NWtl320) FREQ 

IF (MOTION.EO.O) WRITE (NW»1S20) AMPLD 
IF (MOTION. EU.l) WRITE (NW»lb30) AMPLD 
GO TO 160 

130 REaO(NR*90) peaks»pduks 
WRITE <NW» 1310) POURS 
IF (MOTION, EQ.O) WR I TE ( NW ♦ 1S20 ) PEAKS 
IF(MOTION.EQ.l) WRITE(NW,1530) PEAKS 
GO TO 160 

140 READ (NR*90) PEAKH*PDURM 
WRITE (Nw* 1300) PDURh 

IF (MOTION.EO.O) WR I TE ( NW ♦ 1 S20 ) PEAKH 
IF(MOTION.EU.l) wRITE(NW»1530) PEAKH 
60 TO 160 

150 WRITE (NW» 1260) 

IF (MOTION.EO.O) WR I TE ( NW ♦ 1 26 1 ) 

IF(MOTION,EQ,l) WR I It ( NW » 1262 ) 

C ABOVE IS HEADING FOR SHOCK PULSE TABLE 

C NOW READ THE NUMBER OF INCREMENTS IN THE PULSE AND THE INCREMENTAL 
C TIME INTO WHICH PULSE IS DIVIDED 
REA0(NR.154) MINCPtTINCP 
154 FORMAT ( I5*F10. 4) 

C INITIALIZE FIRST Yd VALUE 
INCPO=NINCP+l 

READ(NR,90) (PUVAL ( I ) ♦ 1=2. INCPO) 

151 PUVAL(l) =0.0 
TP(1)=0,0 

153 DO 159 I=2.INCPn 
EYE=I-1 

TP(I)=EYE*TINCP 

IM1=I-1 

IPlsI-^1 

159 IFI.NOT, (I.EQ.I.OR.I.EU.INCPO) ) PULDOT(I)s(PUVAL(IPD- 
IPUVAL(lMl) )/(2.*TTNCP) 

DO 158 I=l,INCPO 
MAPGE=(I/45)*45-I 
IF(MAPGE.EQ.O) WRIT£(NW»1531) 

158 WRITE (Nw» 1265) TP ( I ) ♦ PUVAL ( I ) 
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WRlTE(Nw»1500) 

ANINCP=NINCF 

pdurt=animcp*tincp 

WRITE (NW*1510) PDURT 
1265 F0RMAT(13X»2F12.7) 

PULDOT ( 1 ) = (PUVAL ( 2) -PUVAL ( 1 ) ) /T INCP 
PULDOT ( INCPO) = (PUVAL ( INCPO) -PUVAL (NINCP) ) /TINCP 
C FORM A NEW TABLE FOR ThE INTEGRAL OF YB 
IF (MOTION. EQ. 1 ) CALL FORM! 

WRITE(NW*15J1) 

C READ PROGRAM CONTROL CAROS 

160 READ(NR»91) I PSTEP ♦ M , MAX I NC ♦ NCYCLS ♦ TPRINT » TMAX ♦ H» E » HM ♦ BOGUS 
91 FORMAT (4I5»6E10. A) 

WRITE (NW»1500J 
WRITE(NW*1020) 

NCONE=NCYCLS 

IF (NCYCLS, EQ.O) NCON£= NCYCLS+1 

WRITE (NW» 1280) TMAX »M AX I NC »NCONE »h » TPRINT » IPSTEP » E »HM tM » BOGUS 
C READ AND WRITE INITIAL CONDITIONS 

READ (NR ♦ 90) TINIT»Y1IN»Y2IN»Y3IN»V1IN»V2IN»V3IN 
WRITE (NW» 1270) T I N I T » Y 1 1 N » V 1 1 N ♦ Y2I N* V2 IN ♦ Y3INt V3IN 
1010 FORMAT ( I5»75H 

1 ) 

1020 FORMATdH dlRdh*)) 

1030 FORMAT (40X» 15HPROGRAM OPT IONS/20X ♦ 58H0PT IONS SELECTED FOR THIS PAR 

ITICULAR RUN ARE DESIGNATED BY/70X»8H /71 X »6H0PT I ON/70X ♦ 

28H ) 

1031 FORMAT (iOX»25HCH01CE OF VIBRATION MODEL/) 

1040 FORMAT (20X.50H ) 

1050 FORMAT (20X»50H1 = THREE-DEGREE OF FREEDOM SYSTEM USING ISOLATORS) 
1060 FORMAT (20Xf53H2 = THREE-DEGREE OF FREEDOM SYSTEM USING RIGID MOUNT 
IS) 

1070 FORMAT (20X,52H3 = SIMPLE MODEL* ENTIRELY RIGID EXCEPT FOR GAS FILM 

1 ). 

1080 FORMAT (/10X,34HGAS FILM FUNCTIONAL REPRESENTATION/15X* 19HL0AD CHAR 
lACTERISTIC/) 

1090 FORMAT (20X»35H4 = EXPONENTIAL LOAD CHARACTERISTIC) 

1100 FORMAT (20X.38H5 = POWER FUNCTION LOAD CHARACTERISTIC) 

1110 FORMATdOX,22HDAMPING CHARACTERISTIC/) 

1120 FORMAT(20X,35H6 = EXPONENTIAL DAMPING COEFFICIENT) 

1130 FORMAT(20X,38H7 = POWER FUNCTION DAMPING COEFFICIENT) 

1140 FORMAT (lH0*9X»27hTYPE OF EXCITATION FUNCTION/) 

1150 FORMAT (20X, 27H8 = STEADY SINUSOIDAL INPUT) 

1160 FORMAT (20X* 19h9 = HALF SINE SHOCN) 

1170 FORMAT (20X* 44H10 = HAVERSINE SHOCK 

1180 FORMAT(20X» 44H11 = ARBITRARY SHOCK (TABLE) 

1190 FORMAT dhO* 9X*26HIS CALCOMP OUTPUT DESIRED-) 

1200 FORMAT(20X» 8H12 = YES) 

1210 FORMAT(20X* 7H13 = NO) 

1220 format dHO»20X»53MPARAMETERS OF THE THREE-DEGREE SYSTEM USING ISOL 
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1AT0RS//12X, 7HMASS OF* 5X» 14HEFFECTIVE MASS* 4X»7HMASS OF* 

2 6X*9HSTIFFNESS* 6X * 7HDAMPING* 7X* 7HFLEXURE/13X»6HCASING*9X 

3 * 8H0F plate* 7X*6HR0T0R * 10X*2HOF*12X*2HOF* 9X* 9HSTIFFNESS/ 
426X*11HAND FLEXURE*iax* 9HIS0LAT0RS* 5X* 9HISOLATORS/ 

514X*4H(LB) *10X*4H(LB) *10X*4H(lB) * 9X* 7H(LB/IN>* 5X* 1 1H (UB-SEC/IN) 
6* 4X* 7H(LB/IN)//10X*5(G10.3*4X) *E12.4) 

1230 FORMAT (IhO* 9X*23hOESCRlPTlON OF GAS FILM// 

120X*50HTHE EXPONENTIAL LOAD CHARACTERISTIC IS OF THE FORM/ 
225X»22HL0A0 (LB) = A*EXP ( B*H ) /20X * 5HWHERE/ 

325X*27HH = FILM THICKNESS ( INCHES) /25X, 4HA = *E12.4*5H (LB)/ 

425X* 4H8 = *E12.4* 7H( 1/IN)) 

1240 FORMAT (1H0*9X*53HTHE POWER FUNCTION DAMPING COEFFICIENT IS OF THE 
1FORM/20X*21HB(LB*SEC/IN) = C*H**D/10X*5HWHERE/ 

220X*28HH = FILM ThICKNESSS (INCHES)/ 

320X* 4HC= *E12.4*10H LB-SEC/IN / 

420X* 4HD= *E12.4* 7H (1/In)) 

1250 F0RMAT(10X*29H£QUILIBRIUM FILM THICKnESS IS*612.5* 5H MILS// 
115X*56HASSUMING STEADY-STATE AERODYNAMIC FORCES ON THE ROTOR OF* 
2G12.5* 3H L0/ 15X*34HAND A LOCAL GRAVITATIONAL FIELD OF*G12.b* 

315H times STANDARD) 

1251 FORMAT (/10X*85HNON-LIN£AR LOAD AND DAMPING FUNCTIONS ARE HELD CONS 
ITANT AT FILM THICKNESSES LESS ThAN*G10.3*8H INCHES.) 

1260 FORMAT ( 10X*62HINPUT MOTION IS A SHOCK PULSE DESCRIBED BY THE FOLLO 
IWING TABLE) 

1261 FORMAT (/20X*4HTIME»2X.12HDISPLACEMENT/19X*5H(SEC) »4X* 

1 8H( INCHES) ) 

1262 FORMAT (/20X*4HTIME*2X* 12HACCELERAT ION/ 19X*5h( SEC) *4X*5H(G^S) ) 

1270 FORMAT(1hO*9X*29hINITIAL OR RESTART CONDITIONS/ 

120X* 12HIMTIAL T IME * 1 1 X * 1H=* G1 0 .3* 8H SECONDS/ 

220X*24HDISPLACEMENT OF MASS 1 =*G10.3* 7H INCHES/ 

320X* 18HVELOCITY OF MASS 1 *SX* 1H=*G10,3* 7H IN/SEC/ 
420X*24HDISPLACEMENT OF MASS 2 =* G10.3* 7H INCHES/ 

520X* 18HVEL0CITY OF MASS 2 * 5X . 1 H= * G1 0 . 3 * 7H IN/SEC/ 
620X*24HDISPLACEMENT OF MASS 3 =*G10.3* 7H INCHES/ 

720X* 18HVEL0CITY OF MASS 3 * 5X * 1H= * G1 0 • 3 * 7H IN/SEC/) 

1280 FORMAT (/ 10X*26HPROGRAM CONTROL PARAMETERS//15X * 19HDESIRE0 TERMINA 
1TION/20X*26HTIME VALUE WILL NOT EXCEED * F 1 0 .4 * 8H SECONDS/20X* 

248HAS A SAFETY MEASURE WE WILL NOT PERMIT MORE THAN,Ib,6H STEPS/ 
320X*25HCALCULATION WILL CONTINUE* I5*2BH CYCLES AFTER PULSE DURATIO 
4N/15X*9HSTEP SIZE/20X*45HINITIAL STEP SIZE FOR CALCULATION PURPOSE 
5S IS*G10.3*8H SECONDS/20X»28HVALUES ARE PRINTED OUT EVERY ,G10.3* 
68H SECONDS/20X*5HEVERY* IS*38nTH POINT IS PLOTTED IF CALCOMP IS USE 
7D/20X*41HERROR CRITERION FOR INTEGRATION TECHNIQUE*G10.3/ 20X*20HM 
8INIMUM STEP-SIZE =*G12.5*8H SECONDS/20X*48HNO. OF INTEGRATIONS BEF 
90RE STEP-SIZE IS DOUBLED** IS/20X*48HNORMALIZED ERROR LIMIT FOR INT 
lEGRATION TECHNIQUE*F10,4) 

1300 FORMAT (10X*37HINPUT MOTION IS HAVERSINE SHOCK PULsE/ 

215X* 17HPULSE DURATION IS*F10.4* 8H SECONDS/ ) 

1310 FORMAT ( 10X*37HINPUT MOTION IS HALF-SINE SHOCK PULSE/ 

215X* 17HPULSE DURATION IS*F10.4*8H SECONDS/) 
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1320 FORMAT ( 10X,59HU£SCKIPTI0N OF INPUT MOTION IS STEADY-STATE SINE EXC 
1 1 T AT I ON/ 

3l5X»12HFhEQUENCY IS»F10.4, 3h HZ/ ) 

1330 FORMAT ( lOXfSSHTHRuST HEARING uN RIGID FLEXURE (SINGLE DEGREE OF PR 
lEEDOM) /10X» I3hMASS OF ROTOR/ 1 3X (L6) /I OX ♦GIO .3/) 

13A0 FORMAT(10X»55HPARAMErERS OF T riE SYSTEM ON RIGID MOUNTS ARE AS POLL 
10WS//13X»bOHEFFECT IvE MASS EFFECTIVE MASS MASS OF ROTOR, 

26X,9HST IFFNESS,9X,«HFLEXURE /12X»16H0F LOWER HOUS I NG, 2X , 1 9H0F STAT 
30R + FLEXURE, IbX, 16hOF LOWER hUUS ING ,6X , 9HST IFFNESS/ 

4l6X’l?H(MASSl ) , (LB) »4X,1^H(MASS2) , (LB) »7X,12h(MASS3) , (LB) ,6^» 

b 7H (LB/ IN ) , 1 1 X , 7h (LB/ I N ) / 1 6X , G1 0 » 3 , 6X , G1 0 • 3, 9X , 61 0 . 3, 2 ( 5X , F 1 ? 4» ) 
13b0 FORMAT(30X,lAHPRO6RAM OUTPUT) ♦ i u . ^X , 0 . J, 2 ( bX , E 12. A) ) 

1360 FORMAT ( l4X,4HTlMt, 6X, BHGAS FILM, 4X, BriPEOESTAL, ?X,2MY1,10X, 
12HY2,lUX,2HY3,10X,2HVj,10X,2hV2,10X,2HV3/ 


224X, RHThlCxNESS, 2X , 1 UHDEFLECT 1 VE/2bX , 7H(Y3-Y2), bX, 7M(Y2-Y1)/ 
3 12X. 9H(SEC0N0S), 4X,6h(MlLS), 6X»6H(MILS), 6X,6H(MILS), 6X, 

46H(MILS), 6X,6h(MlLS), 4X , 1 OH (MILS/SEC) , 2X , 1 OH ( MILS/SEC ) , 2X , 
510H(MILS/SEC) ) 

1370 FORMAT(10X,9(F10.4,2X) ) 

1380 format ( 2bX , bOHREST AR r CONUlTlONui (FINAL) AND IDENt IF ICA f I ON CODE/ 
114X,4HTIME, 9X,2hYl,10X,2hy2,loX,2HY3,10X,2MVl,10X,2HV2il0X,2HV3, 

2 6X,bHNCASE/ 14X,bH(SEC), 6X,6h(MILS), 6X,6H(MILS), 6X,6H(MILS), 

3 4X, lOH (MILS/SEC) , 2x . 1 OH ( M I LS/SEC ) , 2X , 1 OH ( M I US/SEC ) / 
411X,7(F10.4,2X) ,Ib) 


1235 FORMAT (1P0,9X,27hCHARACTERIST1CS OF GAS FILM// 

120X,b3HThE POWER FUNCIION LOAD CHARAC TER I ST IC IS OF THE FORM// 
225X»18HL0A0 (Ld) = /2C)A , bhWHERE/ 

325X»27HH = FILM THICKNESS ( INCHES ) /2bX , 4HC = ,£12.4, 

418H LB-SEC/IN.** ( 1 +D) / 


525X, 4HD = ,E12.4 ) 

1245 FORMAT ( IHO , 9X,53h 1 HE EXPONENTIAL DAMPING COEFFICIENT IS OF THE FOR 
IM /20X»2bHB (LB«SEC/IN) = C*EXP ( D«h ) / 1 OX ,bHwHE RE/ 

220X,28hh = FILM ThICKNESSS (INCHES)/ 

320X, 4HC = ,E12.4,15H LB-SEC/ ( IN**2) / 

420X, 4H0 = ,E12.4»21H (DIMENSIONLESS UNITS)) 

1400 F0RMAT(/) 

IbOO FORMAT(//) 

1510 F0RMAT(// 10X,16HRULS£ DURATION =, G12.b,4H SEC) 

1520 F0RMAT(15X,17hPEAK AMPLITUDE IS,F10.4,20H INCHES DISPLACEMENT) 

1530 F0RMAT(15X,17HPEAK AMPLITUDE IS, FI 0 .4, 1 7H G»SS ACCELERATION) 

1531 FORMAT(lhl) 

999 RETURN 

END 

SUBROUTINE S IM INI ( S I , FX , A,E , REL,MAX I T ,FK»MRET ) 

EXTERNAL FX 
LOGICAL REL 
DIMENSION A(2) 

AA=A(1) 

BB=A (2) 

H = (BB-AA) / 6.0 
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XMID = (AA+BB) * .5 
Q = 2. * FX(XMID*FK) 

SI = FX(AA»FK) ♦ Q ♦ FX(BB»FK) 
PR s H*(S1*Q) 

IF(MAXIT.EQ.l) GO TO 150 
DO 100 1T=2*MAXIT 
ISTOP = 2** IT - 1 
H= .5*h 
SsO. 

DO 80 M=l*IST0Pf2 
X=AA*3.0*M*H 
80 S=S*FX(X»FK) 

SI=H*(4.*S ♦SI) 

DIFF=ABS(SI-PR) 

IF(rEL) 0IFF=DIFF/ABS(SI) 
IF<DIFF.LE.E) GO TO 160 
PR=S1 
Sl=Sl+S+S 
100 CONTINUE 
RETURN 
150 SI=PR 
160 MRET=0 
RETURN 
END 


SUBROUTINE SOLVE 
REAL M1»M2»M3»C1»K1*K2 

COMMON TPRINT»TMAX*IPSTEP»NCYCLS*M0TI0N»TITLE(15) ♦Y8A*YBINT (200) 
COMMON NR ♦ NW ♦ NCASE ♦ MODEL ♦LOADF»IDAMP»IFOFT» ICC* 

1 Ml «M2*M3*K1 *C1«K2«GL«FA«AL0A0«A0AMP»8L0A0»B0AMP* 

2CL0AD»CDAMP*0L0AU*00AMP»HINIT*1FM7*AMPLD*FREQ*PEAKS*PDURS*PEAKH* 
3PDURH*NINCP*TINCP»TP(200) *PUVAL(200) *PULDOT(200) *TINIT»Y1 IN* Y2IN* 
AY3IN*V1 IN* V2IN*V3IN*T1ME (6000) *Y1*Y2*Y3*V1*V2*V3*N*VEC<6) *0VEC(6) * 
5MAXINC*H»HM*M*E*AK0*NE*T* ELIM*INOX* YB*YBDOT* 

60MEGA*WM0*WY3Y2H*BY3Y2h 

COMMON FILMT (6000) *Y2MY1 (6000) * VAPEAK (6000 ) * YBPLOT (6000) 

COMMON TZ*ITIME*YBPREV 

C THIS SUBROUTINE SOLVES THE EQUATIONS OF MOTION*WHICH WILL BE DEFINED 

C IN SUBROUTINE OIFFEQ 

C 

C INITIALIZE THE VECTORS 

C INITIALIZE STORAGE FOR MAX AND MIN VALUES 
ITIME=1 
IVALsl 
IMAXIM=0 
IFM7* IFOFT-7 
ISTEAD=1 
WRITE(NW*1111) 

IF(M0DEL.EQ.3) GO TO 10 
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WPITE(NW,1112) 

N = 6 

'/ECU ) = vl IN 
Vr.C (2) =V2IN 
VEC ( J) =V3IN 
VEC (4) =Y 1 IN 
VEC(5)=Y2IN 
VEC (6 ) =y3IN 
GO TO 20 

iO WRITE (Nw» 1 1 13) 

N = 2 

VEC(1)=V3IN 
VEC (2) =Y31N 
20 T=TINIT 

^ integrations before step SIZE IS DOUBLED 

C INITIALIZE STORAGE FOR T- PRINT VALUES 
TPRVAL=r*TPRINf 
KELP=i 

00 200 kOunt=1 ,MAaINC 

C SET FLAG OFF FOR PRINT 
IPRINT=2 

C INITIALIZE M AND NE 
M=MOLD 

IF(K0UNT.£Q, 1) NE=0 
IFIKOUNT.NE.I ) NE=] 

TADINC=T*2.«H 

IF ( TADINC .LT. TPRVAL ) GO TO 1200 
TOINC=T*h 

IF ( T0INC.lt, TPRVAL ) GO TO 1190 
C T* OT WILL EaCEED PkINT VALUE 
C SAVE OLD DT VALUE. I.E. SAVE H 
HOl.n = H 
H=TRWVAL -I 
C SAVE NEW 01 VALUE 

hnew=h 

C TURN print flag ON 
IPR1NT=1 

C SET M=10»00U » SO THAT Ul WILL NOT BE DOUBLED 

1190 M=10000 
C RONKU INTEGRATION 
1200 TZ=T 

CALL RONKU (T»VEC»OVEC»h»HM»N.NE»M*E»AKO) 

YBPREV=Yd 

TEST PRINT flag on. 1=0N. 2 = 0FF 

I F' { IPR InT . EQ .2 ) Gu To 1536 
C TEST DT =nEw OT 

IR(H-PN^;w) 1536. 10. 1536 

C RESET 0 T TO OLD D T 
1210 H=hOLD 
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C INCREASE FOR NEXT PRINT VALUE 
TPRVAL=T^TPRINT 

C FORM AN INDEX TO TAG NO. OF ITERATIONS 
INDX=KELP 

IF(MOTION.EQ.O) YBPLOT (KELP) =YB 
IF(MOTION.EQ.l) YBPLOT (KELP) *YBA 
TImE(KELP)*T 
MAPGEs (KELP/45) *4S-K£LP 
IF(MAPGE.EO.O) WRITE(NW*1U8) 

IF(M0DEL.EQ.3) GO TO 1490 
Y2MY1 (KELP)=VEC(4)-VEC(5) 

FILMT(KELP)=HINIT+VEC(6)-VEC(5) 

WRITE (NW. 2222) T.FILMT (KELP) * (VEC ( 10) . 10=4.6) . ( VEC ( 10) . 10=1 .3) 
1Y2MY1 (KELP) .YB 
GO TO 1530 

1490 FILMT (KELP)=MINIT+VEC(2)-YB 

WRITE (Nw. 2222) T .FILMT (KELP) .VEC(2) .VEC( 1 ) . YB.YBA 

1530 IF(FILMT (KELP) .LT. 0.0) GO TO 1531 
GO TO (400.500.600.700) .IFM7 

1531 IFdFoFT.EQ.S) GO TO 400 
WRITE(NW.1115) 

1115 F0RMAT(/5X.47HRUN COMPLETE AT THIS POINT. ROTOR HAS BOTTOMED.) 
GO TO 666 

400 IF (KELP. LT, 3) GO TO 1535 
D1=FILMT (KELP) 

D2=FILMT(KELP-1) 

D3=FILMT (KELP-2) 

IF(D1.LT.D2.AN0.D2.GT.D3) GO TO 410 
IF(D1.GT.D2.AND.02.LT.D3) GO TO 420 
GO TO 1535 
410 IMAXIM = 1 
GO TO 425 
420 IMAXIM=2 
425 VAPEAK(IVAL)=02 
IVAL=IVAL^1 
C TEST FOR MIN OH MAX 

GO TO (1535.430) .IMAXIM 
C TEST FOR STEADY-STATE 

430 IF(ISTEAO.NE.l) GO TO 661 
IF(IVAL.LE.5) GO TO 1535 

ERRC= (VAPEAK ( I VAL-J) -VAPE^K ( IVAL-5) ) /VAPEAK ( I VAL-3) 
IF(ABS(ERRC) .GT,.01)‘ GO TO 1535 
ISTEAD= 2 
GO TO 1535 
C IFOFT = 9 

500 IF (T.LT. POURS) GO TO 1535 
510 IF(KELP.LT.3) GO TO 1535 
D1=FILMT(KELP) 

D2=FILMT (KELP-1) 

D3=FILMT (KELP-2) 


IF (Dl.LT.n?.AN0.D<£.GT.D3) GO TO 515 
IF (D1 .GT.Oa. AN0.D2 .lt. 03) GO TO 515 
GO TO 1535 
515 IMAXIM=IMAXIM+1 

IF (NCYCLS.EQ.O ) GO TO 520 

IF ( IMAXIM.GT. (2-»NCYCLS) ) 60 TO 666 

GO TO 1535 


GO TO 1535 


700 


1535 

1536 
200 


520 IF(IMAXIM.GE.2) 60 TO 666 
60 TO 1535 
IFOFT = 10 
600 IF (T.LT.PDU«H) 

GO TO 510 
ANINCP=NINCP 
IF (T.LT.ANINCP*TINCP) 

GO TO 510 
KELP=KELP+1 

IF(T.GT.TMAX) GO TO 660 
CONTINUE 
WRl TE (NW» H 14) mAXInC 
Go To 6fa6 

C MAKE A HEADING AND 
1111 FORMAT ( Ihl » 19X»43h 

1 20X» 43H 0 0 

2 20X. 43H 0 0 

3 20X» 43H 0 0 

4 20X» 43H 000 

5) 

1112 FORMAT(120H 

lEMENT VELOCITY 
2 120h 

30R OF MASS 
4 120H 

5S) (IN/SEC) 


GO TO 1535 


PRINT RESULTS 

000 U U TTTTT PPPP U U TTTTT / 
UU T PPUU T / 

U U T PPPP U U T / 

UUTP UUT/ 

UUU T P UUU T //// 

TIME FILM DISPLACEMENT DISPLACEMENT DISPLAC 

VELOCITY VELOCITY RELATIVE FORCING / 
THICKNESS OF MASS 1 OF MASS 2 OF ROT 

1 OF MASS 2 OF ROTOR UISP. BET. FUNCTION / 
(SEC) (INCHES) (INCHES) (INCHES) (INCHE 


1113 FORMAT(95HO 


time 


(IN/SEC) 
F ILM 


(IN/SEC) MASS2 ♦MASSl 
DISPLACEMENT VELOCITY 


FORCING 

75H 

function 

75H 

(G*5) 


OF ROTOR 


OF ROTOR 


(SEC) 


(INCHES) 


(IN/SEC) 


THICKNESS 
/ 

(INCHES) 

/) 

1114 FORMAT( /5X.38HMAXIMUM NO. OF TIME STEPS HAS EXCEEDED*I5) 
2222 FORMAT (1X.K7.4»E11.4.3E13.4»5E12.4) 

1116 FORMAT (/5X.42HSTEADY-STATE CRITERION HAS BEEN SATISFIED.) 

1117 FORMAT (/5X,33HSPECIFIED TIME HAS BEEN ATTAINED.) 

1118 FORMAT(lHl) 

660 WRITE(NW»1117) 

GO TO 666 

661 WRITE (NW. 1 1 16) 

666 IF ( ICC.EC^. 12) CALL CALCOM 
RETURN 
END 


(INCHES) ) 
FORCING 
/ 

FUNCTION 


(INCHES) 
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SUBROUTINE D IFEQ ( N ♦ T ♦ VEC ♦ DVEC ) 

REAL M1»M2»M3»C1*K1,K^ 

DIMENSION VEC(6) *0V£C(6) 

COMMON TPRINT*TMAX*IPST£P»NCYCLS» MOT ION* TITLE (15) »YBA»YBINT (200) 
DIMENSION A(2) 

COMMON NR»NW»NCASE»MOOEL»LOADF» IDAMP*IFOFT«ICC* 

1 M1*M2«M3«K1 »C1*K2«6L«FA«AL0A0*ADAMP«BL0AD*BDAMP» 

2CL0AD*C0AMP,DL0AD»D0AMP,H1NIT.1FM7*AMPLD*FREQ»PEAKS*PDURS*PEAKH 
3PDURH,NINCP*TINCP*TP (200) »PUVAL (200) »PULDOT (200) ♦ T INIT » Y 1 IN* Y2I 
4Y3IN»V1IN*V2IN*V3IN,TIME(6000) *Y1*Y2»Y3*V1*V2*V3* 

5NDUM*DDUM (6) »DVUM (6) * 

5MAXINC»H»HM*M»E*AKO»N£*TTDUM» £LIM»INDX* Y8»YBD0T» 

60M£GA*WH0»wY3Y2H*BY3Y2H 

COMMON FILMT (6000) » Y2MY1 (6000) »VAPEAK (6000) tYBPLOT (6000) 

COMMON T2* ITIME,YBPR£V 
external FX 

c this subroutine defines the differential equations sOLvED BY RONKu 

DATA MAXIT»ERR/100».01/ 

WEX (X)=AL0AD*EXP(8L0AD*X) 

WPF (X)=CL0AD«X«*DL0AU 
BEX (X) =ADAMP*EXP (BDAMP*X) 

BPF (X)=CDAMP*X*<»DOAMP 
PI=3. 1415926535 
IFM7=IFOFT-7 
IF (MODEL. EQ. 3) GO TO 20 
DVEC(4)=VEC(1) 

DVEC(5)=VEC(2) 

DVEC(6)=VEC(3) 

20 GO TO (30»60»90»120) tlFM? 

30 0MEGA=2.*PI«FRE0 

YB=AMPL0*SIN(0MEGA*T) 

YBDOTsOMEGA*AMPLD*COS(OM£GA*T) 

YBA=-YB*(2.*PIi»FREQ)**»2/386.069 
IF (MOTION. EQ.O) GO TO 130 
YB=-386.069/OMEGA**2*AMPLO*SIN(OMEGA*T) 
YBDOTs-386.069/OM£GA*AMPLD<*COS(OMEGA*T) 

GO TO 130 
60 YBsO. 

IF(T.LT.PDURS) YB=PEAKS*SlN (PI*T/PDU«S) 

YBOOTaO. 

YBA=YB 

IF(MOTION.EQ.l) GO TO 65 

IF (T.LT. POURS) Y6D0T"PI/PDuRS*PEAKS*C0S(P1*T/PDURS) 

GO TO 130 

65 IFd.LT. POURS) Y8*386. 069*PEAKb*PDU«S/P I* ( T-PDURS/PI<» 

1 SIN(PI»T/POURS) ) 

IF (T.LT. POURS) YBDOTa386.069*P£AKS*POUHS/PI*( l.-COS(PI*T/POURS) ) 
IF(.NOT. (T.LT. POURS) ) YBDOT*7 72. 138*PEAKS*P0URS/PI 

IF(.NOT. (T.LT. POURS) ) YB*386. Ob9*PEAKS*POURS**2/PI ♦ YBOOT* ( T-PDURS) 
GO TO 130 
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90 


YB = 0. 

IF (T.LT.PDUKH) YB=0.S«PEAKH»( l.-C0S(2.*PI*T/PDUPH) ) 

Y8DOT=0. 

Y0A=YB 

IF (MOTION. EQ. 1 ) GO TO 95 

IF (T.LT.PDUPH) YBUOT=PI/POURh»SIN (2.*PI<^T/PDURH)*PEAKH 
GO TO 130 

95 IF (T.LT.PDUrJH) Yb = 386 . 0 69/4 . «-PE AKH* ( T»»2-P0UHIH»»2/ ( 2. <^PI»-»2) » 

1 ( 1 .-CoS (2.*PI*T/PDUPH) ) ) 

IF (T.LT.PDUrH) YBuOT = 386.069/4.«-PEAKM* (2.*T-PDUrH/PI« 

1 SIN (2.*PI<»T/P0UPH) ) 

IF ( .NOT. (T.LT.PDUkh) ) YBU0T = 193. 035«PEAKH*PDURH 

IF ( .NOT . (T.LT.PDUKH) ) V 8=386 . 069/4 . *P£AKH*PDURH**2* YBDOT* ( T-PDURH) 
GO TO 130 
120 ANINCP=NINCP 

pqupt=anincp*tincp 

IF (T.GT.PDURT) YB = 0.0 
IF (T.GT.PDUPT) YBDO r = 0.0 
IF (T.GT.PDURT. AND. MOTION. EU.O) GO TO 130 
IF(T.GT.PDURT.AND.MOriON.EQ.l) GO TO 125 
NP1=NINCP+ 1 

CALL TLU(T»YB.TP.PUVAL»NP1) 

IF (MOTION. £0.1 ) GO TO 125 
CALL TLU(T.YBDoT»TP.PULDOT»NP1) 

GO TO 130 
125 YBA-fB 
1240 A(1)=TZ 
A (2>=T 
IOP=l 

IF (T.GT .PDURT) GO TO 1270 
IF(ITIME.EQ.l) GO TO 1260 
1250 MRET=12b 

CALL SIM INI (YB»FX,A.EWK»0»MAaIT»1.»MRET) 

IF(MRET.EQ. 126) GO TO 126 
Y0 = 386. 069<*YB + Y0PKEV 
GO TO (127»12H0.1265) tlOP 

1260 ITIME=2 
IGR=1 

YBPREY=0,0 

IF(TZ) 1261 » 1250» 1261 

1261 A(1)=0.0 
A(2)=TZ 
I0P = 3 

GO TO 1250 
1265 YBPRE\/ = YB 
GO TO 1240 
C SAVE YB VALUE 
C T GREATER THAN PDURT 
1270 IF(IGR.EQ.3) GO TO 1275 
A (2) =POURT 
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IGR=2 

I0P=2 

GO TO 1250 
1275 Y8=YBPDUR 
GO TO 2126 
1280 YBPUUR=YB 
GO TO 127 

126 WRITE(NW*102) MAXIT.ERR 
MRITE(NW*103) 

127 IF(T.LT.PDURT.0R.IGR.LE.2) GO TO 2127 

2126 YBDOT=OOTDUR 
GO TO 2128 

2127 CALL TLU(T«YB00T»TP»Y8INTfNPl) 

YBOOT=YBDOr*386.069 
IF{IGR.NE.2) GO TO 130 
IGR=3 

ootour=ybdot 

2120 YB=YB+Ya00T*(T-PDOKT) 

130 IF(M00EL.NE.3) GO TO 150 
ARGU=VEC (2)-YB*HINIT 
IF (ARGO.lt. .0002) ARGO=ELIM 
IF (LOAUF .EQ.4) WY3Y2h=WEX ( ARGO) 

IF (LOADF.EQ.4) wmO=wEx (HIM T ) 

IF (LOAOF .EQ.5) wY3Y2H=WPF (ARGO) 

IF (LOAOF .E0.5) WhO=WPF (HINIT ) 

IF(IDAMP.E0.6) BY3Y2H=BEX(ARG0) 

IF ( IDAMP.EQ.7) I:JY3Y2H=BPF ( AHGO) 

GO TO 190 

C SELECT FONCTIONS 

150 ARGU=VEC(6)-VEC(5)*HINIT 

IF (ARGO. LT. .0002) ARGU=ELIM 
IF (LOAOF. EQ. 4) WY3Y2h=WEX ( ARGO) 

IF (LOAOF. EU. 4) wHO=WEX(HINIT) 

IF (L0A0F.E(i.5) W Y 3Y2 h= WPF ( ARGO ) 

IF(LOAOF.EQ.S) WHO=WPF (HIMT ) 

IF ( I0AMP.EQ.6) BY3Y2H=BEX (ARGO) 

IF(IOAMP.EQ.7) 6Y3Y2ri=BPF (ARGO) 

C NOW OEFINE OIFFERENTIAL EoUATIONS 
190 IF(M0DEL-2) 200.300*400 

200 DVEC(1)=(K2/M1)*(VEC(5)-VEC(4) ) - (Kl/Ml ) * (VEC ( 4 ) -YB) -Cl/Ml* (OVEC 

1 (4)-YB00T) 

300 OVEC (2) =- ( 1 ./M2) *WY3Y2h* (1 ./M2) *WHO^ ( 1 ,/M2) * (OVEC ( 6) -OVEC (5) ) 

1 *BY3Y2H-(K2/M2)*(VEC(5)-VEC(4) ) 

OVEC(3)=-(l./M3)*WhO* ( 1 ./M3) *WY3Y2H- ( 1 ./M3) * (OVEC (6) -OVEC (5) ) *BY3Y 
12H 

IF(M00EL.EQ.2) OvtC ( 1 ) = (K2/M1 ) <» (VEC (5) -VEC (4) ) - (K 1/M 1 ) * ( VEC (4) - 
1 YB) 

GO TO 401 

400 0VEC(2)=VEC(1) 

OVEC ( 1 ) =-WhO+WY3Y2h- (0VEC(2)-YB00T)*BY3Y2H 
DVEC(1)=DVEC(1)/M3 

401 return 
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100 

FOh^maT (/ { 3A ^ 3MvECf 1 J ) f 1 hT ♦ 6 ( 3X • 4HD VtC ♦ 13) ) 




101 

FOkMA r (lA<t)L9*3»Fb*4*nE10*H) 




10? 

format ( l^'0,?9HMAAlhU(-t NUMttR OF 1 TEKA T IONS* 13# 25H 
iRr n .4) 

WAS 

REACHED BEFO 


1U3 

FuRMaT ( 5X#6hF0R Yti) 




104 

FORMAT (bA,9rSF0K YhOOl ) 
EM[) 





SUhRUUI inf TLU{A*h*C*U*N) 


TLU 

1 


LINEAR INTERPOLATION ROUlINE 


TLU 

2 


A= iNOLRhNDENT \/ARlAbLt 


TLU 

0003 


S= OEPlNLENT VARIABLE (ANSWER) 


TLU 

0004 


C= INDtPENOENF TABLE 


TLU 

000b 


U= OEPENUENT TABLE 


TLU 

0006 


N= number uf points in Table 


TLU 

0007 


independent table must be sorteu* either ascending 

OK 

descending TLU 

0008 


dimension C ( 1 ) »D (1 ) 


TLU 

0009 


IF (N-1 ) 1 *2* 3 


TLU 

0010 

i 

B = 0. 


TLU 

0011 


GO TO 100 


TLU 

0012 

? 

B=D ( 1 ) 


TLU 

0013 


GO TO 100 


TLU 

0014 

3 

ML = 1 


TLU 

001b 


mu=n 


TLU 

0016 

8 

if (MU-ML- 1 ) ISf lb»9 


TLU 

0017 

9 

M= (MU+ML) /? 


TLU 

0018 


IF (C( 1) -C (2) ) 11*2*10 


TLU 

0019 

10 

IF(C(M)-A) 13*12*14 


TLU 

0020 

11 

IF (A-C (M) ) 13* 12* 14 


TLU 

0021 

12 

B=D (M) 


TLU 

0022 


GO TO 100 


TLU 

0023 

13 

MU = M 


TLU 

0024 


GO TO e 


FLU 

0025 

14 

ML = M 


TLU 

0026 


GO TO 6 


TLU 

0027 

15 

B = D(ML)M0(MU)-D(ML) )^{ (A-C(ML) )/(C(MU)-C(ML) ) ) 


TLU 

0020 

100 

RETURN 


TLU 

0029 


END 


TLU 

0030 


FUNCTION FX(X*FK) 

COMMON TPRiNTtTMAXt IPSIEP*NCYCUS»M0T10 N*TITlE (IS) f YBA, YblNT (^00 ) 
Common nr*nw»ncase* mooel *loaof » ioamp* ifoft * icc» 

1 Ml ,M2*M3«K1 »Ci *K2»GL»FA,ALOAO»AOAMP.BLOAO.BOAMPt 

2CL0AD»C0AMP,0L0AD*D0AMP*HINIT* IFM7»AMPLO»FRE(J»PEAKS*PL)UKS*PEAKHt 
3PDURH,NINCP,TINCP.TP(200) »PUVAL (200) »PULDOT (200) »TINIT*Y1IN»Y2IN, 
AY3IN,VlIN,V2IN,V3IN.TIME(6000)»Yl,Y2*Y3»VlfV2»V3.N»VEC(6) ,UVtC (6) » 
5MAXINC»H»HM,M,E,AK0*NE»T* ELIM»INDX* YBtYBUOTt 

60MEGA«WH0* WY3Y2H»BY3Y2h 
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COMMON FILMT(6000) *Y2MY1(6000) * VAPEAK (6000 )♦ YBPLOT ( 6000 ) 

REAL M1*M2,M3*C1*K1,K2 

K=FK 

NP1=NINCP+1 
GO TO (1*2) *K 

1 CALL TLU(X»F1»TP.Y8INT »NP1) 

FX=F1 

RETURN 

2 CALL TLU(X,F2»TP»PUVAL*NP1) 

FX=F2 

RETURN 

END 


SUBROUTINE FORMT 
DIMENSION A(2) 

COMMON TPRINT»TMAX» IPSTEP»NCYCLS*M0T ION* TITLE (15) * YBA*YBINT (200) 
COMMON NR*NW*NCASE*MODEL*LOADF» IDAMP* IFOFT* ICC* 

1 Ml ♦M2»M3*K1,C1*K2*GL*FA*AL0A0*ADAMP»BL0AD*BDAMP* 

2CLOAD*CDAMP*DLOAD*DDAMP,hINlT* IFM7*AMPLD*FREQ* PEAKS* PDURS*PEAKH* 
3PDURH,NINCP*TINCP*TP (200) »PUVAL (200) *PULOOT (200 ) * T INI T ♦ Y1 IN* Y2IN , 
4Y3IN*V1 IN*v2IN*V3IN*TIME(6000) * Y 1 * Y2* Y3 » V 1 * V2 ♦ V3*N* vEC (6 ) *DVEC(6) * 
5MAXINC*H*HM*M,E* AKO*NE*T* ELIM*INDX. YB*YBDOT* 

60MEGA*WHO*WY3Y2H*BY3Y2h 

COMMON FILMT (6000) *Y2MY1 (6000) * VAPEAK (6000) * YBPLOT (6000) 

DATA MAXIT*ERR/100*.01/ 

EXTERNAL FX 

REAL Ml *M2*M3*C1 *K1 *K2 
NP1=NINCP +1 
YBINT(1)=0.0 
DO 10 I=2*NP1 
A ( 1 ) =TP ( I-l ) 

A(2)=TP(I) 

MRET=5 

CALL SIMINI {TEMP*FX*A*ERR*0»MAXIT*2,*MRET) 

IF (MRET.EQ.5) GO TO 5 
GO TO 10 

5 WRITE(Nw*102) MAXIT*ERR 
WRITE (NW*101 ) 

10 YBINT(I)=YBINT(I-l)+TtMP 
RETURN 

101 FORMAT (5X*49HF0« FORMATION OF NEW TABLE FOR THE INTEGRAL OF YB) 

102 FORMAT ( IhO »29HMAXIMUM NUMBER OF ITERATIONS* I3*25H WAS REACHED BEFO 
IRE EPS =*E11.4) 

END 
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SUBROUTINE FILM 

C THIS SUBROUTINE CALCULATES THE EQUILIBRIUM GAS FILM THICKNESS. 

COMMON TFRINT»TMAX» IPSTEP»NCyCLS» mot ION. title ( lb) »YBA. YBINT (200) 
COMMON NR, Nw » NCAS£» MODEL. LOADF » IDAMP.IFOFT* ICC* 

1 Ml ♦M2»M3,K1 ,C1*K2»GL*FA*ALOAO*ADAMP*BLOAD»8DAMP» 

2CL0AD,CDAMP*DL0AD,0DAMP*HINIT*IFM7»AMPLD*FREQ*PEAKS*PDURS»PEAKH* 
3PDURH,NINCP,TINCP»TP (200) »PUVAL(200) »PULDOT (200) » T INIT ♦ Y 1 IN* Y2IN* 
AY3IN,V1 IN* V2IN*V3IN*TIME (6000) ♦ Y 1 * Y2 * Y3 * V 1 ♦ V2 * V3 * N* VEC ( 6) *DvEC(6) * 
5 MAxINC,H,HM,M,E,AK0,NE,T. ELIM.INDX* Y8*YBD0T* 

60MEGA*WM0*WY3Y2H,BY3Y2h 

COMMON FILMT (60 00) *Y2MY1 (6000) * V APEAK (6000 ) * YBPLOT (6000 ) 

REAL M1*M2*M3*C1*K1,K2 
W=“FA+M3*GL 
IF(LOADF-4) 11,11*21 
11 BH = ALOG(i*)-ALOG(ALOAD) 

B=BLOAD 

IF (ABS(BLOAO) . LT .. 0000 1 . AND .BLOAD . LT . 0 . 0 ) 8=-. 00001 
IF (ABS(BLOAD) .LT .. 0000 1 . AND . BLOAD . GT . 0 . 0 ) B=. 00001 

HINIT=BH/B 
RETURN 

21 CABIN =W/CLOAD 
D=DLOAD 

IF (ABS (DLOAD) . LT . . 0 00 0 1 . AND . DLOAD .GT . 0 . 0 ) D=. 00001 

IF (ABS (DLOAD) .L T .. 0 0 00 1 . AND . DLOAD .LT . 0 . 0 ) 0=-. 00001 

ONeNT=ALOG (CABIN) /D 

IF (ABS (ONENT) .GT.aa.O) GO TO 30 

HINIT=EXP (OnENT) 

RETURN 

30 WRITE (NW* 100) DLOAD 
CALL EXIT 

100 FORMAT (25HODLOAD IS TOO BIG. DL0AD=*E15.5) 

END 


subroutine CALCOM 
DIMENSION CLX (20) *CLY (20) 

COMMON TPRINT *TmAX, IHs7EP*NCYCLS*M0TI0N* title (15) ♦ YBA , YB INT ( 200 ) 
Common nr*nw*ncase*model»loadf*idamp*ifoft*icc* 

1 Ml *M2*M3*K1*C1 *K2*GL*FA*AL0AD*ADAMP*BL0AD*BDAMP* 

2CL0AD*CDAMP,DL0AD*DDAMP,HIN1T* IFM7*AMPLD*FREQ*PEAKS*PDURS*PEAKH* 
3PDURH,NINCP*TINCP* TP (200) *PUVAL(200) *PuLDOT (200)*TINIT*Y1IN*Y2IN* 
4Y3IN,V1 IN*v2IN,V3IN,TIME(6000)*Y1*Y2*Y3*V1*V2*V3*N*vEC(6) ,DVEC(6) * 
5MAXINC*H*HM*M,E* AKU*NE*T* ELIM*INDX* YB*Y8D0T* 

60MEGA,Wh0*WY3Y2H*BY3Y2h 

COMMON FILMT (6000) *Y2MY I (6000) * V APEAK ( 6000 ) * YBPLOT (600 0) 

COMMON IBUF(IOOO) 

DATA XORIG* Y0RIG*AXLEN,AYLEN/2. , 1 .5* 12. *6./ 

WRITE (Nw,202) 

202 format ( I hl .30HROAU MARKERS FOR CALCOMP PLOTS) 

NVALS=INDX 
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NPTS=NVALS/IPSTEP 

I=NPTS*IPSTEP+1 

J=r+IPSTEP 

WRITE (NW*200) NVALS»NPTS»I»J 
CALL FACT0R(.75) 

C CHANGES IN CALCOM. INSERT AFTER CALL TO PLOTS 

C STOMP OUT NEGATIVE FILM THICKNESSES. TAG FIRST VAL FOR ERROR MESSAGE, 
ITAG=0 

CALL PLOT (X0RIG,Y0RIG*-3) 

DO 20 K=1,NPTS 
IF(ITAG) 1S»15.16 

15 IF(FILMT(K) ,LT. 0.) 1TAG=K 

16 IF(FILMT(K) ,LT. 0.) F1LMT(K)=0. 

20 CONTINUE 

CALL SCALE (TIME ♦AxLEN.NPTStlPSTEP) 

WRITE(Nw»201) TIME ( I ) ♦TIME (J) 

IF(lTAG) 22f2?.21 
21 XPAGE=TIME(ITAG)/TIME(J) 

CALL SYMBOL (XPAGEfO. ♦ . 14t7hC0NTACT ♦30. ♦?) 

C SCALE the axis ACCORDING TO ThE OATa BUT FORCE THE MINIMUM FILMT TO 0. 
C DEFINE AN ADDITIONAL FILMT TO INSURE AT LEASE ONE ZERO 
22 FILMT(I)=0. 

C CALL SCALE USING I POINTS^ IT WILL STORE 0. AS FILM(J) AND OELTAV IN 
C FILMT (J+IPSTEP) 

NP1=NPTS*1 

CALL SCALE (FILMT ♦A YLEN^NPl -IPSTEP) 

I I=NP1*IPSTEP*1 
JJ= II+IPSTEP 

WHITE (NW^ 203) FILMT (II) ♦FILMT (JJ) 

C NOW MOVE FIRSTV AND OELTAV BACK WHERE THEY BELONG. 

FILMT ( I ) =FILMT ( j) 

FILMT ( J)=FILMT( J+IPSTEP) 

C SCALE OATA FOR TIME AND INPUT MOTION 

CALL SCALE ( T I ME ♦ AXLEN ♦ NPTS ♦ I PS TEP ) 

CALL SCALE (YBPLOT ♦ A YLEN^ NPTS ♦ IPSTEP) 

CALL PLOT (-XORIG^-YORIG^-3) 

WRITE (NW^ 20 A) YbPLOT ( I ) ♦ YBPLOT ( J) 

C DRAW THE AXIS FOR THE APPROPRIATE INPUT MOTION 

0lF(MOTION .EQ. o *AND. IFOFT ,EQ. 8)CALL AXIS ( i . ♦ i .Sf2AHlNPUT VIBR 
lATION - INCHES^24^AYLEN^90.^YBPLOT(I) ♦YBPLOT(J) ) 

OIF(MOTION ,£Q. 1 .AND. IFOFT .EQ. 8)CALL AXIS ( 1 , ♦ 1 .5*37HINPUT VIBR 
lATION - ACCELERATION IN G9‘St37 ♦AYLENtSO. ♦ YBPLOT ( I )♦ YBPLOT (J) ) 
OIF(MOTION .EQ, 0 .ANO. IFOFT ,GT. 8)CALL AXIS ( 1 . ♦ 1 ,5*20HINPUT SHOC 
IK - INChES^20^AYLEN^90. ♦YBPLOT (1) ♦YBPLOT (J) ) 

OIF(MOTION .EQ. 1 .AND. IFOFT .GT, 8)CALL AXIS ( 1 , ♦ I ,5* 33HINPUT SMOG 
IK - ACCELERATION IN GXS^33fAYLEN^90. ♦ YBPLOT ( 1 )♦ YBPLOT (J) ) 

WRITE (Nw^205) 

C DRAW THE AXIS FOR FILM THICKNESS 

OCALL AXIS{X0RIG^Y0RIGt31H6EARING FILM THICKNESS - INCHES^ 
131^AYLEN^90.^FILMT (I) .FILMT ( J) ) 

WRITE (NW.206) 

C DRAW time axis 

OCALL AXIS (XORIG.YORIG. 14HTIME - SECONOS»-14»AXLEN.O. » 
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ITiMtU) ,TIME(J) ) 

CALL PLUT (XORIG. Y0RIG»-3) 
Wk I T E ( Nw « EO 7 ) 

C PLOT THE graphs 


C 


IK (MPTS.LT. JO) NOPON=l 

IFCNHTS.GE.JO) NOPON=NPTS/30 

CALL LINE (1 IM£,YdPLOT» NPTS» IPSTEP»N0P0N»3) 

CALL L ink: ( TIME ,F 1LMT» iNPT S ♦ I PSTEP , NOPON » 1 ) 

WRITE (NWtPOd) 


OPaw a centerline to REPRESENT 
sa call plot (-xopig»-yorig»-3 

WWI IE (Nw,209) 


equilibrium film 


thickness 


DO 99 K=1,1J 


CLY (N) =Y0RIG+HINI I/F ILMT ( J) 

99 CLX (K) =xCRIG+FlOAT (K- 1 ) /12.*AXLEN 
CLX ( lA) =0.0 
CLX(lS)=l.o 
CLY(14) =0. 

CL Y (IS) =1 . 


CALL CNTRLN(CLX»CLY» 1 J» 1 ) 
WRITE (Nw»210) 


RETURN REFERENCE TO LOWER LEFT HAND CORNER AND LETTER TITLE 


100 call symbol (2.b0»0.50« .14* rnLE»0,» 75) 

fncase=ncase 

call symbol ( 1 1 .5* .30 » . 100 ♦ 1 IhCASE NO, = 
call number (999. *0.30* . 10, FnCASE* 0. *-l ) 
call SOATE(OATE) 

call SYMbOL{11.5*0.15,.10*OATE*0.*8) 

PERFORM Annotation 

OIF (MODEL ,EQ. 1) CALL S YMBOL (4 . 50 , 9 , 65 , , 
VERTICAL ON ISOLATORS 

2) CALL SYMBOL (4.50,9.65, , 
vertical on rigid MOUNTS 

3) CALL SYMBOL (4.50,9.65, , 14, 
VERTICAL ON RIGIU MOUNTS, CASING, 

CALL S YMdOL (4.d0 , 9.25, . 14, 1 9HINPUT FUNCTION IS 
OIEdFOFT .EQ. d) CALL SYMBOL (999. ,9.25, . 14, 


1 S3HSIMULATOR 
0IE(MODt.L ,EO. 
1 S3HSIMULATOR 
OIF (MODEL .EO. 
1 SSRSIMULATOR 


* 0 . , 1 1 ) 


.14. 


CD 

CD 


.14, 


,0.,53) 

,0.,53) 

AND FLEXURE, 0. ,55) 
,0.,19) 


1 40MSTE AOY-STATE VIBRATION 
IFdFOFT.EO. 9 . and. MOTION. EQ.O) CALL 
l40HhALF SINE DISPLACEMENT SHOCK 
IF ( IFOF I .F D, 9. AND, MOT ION. EQ, I ) CALL 
lAOHt-ALF SINE ACCELERATION SHOCK 
IF ( I F OF ( .E(J. 1 0 . AND.MO r lON.EQ. 0 ) CALL 
l4i)HHAvERSlNE displacement SHOCK 
I F ( I F uF f . F 0 . 1 0 . AND .MU T lON.EQ . 1 ) CALL 
l40HHAVtKb I\E ACCELERATION SHOCK 
I F ( IF OF T .F 9 . 1 ) . ANo. MU | I UN.EO. 0 ) CALL 
l9iHSHOC< 'HSRLACt.MtN] VALUES REAU IN 
IF ( IFOF 7 .FO. I 1 . AND.MO r lON.EO, 1 ) CALL 
141HSHUCK ACCELr.RAdON VALUES READ IN 
GO TO (1 Ul .120 d 30, 1 AU ) , IFM7 


,0. ,40) 

SYMB0L(999.,9.25,.14, 

,0.,40) 

SYMBOL (999., 9. 25, .14, 
, 0 . , 40 ) 

SYMBOL (999., 9.25,. 14, 
,0.,40) 

SYMBOL(999.,9.25, .14, 
,0. ,40) 

SYMBOL (999.,9.2b,. 14 , 
As data, 0 ., 41) 
SYMB0L(999.,9.25,.14, 
AS DATA, 0., 41) 
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110 CALL SYMBOL(4.5O*8.05*.14»13HAMPLITUDE IS ♦0.*13) 

CALL NOMBEP(999.*8.8b*.14*AMpLD.O.*5) 

111 IF(MOTION.EQ.O) CALL SYMBOL < 999« » 8. 85» . 14»7H INCHES*0.*7) 
IF(MOTION.EQ.I) CALL SYMBOL (999. *8.85* • 14*7H GitS *0.»7) 

60 TO (112*116*116»140) tIFM7 

112 (iALL SYM80L(999.»8.8b*.14»17H FREQUENCY IS *0.*17) 

CALL NUMeER(999.«8.85«.14tFREQ«0. tS) 

CALL SYMB0L(999,.«.8S*.14»3H HZ*0.»3) 

GO TO 150 

115 CALL SYMBOL(4.SO«8.85t.l4tl4HP£AK VALUE IS •0.«14) 

CALL NUMBER (999. *8.85* .lAtOPEAKtO.tS) 

GO TO 111 

116 CALL SYMBUL(999.*8.85t.l4t24H PULsE UURATION IS »0.»24) 

CALL NUMBER ( 999. « 8.8b «.14»OOOR->0.« 5) 

CALL SYMB0L(999.«8.85*.14*8H SECONDS* 0. » 8) 

60 TO 150 
120 OPEAK=PEAKS 
ODOR=POURS 
GO TO 115 
130 OPEAK=PEAKH 
ODOR=PDURH 
GO TO 115 

140 ODOR=FLOAT (NINCP) *TINCP 

CALL SYMBOL(4.50»8.85».14»25HPEAK VALUE IS SHqWn *0.*2S) 

GO TO 116 

150 CALL SYMBOL(4.57*8.32».14»3»0.*-1) 

CALL SYM80L(5.00*8.25*.14»17H = INPUT FUNCTION, 0 .♦ 1 7 ) 

CALL SYMBOL (4.57, 8. 07, .14,1,0.*-!) 

CALL SYM80L(5.00,8.00*.14,32H = THRUST BEARING FILM THICKNESS, 0. , 
132) 

WRITE (NW,21 1 ) 

CLX(16)=4.50 

CLY(16)=7.82 

CLX(17)=5.2S 

CLY(17)=7.82 

CLX(18)=0. 

CLY(18)=0. 

CLX(19)=1. 

CLY(19)=1. 

CALL CNTRLN(CLX(16) ,CLY (16) ,2,1) 

CALL SYMB0L{5.42,7.7b*.l4,26HE0UILIBRIUM FILM THICKNESS»0. ,26) 
CALL PLOT(-.4,0.,-3) 

CALL DASHPT(0. ,11.3, ,113) 

CALL OASHPT(14. 7,11.3*. 147) 

CALL DAShPT(14.7*0.*.113) 

WRITE (Nw*212) > 

200 FORMAT (/5X,6HNVALS=* I5*6H NPTS=*15*3H I=*I5*3h J=*Ib) ‘ 

201 FORMAT (/5X,8HTIME ( I ) = *E 12.5 * bhT IME ( J) = *E12.5) 

203 FORMAT (/bX* lOHFILMT ( 1 1 ) =,E12.5* lOHFILMT (JJ)=*E12.5) 

204 FORMAT (/SX* lOHYBPLOT ( I ) =*E12.5* lOHYBPLOT (J)=*E12.5) 
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205 FORMAT {/5X*25HDIAG VEKT AXIS-Z DRAWN ) 

206 FORMAT (/bX*25HDIAG VERT- AXlS-2 DRAWN ) 

207 F0RMAT(/bX,25HDIAG hORiZ AXIS DRAWN ) 

208 FORMAT (/bX,25HCDRVES PLOTTED ) 

209 FORMAT (/bX»25H0RIGIN SHIFTED ) 

210 FORMAT (/bX»25HCENTERLINE DRAWN ) 

211 FORMAT (/bX,25HANN0TAT JON ALMOST DONE ) 

212 FORMAT(/bX»2bHANNOTATION DONE ) 

RETURN 

END 


SUBROUTINE RONKU ( X * Y »DY ♦ H.HM,N*NE*M *E * AKO) 

DIMENSION Y(6)»DY(6)«Yl(6)»YD(b)*YP(6*5)*0P(6*5)»A(4)*8(A)*C(4) 
30 IF(NE) 1»1»2 

1 A(l)s0.b 

A(2)= .2928932 
A(3)= 1.7071068 

A (A) si .0/6.0 
B(1)s2.0 
B(2)=1.0 
B(3)=1.0 
3(A)s2.0 
C(l)s0.b 
C(2)=A(2) 

C(3)sA(3) 

C(A)s0.b 

IF (NE) 32»11.2 

11 DO 3 
YP<I»1)=Y(I) 

3 QP(I «l)s0.0 
32 K1=0 

2 X1=X 
HH = h 
K = 1 

12 X50 = XI 
DO 4 Js2.b 
DO 17 1=1, N 

17 YD(I)=YP(I.J-1) 

call DIFEO (N.Xl ♦ YD»DY) 

IF (J-3) 151,lb3flb2 

151 Xl=Xl+0.b*HH 
GO To 153 
lb2 X1=X1+0.5«HH 
153 CONTINUE 

DO 5 1=1, N 

TERM=A ( J-1 ) * (DY ( 1) -B ( J-1) *UP n , J-1) ) 

YP(I,J)=YR(I, J-1 ) +HM*fERM 
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5 QP ( I ♦ J) =QP ( I » J-1 ) +3.0*TEPM-C(J-1)«DY (I) 

4 continue 

XI =X50 

IF (K-2) 6»7»8 

6 DO 61 J=1»N 
Y1 (J)=YP( J»5) 

YP (J*l )=Y ( J) 

61 QP(J»1)=0.0 
K = K+ 1 
HH=H/2.0 
GO TO 12 

7 K=K+1 
X1=X1+HH 

Do 71 j=l»N 
QP ( J* 1 ) =QP ( J»5) 

71 Yp ( J* 1 ) =Yp ( J»b) 

GO TO 12 

8 DO 81 J=1»N 

ri=SURT (Y1 ( J) *<^2*YP ( J»b) <>*2) 

IF(Tl) 18»81»lrf 

Id TEST=AHS ( Yl ( J) -YP ( J»6) ) /T1 
IF (TEST-E) 81»81*y 

81 CONTINUE 
K1=K1+1 
AKO=TFST 

IF (Kl-M) 10,13»1J 

9 K1=0 

IF (ABS (H/2.0)-ABS ( HM ) ) 14 » 1 3 ♦ 1 b 

14 AK0=TEST 
GO TO 10 

15 H=H/2.0 

DO 16 J=1»N 

16 YP ( J» 1 ) =Y ( J) 

GO TO 2 

13 K1=0 
H = H*H 

10 X=X1 +HH 
DO 82 J=1«N 
QP(J»1)=QP(J»5) 

Yp(J»l)=YP(J*5)^(Yp(J»5)-Yl(J))/lb.O 

82 Y ( J) =YP ( 1 ) 

CALL DIFE0(N»X,Y»CY) 

RETURN 

END 



APPENDIX C 


PHOTOGRAPHS OF BEARING SURFACES AND COMPONENTS 


Introduction 


This appendix contains a thorough documenca tion of the condition of the bear- 
ings during various inspections defined by the following sequence of test 
events : 

1. Inspection 1 - Beginning of program before any externally applied shock or 
vibration . 

2. Tests as follows: 

a) axial vibration, rigid mounts, 39,000 rpm 

b) axial vibration, isolation mounts, 39,000 rpm 

c) transverse (Tl) vibration, rigid mounts, 39,000 rpm 

d) transverse (Tl) vibration, isolation mounts, 39,000 rpm 

3. Inspection 2 

4. Tests as follows: 

e) transverse (Tl) vibration, isolation mounts, rotor inert 

f) transverse (Tl) vibration, rigid mounts, rotor inert 

g) transverse (T2) vibration , isolation mounts , rotor inert 

h) transverse (T2) vibration, rigid mounts, rotor inert 

i) axial vibration, isolation mounts, rotor inert 

j) axial vibration, rigid mounts, rotor inert 

5. Inspection 3 

6. Tests as follows: 

k) transverse (T2) vibration, isolation mounts, 39,000 rpm 

l) transverse (T2) vibration, rigid mounts, 39,000 rpm 

7. Inspection 4 

8. All Shock Testing 

9 . Inspection 5 
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At the conclusion of each inspection, except the final one, all bearing surfaces 
were cleaned for reasons which will soon become apparent. The photographs pre- 
sented in this appendix display each significant bearing surface, in turn, in a 
time sequence of events which, more or less, affected the apparent condition of 
that surface. 

Figures 65 through 69 display the apparent condition of the surfaces of the 
thrust plate and runner which form two opposing interfaces with the thrust- 
bearing gas film. 

Figure 65 shows the surfaces at the beginning of the test program. Note, that 
some previous wear is evident around the center of the thrust plate (shown on 
the left). The face of the thrust runner, on the other hand, shows little or 
no indication of wear. 

Figure 66 shows the condition of the thrust-bearing surfaces during Inspection 
2c The radial striations and inner circumferential band shown in part a) on 
the thrust runner are primarily the result of fine debris rather chan wear. 
Verification is presented in part b) which shows the same surfaces after being 
cleaned. Figure 67 shows these surfaces at the beginning of Inspection 3. 

Again the surfaces are dirty but noticeably less so than they were at the 
beginning of Inspection 2. This is because a substantially cleaner air supply 
was used to drive the simulator turbine after Inspection 2. Figure 67 shows 
that some circumferential scratches now appear on the tbi'ust runner which are 
not removable by washing. These were shallow surface markings and did not 
affect bearing performance. Figure 68 again shows a small amount of debris 
accumulation between Inspecticn 3 and Inspection 4. Figure 69 shows these 
surfaces at the conclusion of all the shock and vibration testing. The five 
spots lying in a circle on the thrust runner resulted from suspending the 
runner in a stationary position by feeding pressurized air through the five 
small orifices in the thrust plate. 

Figure 70 shows the condition of the reverse-thrust bearing and opposing runner 
face at the beginning of the program, and Figure 71 presents the same surfaces 
at the conclusion of the testing. The outer circumferential scratch was 
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present to begin with. The difference in the two Figures is attributable to 
debris accumulation between Inspection 4 and Inspection 5. 

Figures 72 and 73 form a bef ore-and-af ter set for the compressor bearing 
journal. Except for occasional debris accumulation, there was no noticeable 
change in the surface condition throughout the program. 

Figure 74 shows the condition of the faces of the compressor journal-bearing 
pads at the beginning of the program. The pad numbering (convention.. is., defined 
with respect to Figure 2 as follows: 

In Section AA (Bearing #1), beginning with the pad containing PF8, 
the pad numbers increase in the direction of rotation shown by the 
curved arrow; 

The pad numbers also increase in this direction in Section BB (Bearing #2) 
starting with the pad containing PF22. 

Figure 75 shows the same bearing surfaces at the conclusion 'of all the shocK.- 
and-vibration testing. Again, the only difference is in the small amount of 
debris accumulation shown in Figure 75* 

Figure 76 shows the outer surfaces of the compressor journal-bearing pads 
along with their corresponding pivots. These photos were taken before 
comniftncement of the testing. Figure 77 shows these same surfaces at the con- 
clusion of all vibration testing. The only significant effect that the vibration 
testing produced was creation of two elongated wear spots on the surfaces of 
the pivot ball and pivot socket. The effect of this v^ar was a slight, but 
noticeable increase in the overall diametral clearance at each journal bearing. 
This did not produce any observable degradation in machine performance, however. 
Although the wear is detectable in Figure 77 , it is shown in better detail in 
the magnified views of Figure 78 . At this point, the pivot balls were 
replaced with new ones and shock testing was begun. At the end of the tests 
the photographs of Figure 79 were taken--showing that the shock testing did 
not result in the same amount of wear as did the vibration testing. 
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Figures 80 and 81 document the condition of the turbine bearing journal 
before and after all the testing. As was the case for the compressor journals 
the turbine journal is virtually ‘unaffected by the testing. 

Figures 82 through 86 are the turbine journal-bearing counterparts of 
Figures 74 through 79 . The turbine journal- pad faces exhibited no signs of 
wear (Figures 82 , 83 ) , while the pivot balls received some wear spots during 
the vibration testing (Figures 84 and 85), were subsequently replaced, and were 
not damaged nor worn by the shock testing (Figure 86). 

In summary, all the primary bearing surfaces showed no deterioration as a re- 
sult of the vibration and shock testing. Considerable debris was picked up by 
the pad faces and thrust-bearing surfaces in the interval between the first 
two inspections. Lesser, but still noticeable debris accumulated between each 
successive inspection interval. The pivot-ball surfaces showed some wear 
resulting from the vibration testing, but their replacements were virtually 
unaffected by the shock tests. Aside from the flexure which fraccured during 
a resonance condition, no damage to the machinery resulted from the test pro- 
gram. The simulator was performing just as well at the end of the test program 
as it was at the beginning. 
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APPENDIX D 


METALLURGICAL REPORT ON FLEXURE FAILURE 


The rotor-bearing system used in the simulator accurately represented the 
turbocompressor rotor-bearing system except for the materials from \^hich 
certain of the parts were fabricated. Because of the long delivery time 
involved in obtaining the specific alloys selected for the turbocompressor 
design, several material substitutions were made for the simulator design. 
These substitutions were permissible since the simulator, unlike the turbo- 
compressor, would be operated at low temperature conditions. One of the 
simulator components fabricated from materials other than the materials 
specified for use in the turbocompressor were the flexures used to support 
the turbine-end journal bearing. The simulator flexures were fabricated 
from AISI 4130 steel instead of Discalloy which was specified for the turbo- 
compressor flexures . 

During the vibration test with the simulator oriented vertically and running 
at 39,000 rpm, the test instrumentation indicated the presence of a resonance 
in one of the turbine-end flexure-supported bearing-pad assemblies. While 
investigating the nature of this resonance, which occurred at a shake table 
frequency of 1365 Hz and an input acceleration level of 7 g's, the journal 
bearing flexure failed. 

The failure was detected from the instrumentation signals which describe the 
position of the journal relative to the casing, and from the signals which 
describe the film thickness between each bearing pad and the journal. At 
the instant of the failure, the signals indicated a sudden displacement of 
the journal and several brief contacts between the pads and journal. The 
simulator was shut down without difficulty. The broken flexure was located 
and removed during a partial disassembly of the simulator. Visual inspection 
of the journal and bearing pads did not reveal any damage to parts other than 


98 



to the flexure; therefore, the simulator was reassembled with a new flexure 
and the test program continued. No further failures of the flexures have been 
experienced during vibration testing. 

The purpose of this appendix is to present a failure analysis of the AISI 4130 
steel flexure. 

Discussion 


The fractured flexure is shown in Figure 87a (side view) and Figure 87b (top 
view) . It is readily apparent that the failure occurred at the radius between 
the beam and shoulder sections of the flexure. This particular section of 
the component is in effect a notch and, as would be expected, is susceptible 
to resultant stress magnification effects. The relevant metallurgical data 
is given in Table III« 

In Figure 88a, the macro-appearance of the fractured surface is given. There 
is a notable absence of fibrosity or shear lip indicating a brittle type of 
failure. Figure 88b is a transverse view of the cracked area in the "as polished" 
condition at 100-power magnification. It is important to note the rather heavy 
oxide layer present in the crack itself. This indicates that considerable 
oxidation occurred before the catastrophic failure. 

Two different areas of the fractured surface are depicted in the electron 
micrographs presented in Figures 89a and 89b, respectively. In both figures, 
the predominance of flat surfaces can be observed; the absence of clear-cut 
dimpling is also apparent. This further substantiates the brittle nature of 
the fracture. 

The grain size associated with the temperature of the flexure at the time of 
quenching (prior austenitic grain size) and the microstructure are shown in 
the optical microphotographs of the failed area shown in Figures 90a and 90b, 

The grain size shown in Figure 90a is extremely fine and is a reflection of 
a low austenitizing temperature. This is a highly desirable grain structure. 
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Table III 


RELEVANT METALLURGICAL DATA 


CHEMISTRY 


c 

Mn 

Si 

Cr 

Mo 

32 

.55 

.29 

.86 

.19 


hardness 

Rc 29-32 
CLEANLINESS 


Satisfactory. Primary inclusions are manganese -■ sulfides and lesser 
amounts of silicates. 

GRAIN SIZE 


Extremely fine, smaller than ASTM-9 


MICROSTRUCTURE 


Duplex - Predominately tempered lower bainite and tempered martensite. 
FRACTURE APPEARANCE 


Brittle 
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particularly from toughness considerations. However, the microstructure shown 
in Figure 90b is duplex in nature, consisting of a mixture of bainite, 
martensite, and carbide. The presence of these micros tructural aggregates is 
a manifestation of the shallow hardenability characteristics of AISI 4130 
steel. Duplex microstructures invariably exhibit lower toughness and greater 
notch sensitivity than either lower hainiie or temper edumartensite. The 
physical properties of toughness and notch ductility at a given strength 
level collectively play an important role in determining the fatigue behavior 
of a given material. The microstructure which develops the optimum level 
of these physical properties must be developed, if possible, during heat 
treatment , 

Metallurgical Observations 

1. The metallurgical cleanliness of the AISI 4130 material is 
satisfactory and representative of good melting and deoxidation 
practice . 

2. The grain size is acceptable. 

3. The microstructure is duplex in nature consisting of bainite, 
martensite, and carbide. This is a highly undesirable micro- 
structure due primarily to its low toughness, its higher transition 
temperature, and its notch sensitivity. 

4. The failure was brittle in nature with very little evidence 
of plastic deformation present. 

5. Considerable oxidation occurred in the crack area before failure. 
Conclusions 


During the quenching of low alloy steels from their austenitizing temperature 
differences in thermal contraction and expansion occur as a result of the aus- 
tenite transformation. These conditions result in large macroscopic 
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(measureable) stresses which must be relieved by plastic yielding. Upon 
occasion, these macroscopic stresses exceed the fracture stress of the material, 
particularly when stress raisers are present, such as the fillet radius at the 
ends of the beam section of the flexure. Where relatively large transitions in 
material volume occur in a given part, the susceptibility to quench cracking 
is greatly increased. 

The most probable cause of the flexure failure was the prop^aga^tioniof a oracdc at 
the abrupt change in section between the flexure beam section and end shoulder. 
Initial formation of the crack is believed to have occurred during fabrication 
of the flexure; specifically, during the heat-treat quenching operation. 
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MATHEMATICAL SYMBOLS 


Arbitrary constant used to obtain curve fit; acceleration amplitude, g ’ s 

Coefficient of viscous damping of thrust-bearing gas film, a nonlinear 
function of arguments indicated in suffixed parentheses, Ib-sec/in; 
arbitrary constant used to obtain curve fit. 

Arbitrary constant used to obtain curve fit; acceleration amplitude, g*s 
Coefficient of linear viscous damping, Ib-sec/in. 

Arbitrary constant used to obtain curve fit; acceleration amplitude, g's 
Base of natural logarithms. 

Net aerodynamic force acting on rotor, lb. 

Net force exerted by thrust-bearing gas film, a nonlinear function of 
arguments usually indicated in suffixed parentheses, lb. 

2 

Local gravitational acceleration, in/sec o 

Thrust-bearing gas-film thickness, in. 

Journal-bearing gas-film stiffness, Ib/in. 

Linear spring stiffness, Ib/in. 

2 

Mass, lb- sec /in. 

Time parameter associated with shock pulse, sec. 

Time, sec. 

Velocity, in/ sec . 

Static component of thrust-bearing gas-film force, a nonlinear function 
of arguments usually listed in suffixed parentheses, lb. 

Displacement with respect to inertial reference frame, in. 

Displacement relative to position of static equilibrium, in. 



Subscripts 


1, 2, 3 = first, second, or third mass elements respectively, 

b = Base 

Superscript Conventions 

= Condition at static equilibrium 
d 

* dt 
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Fig. 1 Schematic of Turbo compress or Simulator Mounted 
On Test Stand (Vertical Orientation) 
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Figure 2 Schematic of Turbocompressor Simulator With 
Capacitance Probe Designation (Horizontal 
Orientation) 











CAPACITANCE 


PROBE 


/ 

HOUSING 


ORFICES FOR 
HYDROSTATIC 
OPERATION 


THRUST PLATE 


Fig. 6 Photograph of the Simulator Forward (Primary) 
Thrust Bearing 



Fig- 7 Arrangement of Journal Bearing Pad 








Vibration Table With Simulator Mounted Horizontally 
•;in Shock and Vibration Isolation Isolators 
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EXCITATION FREQUENCY, HZ 


Fig. 11 Maximum Vibration Test Objectives 
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ACCELERATION, 


SHOCK MACHINE PULSE 



Fig. 13 Comparison Between Actual Shock Pulse and Two Idealized Shock Pulses 
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Fig. 15 Three-Degree'-of-Freedom Axial Response Model of 




imulator on Shock and Vibration Isolators 
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ROTOR 
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GAS FILM 

THRUST PLATE 
G(h,h') and effective 
FLEXURE MASS 
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STIFFNESS 
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HOUSING 
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HOUSING 

STIFFNESS 
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APPLIED 
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T f 
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fx 
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Fig. 16 Three-Degree-of-Freedom Axial Response Model of Simulator on Rigid Mounts 
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MODEL 
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STATIC LOAD, W (LB) 



Fig. 18 Static Load Characteristics of Thrust Bearing 
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FILM DAMPING COEFFICIENT, B(LB-SEC/IN.) 


SPEED: 39,000 RPM 

BEARING AMBIENT: AIR AT I20"F AND 14.7 PSIA 
note: dashed REGIONS OF SOLID CURVES ARE 
ESTIMATED. (DATA POINTS WERE NOT 
CALCULATED IN THESE REGIONS.) 


■FUNCTION USED FOR 
SHOCK AND VIBRATION RESPONSE 
CALCULATIONS B = 9.73h''-^ 


MEAN FILM THICKNESS, h(MILS) 

19 Damping Coefficient Characteristics of Thrust Bearing 




INFERRED FROM MANUFACTURER'S CATALOG 



Fig. 20 Comparison Between Theoretical Transmissibility and Catalog Data for a Single Isolator 
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Fig. 21 Measured Static Load Characteristic of a Single Test Isolator 
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(a) ISOUTOR COMPONENTS 




(b) CROSS SECTKM (c) ASSEMBLED ISOLATOR 


Fig. 22 Construction of Shock and Vibration Isolator 



ROTOR 

THRUST PLATE 
HOUSING 


REGION 

REPRESENTED 
BY THIN 

CIRCULAR PLATE 


(A) 

(B) 

(0 



LOADING 


LOADING . 

u - i 

LOADING 


K = 



32.000 LB/IN 

210.000 lb/in 

340.000 lb/in 


Fig. 23 Representations of Thrust-Plate Housing for 
Stiffness Calculations 
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LOAD, LB. 



Fig. 24 Measured Stiffness of Thrust-Plate Housing 
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INCH 











Figure 29. 


Test Setup to Measure Phase Shift 
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.25 MIL/MAJOR DIV. 



Fig. 31 Measured Modal Behavior at 870 Hz 
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SHAKE TABLE FREQUENCY -HZ 


* 

Fig. 33 Maximum Axial Vibration Levels for Safe Operation 
of the Simulator When Rigidly Mounted on Vibration 
Table (Shaft Rotating at 39,000 rpm) 




PEAK TRUNNION ACCELERATION - g s 




AXIAL VIBRATION SIMULATOR 
MOUNTED ON 
VIBRATION ISOLATORS 
ROTOR SPEED - 39,000 RPM 



riTEST objective: 


-MAXIMUM TEST LEVELS 
ACHIEVED FOR 
SAFE’BEARING OPERATION 


ARBITRARILY DEFINED AS 
0.0001 INCHES MINIMUM 
BEARING FILM CLEARANCE 


SHAKE TABLE FREQUENCY - HZ 


Fig. 34 Maximum Axial Vibration Levels for Safe Operation 
of the Simulator When Mounted on Shock and 
Vibration Isolators (Shaft Rotating at 39,000 rpm) 



RATIO OF THRUST- BEARING FILM EXCURSION TO SHAKER AMPLITUDE 


3.0 


2.5 



1.5 


1.0 


0.5 



— CALCULATED BY 3 DE6REE-0F- FREEDOM MODEL 
“--CALCULATED BY SINGLE DEGREE- OF -FREEDOM MODEL 



• MEASURED VALUES 

I I 
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140 


160 


EXCITATION FREQUENCY -HZ 


Fig. 35 Calculated and Measured Thrust Bearing Film Response 
to Axial Vibration About Region of First Resonance 
(Simulator Rigidly Mounted) 











TIME FROM TABLE RELEASE - MILLISECONDS 


Fig. 37 Calculated and Measured Thrust Bearing Film Response 
to Axial Shock (Simulator Rigidly Mounted, Single- 
Degree-of-Freedom Analysis Model) 
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HERRING FILn THICKNEpo - l.NCMEb __eERR I NG FI Lf1 .THICKNESS - I.NCHES -10 


SIMULATOR VERTICAL ON RIGID MOUNTS. CASING. AND FLEXURE 
INPUT FUNCTION IS HAVERSINE ACCELERATION SHOCK 

PEAK VALUE IS 15*70000 G’S PULSE DURATION IS 0*01460 SECONDS- 


+ = INPUT FUNCTION 

O = THRUST BEARING FILM THICKNESS 

EQUILIBRIUM FILM THICKNESS 



(A) UNDAMPED GAS FILM 


SIMULATOR VERTICAL ON RIGID MOUNTS. CASING. AND FLEXURE 
INPUT FUNCTION IS HAVERSINE ACCELERATION SHOCK 

PEAK VALUE IS 15.70000 G’S PULSE DURATION IS 0.0146Q SECONDS 

+ = INPUT function 

0 = THRUST BEARING FILM THICKNESS 

EaUILlSRIUM FILM THICKNESS 





(B) DAMPED GAS FILM 


Fig, 38 Calculated Response of Thrust Bearing Gas Film to a 
Haversine Approximation of the Axial Shock Pulse to 
Illustrate Effects of Film Damping 
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FILM THICKNESS - MILS 



TIME FROM TABLE RELEASE - MILLISECONDS 


Fig. 39 Calculated and Measured Shock Response During Collapse of 
The Thrust-Bearing Gas Film 


C>J 



DAMPING RATIO (C/2/KH) 



.2 ,4 ,6 .8 1.0 1.2 1.4 1.6 1.8 

FILM THICKNESS, MILS 


Fig. 40 Damping Ratio of the Thrust-Bearing Gas Film Versus 
Film Thickness 







BASE ACCELERATION -6 UNITS FILM THICKNESS-MILS 


THRUST RUNNER CONTACTS 



(a) THRUST- BEARING GAS FILM THICKNESS 



Fig. 41 Calculated and Measured Thrust-Bearing Gas Film Response 
To Axial Shock (Simulator Mounted on Isolators, 

Shafting Rotating at 39,000 rpm) 
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CASING ACCELERATION - G s FILM THICKNESS - MILS 


^THRUST RUNNER CONTACTS 
REVERSE SIDE OF BEARING 



(a) THRUST- BEARING GAS FILM THICKNESS 
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TIME FROM TABLE RELEASE - MILLISECONDS 
(b) CASING ACCELERATION 


Fig. 42 Calculated and Measured Thrust-Bearing Gas Film Response 
to Axial Shock (Simulator Mounted on Isolators, 
Calculated Values Based on Acceleration of Simulator 
Casing) 
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*EACH PAD AND EACH IS VARIABLE IN FREQUENCY, 

ISOLATOR HAS INDIVIDUAL PHASE AND AMPLITUDE 

STIFFNESS AND DAMPING 


Fig. 43 Parameters of Computer Program Used to Calculate Transverse 
Vibration Response 
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SHAFT 

AXIS 


RADIAL MOTION 



PIVOT 

AXIS 


Fig. 44 Definitions of Journal Bearing Pad Motions 









NORMALIZED STIFFNESS COEFFICIENTS, 
NORMALIZED DAMPING COEFFICIENTS, B 



(BASED ON ROTOR SPEED OF 39,000 RPM) 


Fig. 47 Frequency Dependence of Normalized Journal-Bearing 
Stiffness and Damping Coefficients 
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ABSOLUTE TRANSMISSIBILITY 



1500 LB/IN 
140/32 SLUGS 

MEASURED VALUES OF Xrnax^^max 


C= 15 LB-SEC/IN 

C = 20 

C = 25 

-C=30 

-C = 40 

:C = 50 






















FORCING FREQUENCY - HZ 


Figure 48 Method of Determining Isolator System Parameters 
Showing Effect of Damping 
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DIRECTION OF 
APPLIED 
VIBRATION 
(VERTICAL) 


Fig. 49 Measured Orbits of Rotor at 166 Hz During 
Transverse Vibration Testing 
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DIRECTION OF 
APPLIED VIBRATION 
(VERTICAU 


Fig. 50 Measured Orbits of Rotor at 225 Hz During 
Transverse Vibration Testing 









Fig. 51 Calculated Mode Shapes for the Two Rigid-Body 
Critical Speeds of the Rotor-Bearing System 


155 



ROTOR ORBIT AT COMPRESSOR- END JOURNAL 


ROTOR ORBIT AT TURBINE- END JOURNAL 


APPLIED VIBRATION LEVEL =12 g‘s 
FREQUENCY = 166 HZ 



APPLIED VIBRATION LEVEL 
FREQUENCY = 166 HZ 


: 12 g's 


PEAK VALUE 2.77 MILS-Q 
IN DIRECTION OF S 
APPLIED VIBRATION JSi 



MEASURED 

CALCULATED 


□ 0.1 MIL 
0.1 MIL 


DIRECTION OF 

rotation 


DIRECTION OF 
applied VIBRATION 
(VERTICAL) 


I 

I 

I 
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Fig. 52 Comparison Between Measured and Calculated 
Orbits at 166 Hz 
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Fig. 53 Comparison Between Measured and Calculated 
Orbits at 225 Hz 
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Fig. 54 Calculated Transverse Frequency Response Showing 
Region of First Two Rigid-Body Critical Speeds 
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Fig. 56 Effect of Journal Bearing Damping on Calculated 
Transverse Frequency Response 
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Fig. 58 Maximum Transverse Vibration Levels in the T1 Direction for Safe Operation of the 
Simulator When Rigidly Mounted on Vibration Table (Shaft Rotating at 39,000 rpm) 
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Fig. 59 Maximum Transverse Vibration bevels in the T1 Direction for Safe Operation of tlie 
Simulator UTien Mounted on Vibration Isolators (Siiaft Rotating at 39,000 rpm) 
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Fig. 60 Maximum Transverse Vibration Levels in tlie T2 Direction for Safe Operation of Che 
Simulator IVhen Rigidly Mounted on Vibration Table (Shaft Rotating at 39,000 rpm) 
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Fig. 61 Maximum Transverse Vibration Levels in the T2 Direction for Safe Operation of the 
Simulator l^Then Mounted on Vibration Isolators (Shaft Rotating at 39,000 rpm) 
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Fig. 62 Maximum Transverse Vibration Levels in the T1 Direction for Safe Operation of the 
Simulator With Nonrotating Shaft, and With and Without Vibration Isolators 
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Fig. 63 Maximum Transverse Vibration Levels in the T2 Direction for Safe Operation of the 
Simulator With Nonrotating Shaft, and With and Without Vibration Isolators 
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Fig. 65 Forward Thrust Plate and Runner, Inspection 1, After Cleaning 


(b) After Cleaning 


Fig. 66 Forward Thrust Plate and Runner, Inspection 2 






{a) Before Cleaning 
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(b) After Cleaning 



Fig. 67 Forward Thrust Plate and Runner, Inspection 3 
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Fig, 68 Forward Thrust Plate and Runner, Inspection 4, Before Cleaning 
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PAD 4 


BEARING -I PAD:^3 
(b) 


Fig, 76 Compressor Journal Bearing Pads, Rear View with Pivots, Inspection 1 





BEARING#! PAD# I 


BEARING#! PAD #2 


(a) 



Fig. 77 Compressor Journal Bearing Pads, Rear View with Pivots, Inspection 4 



(b) Ball 


Fig. 78 Compressor Journal Bearing Pad Pivot Surfaces, Inspection 4 
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. Kear View with Pivots Tv , 
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BEARING #2 PAD #4 BEARING #2 PAD #3 

(b) 

Fig. 82 Turbine Journal Bearing Pads, Face View, Inspection 1 
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BEARING #2 PAD #4 BEARING #2 PAD #3 

(b) 

Fig. 83 Turbine Journal Bearing Pads, Face View, Inspection 2 
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(b) 


Fig. 86 Turbine Journal Bearing Pads, Rear View with Pivots, Inspection 5 







CaJ Macro View - 2X 



(b) Transverse View - lOOX 

Fig. 88 Magnified Views of Fractured Flexure Surfaces 
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(a) Prior Austenitic Grain Size - Etchant 2% Picral - lOOOX 



(b) Typical Microstructure - Etchant 2 % Nital - lOOOX 
Fig. 90 Grain Structure Micrographs 
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